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CONTROL OF A FIRST ORDER DELAYED SYSTEM BY MEANS OF 
AN ASTATIC REGULATOR AND NONLINEAR CORRECTION * 


S. V. Emel*yanov 


(Moscow) 


A technique employing nonlinear converters which alter the structure of de- 
layed dynamic systems during the transient process by means of two coordinates 
(controlled parameter and its derivative) is proposed as a method of stabilizing 


and improving delayed dynamic systems described by second-order differential 
equations. 


INTRODUCTION 


Often it is possible for the purpose of a concrete technical problem to approximate a motion equation of 
a high order by an equation of the first or second order with a delay [1-3]. 


It is known that existing automatic control systems which contain delay elements possess a small critical 
gain, Attempts to increase the critical gain by introducing a reaction “om the derivative of the controlled co- 
ordinate usually do not produce the required results [4]. The present article examines the possibility of applying 
the method of correction described by the author in [5, 8] for the stabilization and improvement of the quality of 
transient processes in automatic control delayed systems, represented by a second order differential equation, In 
the article a technique of determining the required linear control law , which would provide a high quality tran- 
sient process, is suggested. 


A. A. Bashkirov's analytico-graphical method [6] was used as a means for calculating the parameters of the 
correcting device, 


Movement Equation 


The movement equations of the above system can be written in the form: the equation of the object 


Tat + px = pe + f(t), 
pee 


(with f (t) = 0, f*(t) =# 0, the delay pertains to the object; with f(t) #0, f*(t) = 0, the delay pertains to the 
regulator); the equation of the controlling device: 


(1) 


g(t)—z=¢9; (2) 
the equation of the correcting device: 
E=9 (9.9) (3) 


* The work was carried out under the direction of B, N. Petrov, 


the equation of the servomotor: 
(4) 


Here x 1s the relative deviation of the controlled coordinate, ¢ is the relative variation of the error signal, 
€ is the relative deviation of the output variable of the correcting device, * is the relative deviation of the 
controlling element,T,, p, Ts are the constant coefficients which characterize the dynamic parameters of the 
system, r is the delay constant,@(¢, %} are the nonlinear gains of the correcting device, f(t) and f*(t) are ex- 
ternal disturbances and g(t) is the master effect. 


The functional schematic of the control system is shown in Fig. 1. 


Determining the Structure of a Correcting Device in a Delayed ACS 


In closed automatic control systems with delay (Fig. 1) in the presence of disturbances f(t), f*(t) and g(t) 
it is impossible to reduce transient processes to zero at any values of the closed system gain or any (linear and 
nonlinear) control relationships. In such systems the control element compensates the external disturbance after 
a certain interval of time 7, hence even in an ideal control sys- 
z wt) tem, i.e,, without any limitations and with an actuating mechanism 


/ 
pp aa possessing infinite power, deviations of the controlled variable 
from its set value will arise. 


flu 


In article [7] it has been shown that the time of the transi- 
1g 9) ent process in such systems can at the best be equal to r. 


The shapes of the best transient processes in an ideal ACS 
for sudden disturbances g(t), f(t) and f *(t) are shown in Fig. 2. 


a) é 
2 Let us determine the structure of the correcting device 
which would ensure a transient process approaching the ideal inan 
Fig. 1 ACS with a finite gain of the closed circuit and without consider- 
ing the coordinate limitations, Let coefficients T,,p, and T; 
which characterize the dynamic and static properties of the control system be given, It is required to determine 
the control law which would ensure a stable system with a sufficiently high gain in the closed circuit, and a transi- 
ent process not requiring readjustment and providing with a given gain a minimum adjusting time. 


If the regulator is affected only by the error signal (1/5), where 1/6 is the constant gain, the control sys- 
tem will become oscillatory even at a small value of 1/5, moreover, the amplitude of the first half-wave due to 
disturbance f(t) decreases with a rising gain 1/6. 


Let us choose a control law ¢(t) in such a manner that the transient process corresponds approximately to 
the first half-wave of the process, which occurs in the system without a correcting device, The gain of the sys- 
tem without a correcting device is chosen in such a way that the first half of the process satisfies the given tech- 
nical requirements, Then the quality of the transient process in the system with a correcting device will depend 
on the gain 1/5 (with the remaining constant coefficient of the system being given), 


To ensure the stability of the process, and render the control system stable, the controlling element must 
counteract the disturbance at a time r before the end of the transient state, and must remain in this condition 
to the end of the transient state. With certain values of gain the compensation may occur earlier than time 
T before the end of the transient, and at that time the controlling element will deviate by a greater value than 
required for the compensation of the disturbance, Hence at any time from the beginning of compensation to time 
t before the end of the transient it should be possible to return the controlling element to the new stable state 
and hold it there till the end of the transient process, If this occurs at a time r before the end of the transient 
process, the controlling time will be at its minimum (the action of the controlling element will be utilized to 
the full), If this occurs a little earlier (but not before compensation begins), the transient process will be pro- 
longed, yet will not require any readjustments. 


Thus, in the presence of disturbance f(t) the regulator must displace the controlling element in the direction 
of the disturbance compensation according to the control law ¢€ (t) = (1/5) , and then reverse its movement, 
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Fig. 2 


bringing the controlled element to a position corresponding to the balanced condition of the system, Moreover, 
the control law must change by €(t) = —k,¢y, where k, is a constant coefficient, Further the regulator must hold 
the controlling element in this position to the end of the transient state [€(t) = 0}. The instant at which the di- 
rection of the controlling element must be reversed is determined from the condition that the controlling ele- 
ment must have reached its balanced position time r before the end of the transient state in the ACS, 


In other words, the regulating system must change during the transient state, the structure of the closed 
circuit three times: at first the controlled member must have a negative integrating feedback, next, it must 
have positive feedback and, finally, from time rt before the end of the transient state, it must remain open, If 
by the end of the transient state there are no new disturbances, coordinate ¢ will reach its zero value, since quan- 


tity f(t)— M(t) = 0. 


Subsequent disturbing pulses will be compensated in a similar way. It should be noted that disturbances of 
the nature of a single jump can be represented as certain nonzero initial conditions, Hence, if the control sys- 
tem is efficient at any combination of initial conditions, it will be efficient in the case of nonzero initial condi- 
tions, This corresponds to changing the value of the disturbance at the instant when the transient process in the 
system due to the previous disturbance has not yet been completed. 


The efficiency of the system at any initial conditions can be easily determined in the course of subsequent 
explanations, 


In designing the correcting device it is necessary to distinguish two cases: 


1) When the closed circuit gain is such that quantity t* is larger than r (t* is the time interval between 
the maximum deviation of the controlled coordinate and the end of the first half-wave of the transient process 
without any additional stabilization), 


2) When the closed circuit gain is such that t* is smaller than r, Let us synthesize the correcting device 
for the first case, 


Let a disturbance in a controlled system with a control law of the type €(t) = (1/5) g(t) produce a transiént 
process, which can be either stable or unstable (Fig, 3). Let us plot from the end of the first half-wave of the tran- 
sitional process ¢(t) a length equal to the quantity rt. This fixes the instant at which the controlling member must 
occupy a place, corresponding to balanced position of the system, In time interval t-r, where t is a duration of 
of the first half-wave g(t) , let us mark the instant at which the system structure changes, In general it is imma- 
terial at which instant in interval t,, this change of structure will take place; however, the nearer to the end of the 
time interval t-r this occursthe less will the transient process differ from a half-wave, The transient process will 


ch 
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otherwise be extended with respect to the first half-wave to the extent shown in Fig, 3a. 


Thus, the control law € (t) must be the following: in the first interval (t,) it should be €(t) = (1/5) ¢, in 
the second interval (t,), €(t) = —k, ¢, and in the third interval (ts), €(t) = 0 (Fig. 3b), This control law can be 
formulated as follows. 


Let one of the components of € (t) equal (1/5) ¢. At instant t, let us subtract from it component Kj; ¢(t), 
where kj, = ky + 1/6, and at instant t, + t, let us add to it component kj, ¢, where Ki, = ky; 


The second and third components must have in the interval t, and t, + t, respectively zero values of the con- 
stant coefficients and in one jump change these values to quantities kj, and Rig. 


These components can be obtained from the nonlinear functions ¥ described in [5, 8]. Thus, for instance, 
the second component can be described in [5, 8]. Thus, for instance, the second component can be represented 
in the form of a nonlinear function Vai (%, 9) ¢. It follows from [5, 8] that nonlinear function ¥* jaja? #) is 
equal to 


Fis = at gp>0 (5) 
0 at | < J ins 


where J% = kj,/T*j, and K*j, and Tf, are constant coefficients. 


The third component can be represented in the form of function ¥* nik ¢, %)¢. The nonlinear coefficient 
9, 9) is equal to: 


Fis, = f at > 0 
0 at 


< Jj, 


(6) 


where J* jg = T* jp, and,ke j2 and T* are constant coefficients, 


Thus, the control law which provides a transient process without readjustments in a delayed ACS has the 
form (Fig. 3b and c): 


E(t) = (9. 9) = — 9)¢- 


| | 


Fig. 4. 


The circuit of the correcting device which produces this law is shown in Fig, 4, It consists of two correcting 
devices of the “key"type, which produce nonlinear functions ¥, an amplifier with a gain 1/6 and and adder of 
those three quantitites, 


The following are the equations of the correcting device elements: signal amplifier equation 


1= @) 
¥-cell 1 equation 
X= — 9) 9% (9) 
¥-cell 2 equation 
X= Fis, 9) 93 (10) 
adder equation 
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where np is the relative deviation of the amplif‘er output coordinate, x , and x 9 are the relative deviation of the 
¥-cells output coordinates, 


When Tr > t*, the system will not work under tracking conditions, hence it is advisable to be limited by a 
control circuit gain which does not involve quantities t exceeding the value of t*, 


A Technique of Adjusting the Nonlinear Correcting Device Coefficients 


Coefficients T,, p, and T, of the automatic control system depend on its design and dynamic parameters, 
Delay r can also assume various values depending on the physical characteristics of the system, Thus arises 


the problem of determining the correcting device coefficients for given constant coefficients of the system and a 
given delay. 


Bashkirov's [6] analytico-graphical method is taken as the basis for determining the adjustment of the cor- 
recting device, Let us examine the adjustment of the correcting device coefficients for the case when Tr < t*, 


Let the values of the constant coefficients of the system be: T, = 1 sec, T, = 1 sec,p = 1 and r = 0.5 sec, 
The gain 1/6 is chosen for the condition that t* is larger than r (1/6 = 5), Let us plot by means of Bashkirov's 
method the first half-wave of the transient process gt) which has no linear correction and is due to disturbance 
f (t) (Fig. 5). Let us plot from point c (the end of the first half-wave of the transient process) length cp equal to 
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Fig. 5 


quantity r, Ata certain distance from point p along length ep , let us mark the instant of the structural change 
of the system, This instant can be chosen at any point along the length including points e and p. It is desirable 


| 
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to pick a point near to p since in that case the transient process wil be near to the first half-wave of the transient 
process of a linear system. In approaching this instant to p, however, it is necessary to increase coefficient kj, 
and kjz. Hence it is desirable to consider the requirement of improving the quality of the transient process in 
conjunction with the possibility of obtaining the required values of gain kj, and Rj. It is convenient to select 
the instant of switching at point b at a distance A (the integration step) from point p. 


In Fig. 5 coordinates u(t) and ({t) are also plotted up to point p, Let us now draw a vertical line through 
point b up to its intersection with curyes g(t) and ((t). Coefficients * j, and T° ji are determined on the basis 
of [5, 8) as Re T* jy = Hence be, = ab (Fig. 5). Thus: 


be 1 
~ ab (12) 


From this it is possible to determine the value of one of the coefficients k* j: OF T*j, if the absolute value 
of coefficients kj and T*j, are chosen with respect to the insensitivity of the relay element, For an ideal charac- 
teristic of a relay element. For an ideal characteristic of a relay element the values of coefficients k* ji and Tj 
can be made sufficiently small. 


For determining coefficient kis let us find the required value of coefficient k, at which the controiling mem- 
ber will be displaced from point © to a position corresponding to the new balance condition of the system. 


This value of coefficient k, is found from the integrating link equation in which it is convenient to use the 
method of adding the mean ordinates [6}: 


where p4(t) is the value of the control element coordinate at point d, g(t) is the value of the control element 
coordinate at point $, p(t) is the mean value of the controlled coordinate over interval bp, which is equal to 
the value of the integration step. 


Hence: 


(14) 


Since the gain kis is known and equal to k; + 14, we can obtain by substituting the value k, from formula 
(14) the expression: 


+ 1 
(15 
hs, ) 
or 
(16) 


For determining coefficients k%, and T5: let us draw through point p a vertical straight line up to its cros- 
sing with curves ¢(t) and Let us determine from relationship = K*j9/ = 9/¢ coefficient 


ja = 17 
= = SP ~~ 0.66. (17) 


Hence, taking into consideration the relay element dead zone,we obtain coefficients k%z and T%,. It is 


Ti, 
e 
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easy to determine coefficient kjo, knowing the value of coefficients ig and 1/6; 


Thus, all the coefficients of the correcting device for the case when r < t* have been determined. 


On the basis of the above, it is possible to formulate a rule for adjusting the correcting device for the auto- 
matic control system under consideration, 


he 


i\ 


WA 
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| Fig. 6 


1, We choose the regular gain kp in such a manner that the first half-wave of the transition process due 
! to one of the types of disturbances satisfies the control quality requirements and that r does not exceed the value 
of t*, 


2. The correcting device coefficients are calculated in the following order: 


a) Assuming the value of 1/5 equal to Rp, we draw a line equal to quantity r, from the end of the first 
half-wave in the direction of decreasing time. 


b) In the interval between ¢ = 0 and the end of the line (=T ), we select a point from which we draw a 
vertical line to the intersection of curves ¢(t), ¢(t) and (t), and obtain the value of coefficients k* jy, and Th» 
a from formula (12), 


c) From the end of that line (= r) which is nearest to the origin of the coordinates, we draw a vertical 
line which extends to the intersection of at) and ,(t), and determine coefficients Ri. and T* je from formula 
(17). 


| 1 164, (18) 
| 


e) From formulas (15) and (16) we determine coefficient kya. 
f) From formula (18) we determine coefficient kis. 


So far we have examined automatic control systems working in a stable state, i.e., transient processes due 
to disturbance f(t) were investigated, It should be noted that it is possible to obtain similar results by approach- 
ing the problem from the point of view of its structure and the technique of adjusting the correcting device, by 
examining the effect produced on the. system by disturbance f*(t) (i.e., with delay in the controlled member), 


Let us now examine the dynamic properties of a control system with nonlinear correction working in a track- 
ing condition, 


For this purpose let us return to the previously examined system with linear correction (the case of r < t*) 
and, leaving the adjustment of the correcting device the same as before, let us change the point of application 
and the value of the disturbance, In other words, let us change in one jump the value of the controlling effect 
g(t) and with the help of the previously calculated coefficients let us construct the transient processes in the system, 
It will be seen from Fig. 6 that the nature of the transient process remains the same, i.e., that it will proceed 
without the necessity of readjustment, Thus, for the regulator gain k, at which < t*, the correcting device does 
not require readjustment for different scales and points of application, In other words, the correction efficiency 
is the same both for working in a stable state and tracking conditions, 
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OPTIMUM TRANSIENT PROCESSES IN SYSTEMS WITH 
A RESTRICTED THIRD DERIVATIVE 


A. A, Pavlov 
(Moscow) 


Optimum transient processes in a system with a restricted third derivative of 
the controlled variable are examined. 


The form of optimum transient processes due to step disturbances is deter- 
mined and a possible method of synthesizing the optimum controlling part of the 
system is given. 


INTRODUCTION 


This article deals with an optimum relay automatic control system of the third order, with a linear part, 
which consists of three series connected integrating links. Figure 1 shows the block schematic of this system. The 
controlled object consists of two integrating links.* 


A constant speed servomotor controlled by a relay is used as the actuating mechanism in this system* *. 
The third derivative of the controlled variable in this Automatic Control System (ACS) is restricted. This restric- 
tion is due precisely to the presence in the system of a contant speed servomotor. 


The determination of the forms of optimum transient processes under any initial conditions is the first and 
most important problem of this paper, since the solution of this problem will provide the limiting speed of opera- 
tion of the system. It is however, insufficient to know the form of the optimum control process; it is also neces- 
sary to know how to produce such a process. Thus, arises the second most important problem, namely, the synthesis 
of an optimum controlling part of the ACS under discussion, The solution of the second problem is given in this 
article only for a narrow range of initial conditions which correspond to the effect on the system of step distur- 


bances g(t), f(t), and f(t). 


Let us note that the problem of synthesis is solved directly on the basis of known formulas of an optimum 
transient process, This circumstance determines the order of presenting the contents of the article. 


The behavior of an optimum ACS in a near balance condition is not examined, since this constitutes a spe~ 
cial problem and requires special investigation. In addition, it is possible, if required, to provide the system with 


_the means of switching over when approaching balance from the optimum control relationship to any other rela- 


tionship, which would provide local stability or the required frequency and amplitude of self-oscillations. 


At present there exist a large number of works on the theory of optimum systems of the third and higher 
orders, In these works, especially in [2}-[19] a number of very important conclusions on the synthesis of such sys- 
tems has been produced, A, A. Fel'dbaum [3] was the first to synthesize an optimum ACS with a restricted third 
derivative of the controlled variable. He obtained an expression for a transition surface, which ensured an optimum 


* An airplane with a set course can serve as an example of such a controlled object [1] where ft) = 0. 
** For instance, a hydraulic servomotor. 
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Fig. 1. Block schematic of the system. I) therelay; II) the 
servomotor; III) the object of control; C) the computer. 


movement for a representative point in phase space to the origin of the coordinates with arbitrary initial condi- 
tions. In this article a different approach to the solution of the synthesis problem is taken, A simpler expression 
for the law of optimum control is obtained which provides for an optimum operation with step disturbances only. 
With disturbances of other types such as linear, quadratic or cubic, the system will not provide optimum operation, 
although the quality of transient processes will remain high. Considering the simplification of the controlling 
part of the ACS it should be considered that this solution of the synthesis problem is completely justifiable in cer- 
tain cases. 


This work was carried out in the Institute of Automation and Remote Control of the Acad, Sci, USSR under 
the guidance of corresponding member of the Acad, Sci. USSR, B. N, Petrov. 


Systems Equation, Shape of the Optimum Transient Process 


The Automatic Control System under consideration (Fig. 1) is represented by the system of equations: 


Top = halt), r= g(t)—9, 


+1 at «>0, 
=—pt+h(t), = Oar o=0, (1) 
—iat «<0, 
= 8K,, o = @(z, x, 2). 


The following notations have been adopted here: ¢ is the controlled coordinate (the output coordinate of the 
object), n is the intermediate coordinate of the object, mu is the coordinate of the controlling member, Ky is the 
absolute value of the output relay coordinate, x is the error, o is the output coordinate of the computer, © (x ,x,z) 
is the control law, which represent the unknown function of the error x and its derivatives X and X, and is to be 
determined; Ty, Ths and T,, are the time constants of the object and the servomotor, g(t) is the controlling dis- 
turbance, and f ,(t) and f(t) are disturbances, 


Excluding the intermediate coordinates from system (1), we obtain for error x the following equation: 
where 
It is easily seen from (2) that in a stable state error x is always equal to zero for condition: 
K, > |G(?)|. (3) 


This condition must hold from a certain instant, to the end of the transient process and up to the beginning 
of the next transient process, due to a new disturbance, 
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Thus, inequality (3) provides for the given ACS the type of “permissible "disturbances.* 


Condition (3) imposes restrictions on the values of the external disturbances’ derivatives A(t), iT) and 
§(t). Let us restrict ourselves to disturbances known to be permissible for which G(t) = 0. In this case we can write 
instead of (2): 


z= —8M, (4) 


where 


Introducing notations x =y and X = z and integrating Equation (4), we obtain an expression for the projections 
of the phase trajectories onto the coordinate planes of the phase space y, 2; x, y and x, z, 


_ It is known from the works [2, 6 and 7] that the optimum transient process of a third order relay ACS , the 
roots of whose characteristic equation linear part are real, consists of no more than three intervals and that the - 
signs of the relay functions of neighboring intervals alternate. Accordingly, Table 1 gives equations of the phase 
trajectory projections in intervals, 


Let us note that in Table 1 and henceforth we have adopted the notation 4 9 = sign[{®(X9, yo, Z9)]. Thus by 
substituting 5» = +1 or 59 = —1 (depending on the initial conditions) with corresponding Cy», Cg, and Cy, Cy 
in formulas of Table 1 or other formulas, we shall never obtain negative expressions under the sign of the radical, 
Constants Cy and Cg. are determined by the initial conditions and constants Cy, and Cy by the coordinates of the 
switching point N,(x;, y;, 2), but they can also be expressed by means of the initial conditions, since there is a_ 
direct connection between Cy; and Cyp, Cgg as = can be seen from Table 1. 


TABLE 1 
my 
pe Phase intrajectory projection Remarks 
‘ocess 
¥=— Bear + Cro. é 2? 
= Vo 
Cu = + 


z= + Cw. 


2C Cw— Cn 


+ Cn ar + On, 


— 


* Using A, A, Fel'dbaum's terminology [2]. 


K, 
| 
| 2 
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TABLE 2 


t 0 ty | ty te t+ te + ty 
Cc Cc 
Yo yy 0 
5Cy+C 

| — 2, | | 


The equations of the phase trajectory projection express the properties of the phase trajectory, “flows™ to 
"compress", with the movementof the representative point towards the origin of the coordinates. Hence irrespec- 
tive of the initial conditions thé representative point moves in the third interval along one of the two phase trajec- 
tories (according to the two possible values of 5) = +1 and 5,= ~1) which run to the origin of the coordinates. 
Such trajectories, which pass through the origin of the coordinates we shall call henceforth “zero” trajectories. 


The shape of the optimum transient process with given initial conditions can be determined if it is possible 
to find the coordinates in the phase space of switching points N,(x;, yy, 21) and Ng (Xg, yg, Zg). Table 1 provides 
the required data, In fact by solving simultaneously the phase trajectory projection equations for the first and 
second intervals and then for the second and third intervals, we shall obtain expressions of the required switching 
point coordinates, Neither is it difficult to find the relationships for calculating the duration of intervals of the 
optimum transient process t,, t, and t, and thus for the entire control time top = i + tet ty. 


All the data required for plotting the optimum transient process with any initial conditions are summarized 
in Table, 2, 


The quality of a system is normally judged by its transient process due to a disturbance of the nature of a 
single jump. By studying separately the reaction of the system to disturbances: 


g(t) = Ag, (t) = By, f2(t) = D, (5) 


and using the data of Table 2 it is possible to obtain the required relations for determining the shape and time of 
the optimum transient process due to the effect of such step disturbances, 


These relationships have been summz-rized in Tables 3-5. 


It should be noted that all the expressions in Table 4 are approximate, since they have been obtained by 
an approximate solution of equation 


— Cu = — (6) 


An analysis of the equation shows that such an approximate solution when the system is affectedby distur- 
bance f = By is 


Cy = Ci. (7) 


In fact, if in this instance we substitute solution (7) into Equation (6) expressing Cyg and Cy, by means of 2», 
z3 
we shall obtain for the left hand side expression -<2 s 0,993 and for the right hand expression 
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Fig. 2. Projections of optimum phase trajectories when the system re- 
sponds to g(t) = Ag = 1.0 with M = 1.0 1/sec*, 
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Fig. 3. Optimum transient process when responding to the con- 
trolling effect g(t) = Ag = 1.0 with M = 1.0 1/sec®, 


3 3 
- s7iya © ~0-353 ch. Thus the difference between the left hand and right hand side amounts only to about 


3 
0.02 rv . It means that solution (7) is practically equal to the exact solution, Calculations of Cy, for a wide 
range of By changes confirms the validity of above consideration, 


On the basis of the data in Tables 3-5, projections of optimum phase trajectories and corresponding optimum 
transient processes,when the system responds to step variables of the type (5),have been plotted in Figs. 2-7. The 
character of phase trajectories and the shape of the optimum transient process does not change with the variations 
os the scale of step disturbances, 
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Fig. 4. Projection of optimum phase trajectories when the sys- 
tem responds to disturbance f,(t) = By = 1.0 with M = 1.0 1/sec®, 


Fig. 5. Optimum transient process of response to disturbances f ;(t) = By = 1.0 with 
M = 1,0 1/sec*, 
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Fig. 6. Projections of optimum phase trajectories when the system responds to dis- 
turbances f)(t) = Dy = 1.0 with M= 1.0 1/ sec’, 
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Fig. 7, Optimum transient process when responding to disturbance 
f = Dy with M = 1.0 1/ sec’, 


Synthesis of an Optimum Controlling Part of an ACS under the Effect of Step Dis- 
turbances 


Below a method of synthesizing an optimum control partof an ACS is described, a method which differs from 
those described in the literature [3-19], Above method has already been mentioned in the introduction and will 
now be described in greater detail. 


Let us examine the piiase picture of a second order ACS (Fig 8), It is known that the transients in a second 
order relay ACS will be optimum, i.e., will be complete in two intervals if on the phase plane of the error x and 
its derivative y = X the transition line coincides with the corresponding branches of the zero phase trajectories for 
the values 6 = +1 and 6 = —1, If the equation of such a transition line has the form y + F(x) = 0, the optimum 
control law will, obviously, be o = y + F(x), 
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Fig. 8. Phase picture of the optimum second 
order relay system, whose linear part consists 
of a series connection of two integrating links. 


It follows from Fig. 8 that it is also possible to attain 
in a second order system optimum transient processes for any 
initial condition by other means. In fact by assuming the 
line y + a,x = 0 as a transition line for one pair of initial 
conditions, N; Yo) and N°, ye) and the line y + 
+ a,x = 0 as a transition line for any other pair of initial 
conditions Ny (Xo, yo) and N*, (X9, yo), we shall obtain in the 
system the same optimum transient processes, which are 
completed in two intervals, as in the case when the transition 
line equation coincided with the zero phase trajectory 
equation. 


It is easy to show that there is a close connection 
between the coefficient a in the transition line equation 
y + ax = 0 and the initial conditions xg and yp. Hence it 
it possible to find the optimum control law o = y + ax, 
which would ensure an optimum transient process in the 
system for any initial conditions with variations of a 
corresponding to the variations of the initial conditions, 
Since, however, for second order relay ACS is usually pos- 
sible to obtain a fairly simple expression for the zero phase 


trajectories equations, above method has not been adopted for these cases, For optimum systems of the third and 
higher order, however, when the nonlinear equation of the transition hypersurface becomes very complicated, 
above method of determining the optimum control law may become useful. 


For third order ACS with given initial conditions x9, yp and 2» it is always possible to find in the phase space 


a transition plane,* 


z+ay+bz=0, 


a>0, b>0, (8) 


which would provide an optimum movement for the representative point to the origin of the coordinates (in two 
transitions) when the coefficients a and b must be chosen to equal: 


Tike 
Yi22— 


— 
Vite — (9) 


where and Xg, yg, Zg are the transition point coordinates. 


Thus, the optimum control law for any third order relay ACS, the roots of whose characteristic equation are 


real, is written in the form: 


o=r+ 


10 
9i22— y — (10) 


For above ACS with a restricted third derivative the data of Tables 3-5 provide very simple expressions for 
coefficients a and b and hence for the optimum control law (10), Each of these tables contains corresponding 
expressions for the optimum control laws whose application provides an optimum response to step effects (5). 


A synthesis of an optimum control part of anACS is theoretically quite possible by applying the control 
laws given in Tables 3-5, but the values of Cy) or Cy) must be known, i.e., it is necessary to remember the ini- 
tial conditions, which is not always possible. However, the relations between Cay,,Cyq and Cy, are also given in 


Tables 3-5 and the relationship 


42M | yl 


- © dt is shown in the supplement that the transition plane represented by (8), when the system responds to step vari- 
ables, is in fact crossed by the optimum trajectory at three points only (including the origin of the coordinates) 
so that in applying the control law (10) there will be no false switching. : 
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ari- 
ates) 


which holds in the second interval of the optimum transient process (see Table 1, Figs, 2 and 4) provides by means 
of direct substitution for all the three cases of application of step effects to the system, an optimum control law 
in the form: 


2? + 2M + + 2M (12) 
Thus, expression (12) should be considered as final and the expressions for o given in Tables 3-6 as inter- 
mediate on the basis of which the optimum control law in the form of (12) has been obtained, 


Let us stress once more that the controlling part of an ACS designed according to the control law (12) can 
provide optimum response for step effects only. Let us note also that the numerical values of coefficients a; and 
8; depend on the point of application of the disturbance, Thus, in the response to the controlling effect g(t) = Ag 


we have: a = = » By= or. In the response to disturbances in loading f,(t) = By we shall obtain, respectively, 


Simulation of the System 


Figure 9 shows a block schematic* of an optimum ACS with a control law (12), The optimum ACS was 
simulated according to this schematic by a electronic simulator of the type EMU-3 of the Institute of Automa- 
tion and Remote Control of the Academy of Sciences, USSR, with the help of a set of nonlinear units KNB SAM.* * 
Polarized relay RP-5 U 1722001 was used as a relay element, 


“hg 


Fig. 9. Block schematic of an optimum system with restricted third derivative and o = 
= x +ajvZ,+ OM] yl y + + 2M| y|)z. adders; NC-1, NC-2 and NC-3are 
nonlinear converters; M-1 and M-2 are multipliers, 


The simulation has shown that it is necessary to adjust the conwolling part of the ACS by changing the 
numerical values of coefficients a; and 84 from those obtained by calculation, This is due in the first place to 
expression (12) being equal to zero not only at the points of transition Ny(x;, yy, 21) and Ng(Xg, Yg, 2g) but also 


* If disturbances f,(t) and ft) can be measured directly, it is possible to dispense with differentiating links, when 
determining y and z, since 


**SAM = Computing and Analytical Machine. 
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Fig. 10. Oscillograms of the transient processes, a) is the response to the controlling 
effect g(t) = Ay = 1.0 (calculated time t,,, = 2.516 sec), b) is the response to disturb- 
ance f;(t) = By = 1.0 (calculated time t,,,, = 2.21 sec); c) is the response to disturban- 
ces (calculated time topt = 2.92 sec), The system parameters: T g = 0.5 sec, T, = 

= Ty = 1.0 sec, K, = 1.0, and M= 2 1/sec’, 
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remaining identically equal to zero* along the entire third interval of the optimum transient process (see Figs. 

3 and 5). Hence in the third interval the controlling effect o does not reach the relay (in practice o #0). This 
leads to a blurred operation of the relay at the transition point NgX», yz, 22) which produces considerable distor- 
tion of the shape of the optimum transient process in the third interval. Imperfections in the relay and errors in 
the nonlinear units tend to produce self-oscillations in the near-balance position. Changing of the coefficients 

a; and Bj and the control law (12) by about 3-5%,as compared with their calculated values,provides a clean 
operation of the relay at point NXg, y2, Zg), without causing any appreciable distortion of the shape of optimum 
transient processes for a wide range of initial conditions, Figure 10 shows the oscillograms of the transient proces- 
ses due to step effects, It will be seen from these oscillograms that the transient processes obtained are very close 
to the optimum values, although in the near-balance position the system is in an oscillatory condition, 


Oscillogram a shows the response of the system to a step controlling effect of the type g(t) = Ag with A, = 
= 1,0.It will be seen from this oscillogram that there is hardly any difference between the duration of the transi- 
ent process recorded on the oscillogram and the optimum duration obtained by calculation t,.. In fact the dif- 
ference amounts only to 0.004 sec or 0.16% of t,,,. Oscillogram b shows the response to disturbance f;(t) = By 
with By = 1,0, In this case the difference between the actual duration of the transient and its calculated optimum 
value amounts to 0.21 sec or about 9.5% of topt. In oscillogram c, which shows the transient with a system res- 
ponse to a step effect f(t) = Dy with Dy = 1. O the difference between the actual and optimum calculated regu- 
lation time amounts to 0.39 sec or about 12% of t,,,. This difference between the actual duration of regulation 
and its optimum calculated value is due to errors in the nonlinear units of the computer, which shapes the opti- 
mum control law o, With different points of application of disturbances similar in scale,different nonlinear unit 
characteristic sections are used which have different errors and different effects on the final value of the control- 
ling signal o. As the result of above circumstances, the imperfect operation of the relay element and other un- 
accounted simulation errors produce the difference (in percent) between the actual and the calculated optimum 
regulation time, 


Simulation has shown that the frequency of self-oscillations increases and their amplitude decreases with 
decreasing constant T,,, i.e., with raising the maximum value of restricted third derivative IX max |= M. Simu- 
lation has also shown that a system with a control law (12) will respond well not only to step effects, but also to 
effects of any other shape, for instance, linear or parabolic, Then, however, the transient processes will no longer 
be optimum, although their quality will remain high, 


SUMMARY 


1, In this article the shape of an optimum transient process in a system with a restricted third derivative 
of the controlled variable is obtained and relations which provide the duration of regulation with any initial con- 
ditions are found. Above relations provide an evaluation of the difference between transient processes in a non- 
optimum Automatic Control System (ACS) and those in a similar optimum one, 


2. Expressions for an optimum control law are given and a block schematic of an optimum ACS provided, 
in which optimum transient processes are attained by means of relatively simple means when the system is sub- 
jected to step disturbances. 


3. Simulation of an ACS produced good agreement between the results thus obtained and those calculated 
and confirmed all the theoretical conclusions. Moreover, simulation provided the possibility of studying the sys- 
tem's operation when the disturbing effect is linear or parabolic, Although under the latter conditions,the transi- 
ent processes in the system are no longer optimum,their quality remains high. 


Supplement 


Let us prove that for the cases examined above of applying to the system step effects of the type of (5), 


*Equality o = 0 in the third interval of the optimum transient process for the particular cases of application of 
step effects examined above, namely g(t) = Ag and f ;(t) = By is easy to prove if in expression (12) with appropriate 
a, and 6, one substitutes from Table 1 the following values for y and x: 
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the transition plane 
«+ ay-+bz=0 3) 


is crossed by the optimum phase trajectory in three points only, t,e., there is no additional switching in the sys- 
tem, 


It will be seen from Table 1 that the projections of the optimum phase trajectory onto the coordinate planes 
y, 2 and x, z can be represented in the form of y = yj(z) and x = $,(z) where i is the number of the optimum pro- 
cess interval, By substituting these expressions into the equation of the transition plane (I) we obtain for each in- 
terval a cubic equation in terms of the coordinate z: 


+ 9542" + + 9. (Il) 


Since the transition plane (I) is drawn through the origin of the coordinates 0(0, 0, 0) and the transition 
points Ny(x;, yy, Z,) and Ng (X3, yg, Zg), Certain roots of the Equation (II) are known in advance for each interval, 
Such roots are: for the first interval z = z;, for the second interval z = z, and z = Z, and for the third interval 
Z = and z = Zp. Since the transition points N,(x,, y;, and yg, Zz) serve as boundaries of the optimum 
transient process intervals, z = z, and z = 2, are roots of Equation (II) in the first and second and the second and 
third intervals respectively. 


Since part of the roots of Equation (II) is known in advance, it is easy to find the remaining roots, 


Table 6 shows the solution of Equation (II) in intervals for various cases of applying step effects to the sys- 
tem. In each interval we obtain for the cubic Equation (II) three roots, each of which corresponds to a crossing 
point of the optimum phase trajectory with the transition plane. The analysis of roots z';, z";, Z':, and z", shows 
that the representative point moving along the optimum phase trajectory could cross plane (I) in points N',, N";, 
N', and N”, only if its movement were in an opposite direction to its actual direction of movement which is de- 
termined by the sign of 5 = 41 in each interval, All this can be demonstrated in Figs. 2, 4 and 6, which were 
plotted for initial conditions corresponding to 59= +1. 


Thus, the transition plane (I) with the corresponding values of coefficients a and b does provide an optimum 
response to step disturbances, 
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AUTOMATIC CORRECTION OF COMPENSATOR PARAMETER 
IN FEEDFORWARD SYSTEMS 


S. A. Doganovskii 
(Moscow) 


The problems of applying an automatic optimizer for correcting compensator 
parameters in a feedforward system are examined. 


The optimizer circuit and examples of its use are given, 


Posing the Problem 


It is known that disturbances affecting a system can be balanced out by feedforward adjustment and by com- 
pensation for controlled variable deviations from a set value [1]. In a feedforward system the regulating element 
responds to disturbances before the controlled variable and in theory any error, including zero error, can be ob- 
tained, i.e,, conditions providing an invariable controlled quantity can be attained (see [2]). It is difficult, how- 
ever,.to attain these conditions in practice, owing to changes in the characteristics of the controlled member and 
inaccuracies of adjustment in the regulator compensating devices, These are the reasons for a relatively restric- 
ted application of feedforward systems. : 


Automatic correction of the compensating device parameters for 


oe ae a changes in the controlled member characteristics is a new departure in 
se | the development of automatic control systems which provides for a much 
D } z lien I y wider application of feedforward systems [3, 4]. The problems of de- 
iC = signing such systems, however, have hardly been touched upon in litera- 
| ture; the object of the present article is to examine these questions in 
connection with controlling the strip thickness in a cold rolling mill, 


Fig. 1 2. General Principles of Designing a System for Auto- 
matic Correction of Compensator Parameters 


Figure 1 shows the block schematic of an automatic system for correcting compensating device parameters, 
This device in general consists of analog computer C,, which simulates the controlled member CM characteris- 
tics, Disturbance D is fed to the input of C, and at its output the controlling x is obtained, and applied through 
the second input to the controlled member. The correction of the compensating computer C; parameters is car- 
ried out by means of the optimizing (correcting) computer C, whose input is fed with the controlled variable y. 


When the characteristics of the controlled member change, and thus the conditions of invariance are dis- 
turbed, the optimizing computer C, reacts on the parameters usually of both the linear and nonlinear parts of 
computer C,, thus correcting its operation. This correction of C, parameters is carried out on the basis of a cer- 
tain mean value of the controlled variable, for instance, f[y(t)t, f [y(t)'dt or f y(t)B(t)dt. 


Such a device is in fact a self-adjusting system and the basic element of C, consistsof a regulator for selec- 
ting automaticallyC, parameters which provide minimum variations in the controlled quantity at the output of 
the controlled member, Such a regulator has been named an automatic optimum regulator (optimizer [5-8]). 
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a The choice of the required parameters for C, is made by the designer 
on the basis of a known approximate description of the process occur- 
ring in the controlled member and accurate correction of these param- 
eters is achieved by C, on the basis of gradient scanning, In this 
instance the gradient scanning method provides the smallest stable- 
state mean deviation from the minimum and the best tracking of 


the slowly changing controlled member characteristics for the pur- 
LU pose of obtaining minimum variations in the controlled variable [6]. 
{ 3. Development of the Optimizing Computer Cy, 
CG Circuit 


The circuit shown in Fig, 2 is a possible version of optimizer 
Fig. 2 A, which scans by means of the gradient method, The circuit con- 

sists of the following units: a memory device MD, actuating (inte- 
grating)device AD,which operates the controlled member CM, cycle 
generator CG (of the control unit) and a device for obtaining a 
difference between two quantities (of the unbalance unit). The con- 
trolled member static characteristic, which has but one extremal 
point, is shown in Fig. 3. 


Let the condition of the system correspond to X», Ye on the con- 
trolled member characteristic (Fig. 3), The circuit starts operating 
when CG closes contacts K, and Ky (leaving Ks and K, open), thus 
feeding to the input of CM a test step Ax and making MD record 
the newly established value y = yy + Ay (point 2), The second (opera- 
tive step causes the opening of contacts K, and K,and closing of Ky 
and Ky, thus producing at the output of the device an unbalanced Ay, 


Fig. 3, 


Ay = (y, + Ay) — yo. (1) 


which is fed to the input of AD, producing a change in the output value of x. The rate of change of x being pro- 
portional to Ay. 


Let AD in the above case provide a decrease in x; Then the system will move towards the extremal until 
(point 3) after the completion of a cycle CG again closes contacts K, and Kg and opens Ks and K,, feeds a test 
step Ax to the input of CM and makes MD record the new established value of y (point 4), In the next cycle of 
operation of the CG circuit a new value of Ay is provided which is smaller than the preceding one,and hence the 
rate of change of x is also decreased, The movement of the system in the direction of the extremal will proceed 
until point A is passed and a negative Ay increment appears, which reverses the direction of scanning, The rate 
of change of x drops as the system approaches the extremal since Ay is decreasing. This ensures a smooth ap- 
proach and stabilization of the system at the extremal, Any change in the characteristic of the object will dis- 
place point A, which will produce a renewed scanning and stabilization, A change in the operating condition 
(scanning of the maximum or minimum) is taken care of by the changing of the sign of the fractional step Ax. 
Scanning of an extremal with n input variables of CM can be achieved by a consecutive connection of the un- 
balance unitoutput (point a in Fig. 2) to the input of AD following every two CG cycles of operation, This is a 
simplified gradient scanning method of operation and has important advantages with respect to the simplicity of 
the equipment required. 


Figure 4 shows the trajectory of the system's movement to the extremal by the usual (dotted line) and sim- 
plified (full line) gradient methods, It will be seen that in the latter case after the determination of the partial 
derivate 8Q/0x, at point My the system moves to point M, with the coordinates Q(Xy9 + 4X9, Xgp) and the values 
of the step is proportional to the derivative, i.e., 


(0Q /O2;)y,» 
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where m is the constant coefficient, At point M, the partial derivative 
8Q/0x; is evaluated, and is made a step to point Mg with the coordi- 
nates Q(Xi9 + AX4o9, Xg9 + AXg9), and the value of the step is proportional 
to 8Q/ axg. At point Mg the derivative 8Q/ dx, is again determined 
and a step made, whose value is proportional to 0Q /@x;, thus bring- 
ing the system to point mg, etc, The movement of the system will 
stop when each of the partial vatives becomes smaller than a pre- 
determined small quantity, i.e., when the system is stabilized inthe 
neighborhood of the extremal of function Q, One of disadvantages of 
Fig. 4 this method consists in the prolongation of scanning (approximately by 
the ratio of two as compared with the ordinary method). However, 
for a small number of variables (n = 2) the simplicity of the equipment required is often the deciding factor. 


When the scanning of extremals is made in an object whose control can be attained without the use of elec- 
tromechanical devices, the actuating element can consist of an electronic integrator and the circuit will assume 
the form shown in Fig. 5. Here capacity Cy is used as a MD and an electronic integrator as an AD, which is based 
on an operational dc amplifier (A,). The well-known relay circuit whose time constant RyC, determines the 
duration of the cycles T,, is used as the CG, Contact Ky combines the functions of contacts Kg, Ks, and K, (Fig. 2), 


Let us now find the basic relations which determine the time and accuracy of scanning for the extremal of 
CM, The stable error at the input of the member will obviously amount to $Ax. = The scanning duration will 
depend on the form of the CM characteristic, the initial deviation of the input variable of the system from the 
extremal X9, the value of the test step x, the gain k,, the time constant RC of the integrating (actuating) device 
and the duration of the test Tyg and operating step T,, of the cycle generator,* The scanning duration can be 
found from the formula: 


To =n(T,, os )s (3) 
where n is the number of scanning steps during which the system will enter zone Ax /2, 


Fig. 5 Fig. 6 


Let us determine n for the case when the CM characteristic has the form: 
y = ko? + ky. 
In this case the increment of the output voltage of AD is determined by the expression: 
Tts 


8, \ Aydt = = On, (5) 


where Ay; is the increment of the function at the i-th point. 


* If the object has inertia it is necessary to select the duration of cycles in such a manner that the transients in 
CM have enough time to be completed. 


1010 


| 
, 


Starting with the initial point x, (Fig. 6) the values of 5x; will be: 


bx, = 
bz, = = 9 (x, — 824) = Oz, (1 — 4), 
= Ox, = 6 (x, — = Oz, (1 — 6) = (1 — 8)’, (6) 


Oty, = 0 — = (4 — 0) = Ox, (1 — 6)", 


Since the n-th step takes the system into the zone Ax / 2, the equality 


==, (1) 
holds and serves to find the value of n *: 
lg 
2x0 
ig — 6) (8) 


When the controlled member characteristic has the form: 


(9) 
number n is determined from the formula: 
— (142) 
n=(1 =) (10) 


where Q = <8 Tts k Ax. 


Fig. 7 


For instance with xX) = 38 v, x = 2 v, k = 0,01, RC = 6 sec, T,, = 1 sec, T, = 4 sec, k, = 23.2, we obtain the 
number of steps n = 21 and Ty, = 85 sec, The oscillogram shown in Fig, 7 represents precisely this condition of 


operation, It will also be seen that the number of test steps made by the system during scanning is close to the 
calculated number. 


* See also [6]. 
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The application of this optimizer system for correcting the C, parametes requires a reliable filtering of 
the controlled output variable which forms the correcting signal g(t), The use of a line filter Leads in practice 
to large increases in the duration of scanning Ty since the inertia of the system is increased. 


A pulse filter in the form of an integrating device can be used for filtering out the oscillations of the con- 
trolled variable, A preliminary filtering of the oscillations, whose frequencies exceed the operating range of the 
controlled variable oscillations is achieved in the operational dc amplifier operating as an inertia link. 


The scanning proceeds as described previously, except that before each step the capacity of the pulse filter 
is discharged. 


Figure 8 shows an oscillogram illustrating the procedure of correcting the compensating C, computer gain 
by means of the optimizing computer Cy described above, Here the output controlled variable y(t) due to a sinus- 
oidal disturbance B(t) is shown in the left hand side of the diagram, and the gain of the compensating C, computer 
is close to zero, Next C corrects the gain of C,. First dy /éx, is determined which is reflected in the change of 
the y(t) amplitude, due to the setting of the test step and the change in the gain of C,, This is followed by the 
operational step Axy9 = —m(dy/0x,), which decreases the amplitude of y(t). The determination of y/ x, again 
changes the amplitude of y(t),and a new operational step leads to its further decrease, etc. The right hand side 
of the oscillogram shows the values of the controlled variable y(t) obtained at the end of the correction of C, gain, 
In this system of correction of C,; parameters is carried out by means of the variance of the controlled variable 
oscillations, i.e., the correcting signal is shaped according to the expression: 


+ 


sec 
Fig. 8 


It should be said in passing that sometimes (when a number of uncontrolled disturbances are affecting the 
system) it is better to correct according to correlated function of the output and input variables of the controlled 
member, which provides a correction of C, only by the controlled disturbance, thus avoiding “hunting” due to 
uncontrolled disturbances (see [4]). In actual systems there often exists one main disturbance and a compensa- 
tion of oscillations over a wide frequency range is required. 


4. An Example of the Use of an Optimizer for Correcting Compensating C, Com- 
puter Parameters 


Compensation of rolled metal thickness variations in a continuous cold rolling mill by means of a com- 
pensating C, computer [11] is a good example of applying an optimizer for solving the problem of correcting 
parameters by means of a computer, .Figure 9 shows, schematically, four stands of the mill rolling a strip of metal. 
The strip thickness H, at the output of the mill can vary both due to changes in the operation of the mill (stresses 
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Fig. 9 


in the rolls, their heating, changes in lubrication, etc,) and to the properties of the rolled metal, especially due 
to changes in the thickness H’, of the strip at the input. Changes in the so-called feeding thickness H", are usually 
due to errors in the operation of the hot rolling mill, which supplies the rolled strips. 


Large amplitude but low frequency feeder-strip thickness variations can be eliminated in the first stand by 
means of a thickness regulator, which operates the pressure device PD, on this stand, Thickness gage M, follows 
the first stand at a distance J, , If, for instance, the thickness of the rolled metal is larger than the required value, 
the drive of PD, is actuated and decreases the distance between the rolls, thus reducing the thickness to the re- 
quired value, This control system works slowly, due to the dead zone, time delay, inertia of the pressure device 
and low speed of roller displacement, Feeding-strip variations of a higher frequency and lower amplitude which 
remain after the pressure device regulation, can be eliminated by controlling, for instance, the tension F, between 
the second and third stands, The tension variation can be attained by changing the angular velocity of the third 
stand rolls (more precisely by changing the relative speed of the second and third stand rolls), Since the tension 
of the strip between the stands, the torques, etc,, must not exceed certain limits, it is possible to compensate by 
this method only some of the thickness variations which lie within certain limits, In order to vary correctly veloc- 
ity Qs it is necessary to measure the strip thickness by means of micrometer M, and suppy the measurement 

results to the compensating computer C,, which calcu- 
p Mn lates the correction increment in y for the set value of 
Qg9. Computer C, must reckon with both the static and 
dynamic processes in the rolling mill, which are very 
6P3 complicated. An approximate representation of this proc - 
y ess by means of a system of linear equations, which 
were derived comparatively recently, (see [9, 10]), pro- 
MV vide data for designing a rolling mill model, using it 
for simulating the required mill characteristics and de- 
veloping the C, computer [11]. Since the processes in 
the mill are represented by a system of linear equations 
the following expression will hold; 


! 


hour = K hy (jo) hy + Ky (o) y, (12) 


where hy, are strip thickness variations measured by micrometer My, hoy, are the strip thickness variations at the 
output of the mill, Knyy(jw) and Ky(jw) are the amplitude and phase characteristics (a.p.c.) of the mill simulator 
with respect to inputs hy, and y. On the basis of complete compensation of strip-thickness variations at the mill 
output we find: 
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Fig. 11 


Ky. (jo) 


where Ko(jw) is the a.p.c, of the compensating C, computer. The synthesis 
of the structure of such a computer can be simplified following an analysis 
of the a.p.c, obtained. It can be shown that for 


ky 
Kn, (jo) = (14) 
and 
ke 
Ky (jo) = ar, (15) 
we have: 


K, (jo) K (jw) = 


and that computer C, consists of a series connected delay unit, simulating 
the factor e)“(7~ 79 and link K(jw). 


Thus, C, provides the possibility of making the third stand drive 
“prepare” in advance for the feeding-strip thickness variations, measured 
by M, and to compensate these variations in the approach to the third 
stand, by changing speed Q3). Such a computer provided compensation in 
a three-stand simulator of a rolling mill in the frequency range of 0 to 
2 cps with a limiting amplitude of 25 w [11]. This system, however, is 
an open loop one and changes in the mill characteristics can, therefore, 
upset compensation, Moreover, an electronic model of the mill only ap- 
proximately simulates the rolling process, and in practice C,; cannot com- 
pletely eliminate the variations of the outputstrip thickness, In such a 
case the optimizing computer C, is used, based on the circuit described 
in section 3, which provides the required correction of some of the C; 
parameters, Such parameters may include the gain k and the delay time 
t of the computer. The circuit of the optimizing computer C, is shown 
in Fig. 10. Computer C, is conventionally divided into two units; the 
automatic optimizer A proper and the minimized variable shaping unit 
MVSU. The correction of C, parameters is achieved by means of non- 
contact actuating devices and electronic integrators, consisting of opera- 
tional dc amplifiers A, and Ay, The outpui voliages of the integrators con- 
torl changes in k and r of the computer C, circuit. The test step (Ak and 
Ar) is set by means of voltage source €, 


The correction of C, parameters is carried out on the basis of the 
variance of the output strip thickness, For the shaping of this signal a 
multiplying unit MU», working as a squaring device, is used, * as well as 
an averaging unit F, which consists of a pulse filter previously described. 
The cycle generator which produces the required sequence of steps is not 
shown in the circuit, 


* The dotted line in Fig, 10 shows the connection of the MU, for correla- 
ted function correction. 
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Fig. 11 shows an oscillogram illustrating the process of correcting C, parameters by means of the optimizing 
computer Cg, described above, when the characteristics of the rolling mill simulator are changed. Here devi- 

(13) ations of the output strip thickness with sinusoidal thickness variations are shown in the left hand part of the di- 
agram where the adjustment of the compensating C, computer is not completely accurate, Next, due to the 
changed characteristics of the mill the amplitude of the strip thickness deviation oscillations increases, Follow- 

thesis ing this, C, begins to correct parameters k and r of computer Cy, At first the value of dhpy,/0k is determined, 
nalysis which is reflected in a changed amplitude of oscillations h,,, (point a in Fig, 11). Then the operational step is 


made Axio = —m, “Rout which decreases the amplitude of oscillations hoyt (point b). The determination of 


(14) dhoye/OT again changes the amplitude of oscillations ho,,,, but the operational step changes it little, due to 
closeness of rT to its set value, Next dhoy,/dk is again determined, etc, The value of the deviation of the out- 
put strip thickness is shown in the right hand side of the oscillogram and obtained from consecutive correction of 
the C, parameters k and r. It will be seen from the osicillogram that in correcting the parameters of C,, com- 
puter C, minimizes the strip thickness deviation at the output of the mill simulator with changing characteristics 

(15) of the latter, making the C, adjustment more accurate than that obtained previously manually in the absence of 


SUMMARY 


1, Automatic correction of the compensating device parameters, with changing characteristics of the con- 
(16) trolled member, is a new and promising tendency in the development of automatic control systems, 


2. The developed optimizing computer is relatively simple and capable of correcting several compensat- 
ulating ing device parameters. 


3. Experimental investigation of automatic correction of compensating device parameters by means of an 


ive electronic simulator of a rolling mill confirmed the correctness of the proposed calculating technique and synthesis 
asured of computer circuits. 
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ON DESIGNING CONTROL CIRCUITS FOR OBJECTS WITH 
PURE LAGS* 


E. L. Itskovich 
(Moscow) 


The paper considers the question involved in the design of control circuits for 
objects with lags which are acted upon by random disturbances, Circuit design to 
meet given quality requirements is based upon a study of the properties of the cor- 
relation function of the quantity to be controlled, A description is given of an ap- 
plication of the method presented to the designing of a control circuit for one of the 
output quantities of a rotating-kiln cement-roasting process. 


In the determination of the transfer functions for actual technological aggregates, there sometimes appear, 
in the block schematics, links with pure (transport) lags. Lags may be inherent either in the object itself, or in 
the transducers measuring the quantities to be controlled, Independently of the exact locus of the element de- 
laying the transmission of a disturbing stimulus, it is possible in a single-oop circuit, to consider the lag as a fac- 
tor-inherent in the given object, In this case, design of the control system must be based on the presence of a 
link with delays in the invariant (given) portion of the system, 


Below, we shall consider objects, which in actual production conditions are subject to the action of disturb- 
ing stimuli which are random functions of time. The majority of technological aggregates fall under this head- 
ing. We may cite, for example, such a ubiquitous technological aggregate as a pipe ball mill. 


The output quantity which characterizes the mill operation is the fineness of the material milling. It is 
defined by the material output by the mill and has a transport lag with respect to the controlling action, namely, 
the supply of raw material to the mill. This lag is occasioned by the translation of material in the object, and 
comprises some tens of minutes, The basic disturbances occurring in the raw material, such as variations in coarse- 
ness, hardness and other physical properties of the material, as well as variations in the rate of flow of material 
into the mill, are random functions of time and, under actual factory conditions, the temporal variations of these 
disturbances cannot be predicted. 


Other equally widely distributed objects are stoking boilers or furnaces, The output quantity which char- 
acterizes the process in a stoker, the oxygen content in the flue gases, has a transport lag of the order of tens 
of seconds, A significant portion of this lag is concentrated in the measuring device, and is determined by the 
time taken to pass the gas selected for analysis through the gas analyzer's gas-screening device. The fundamen- 
tal disturbing factors, namely, the variations in physical and chemical properties of the fuel, its make-up, the rate 
of flow of fuel and air, all these have a purely random character. 


Regulators designed for the objects cited frequently do not possess the given quality of control, despite the 
satisfactory system reaction to unit disturbances of arbitrary form, Sometimes, under the action of actual disturb- 
ing factors, the control of objects with lags is not even stable, even though the system constructed satisfies the 
stability criteria, 


* This work was reported at the seminar on probabilistic methods in automatic control, held at the IAT AN SSSR, 
on December 15, 1958. 
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Below, we shall consider the causes of the phenomena we have mentioned, and we shall determine the 


conditions which, if met, make it possible in practice to construct control systems of given accuracy for objects 
with lags acted upon by random disturbances, 


In this paper we shall consider disturbances in the class of stationary random functions of time, thus allow- 
ing ourselves to use for their analysis the concept of the correlation function, which is determined from experi- 
mental recordings of the process, However, all of the subsequent analysis also applies to thosenonstationary dis- 
turbing factors which satisfy the conditions necessary for the experimental determination of the correlation func- 
tion by smoothing one realization of the random function [1]. 


A simplified block schematic of the control system to be considered is shown in Fig. 1. At the input to the 
invariant portion of the system, made up of the links with transfer functions W,(p) ande~P7 , there act a stationary 


random disturbing stimulus, n(t), and a controlling stimulus, g(t), which in cases of controlled parameter stabili- 
zation, is constant, 


We now show that, under certain definite conditions, construction of a control system in accordance with the 
scheme in Fig. 1 is impossible. For this, we consider one of the most highly perfected methods for constructing 
control systems by means of a physical realization of the optimum transfer function, 


Let the correlation function of the disturbing stimulus have the form: 


The spectral density of this process is then defined by the expression 


Sa(w) = (2) 


For the disturbance of the form shown, and for a constant controlling stimulus, g(t) = const, the optimal 
transfer function for the system is written in the form (cf., for example, paper [2], example IV) 


(p) = 4+ K + , (3) 


where T is the time of observation and A, K and D coefficients which depend on T, if 


By knowing the transfer function of the entire system (Fig, 1) and its expression in terms of the transfer func- 
tions of the individual links 


W, (p)e—?* 
1—Wr(p) 


 (p) = 


(4) 


we can find the necessary transfer function of the controlling device: 


e?t 


By substituting (3) in (5), we find that 


Ws (P) Wil) A (i —e~Pt) +K + De~?T 
P 


(6) 


For a physical implementation of the controlling device, it is necessary to realize the transfer function in 
(6). Ways of realizing the second term of the transfer function by its logarithmic frequency characteristic are 
cited in [2]; to realize the first term, it is necessary, in particular, to implement a leading link (cf., for ex- 
ample , [3]). 
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For a given correlation function of the quantity subject to leading [in our case, the quantity x(t) }, the 
leading portion is the lower the greater is r, i.e,, the lower is the correlation of the values of the quantity sepa- 
rated by time r. Correspondingly, with an increase of r, there is also an increase in the lead error. Since the 
lead time, r, ts the lag time in the object, it is possible to consider the lead error as a function of the lag time, 
For sufficiently long lags, the values of the quantities separated by time r becomenoncorrelated for all practical 
purposes, In this case, construction of a leading link becomes impossible. Consequently, it becomes impossible 
in theory to realize physically the transfer function of the regulating device given in (6) for sufficiently large r, 


The case we have been considering means that the correlation function of the output quantity x(t) falls to 
a value close to zero during time ty (decay time of the correlation function) which is less than, or equal to, the 
lag time in the object: 


(7) 


If condition (7) holds, the very concept of controlling the process in the object by its output quantity x(t) 
becomes undefined, since the value of the quantity entering the controlling device at a given moment of time is 
in no way correlated with the simultaneous changes in the supposedly controlled process in the object. Conse- 
quently, the holding of condition (7) completely excludes the possibility of constructing a control system for the 
quantity under consideration, 


Only in the case when a comparison of the correlation function's decay time with the object lag provides 
the condition 


(8) 


does it become possible to construct both a leading link and a control system for the output quantity. The hold- 
ing of inequality (8) may serve as an existence condition for a control system for an object with lags. 


In the design of a control system for such an object, it is necessary to provide, in addition to stable system 
operation, a definite quality of control, the requirements on which are imposed by the technological conditions 
for the given object. Independently of the method of designing the regulating device, optimal tuning of the 
regulator, even if condition (8) holds, may give a mean-square error of the quantity being controlled which 
does not satisfy the requirements for control quality. Physically, this will be defined by an insufficient degree 
of correlation between values of the output quantity separated by time t= r. With this, even a regulator de- 
signed by the method of physical realization of the optimal control system cannot provide the given quality, 

_ since the lead error obtained from the leading link will be too large. 


Consequently, to obtain the necessary quality of operation of a control system, even before designing the 
controlling device it is necessary to provide a definite degree of correlation between values of the output quan- 
tity separated by time t = r, or to make sure that in this time the variations of the output quantity do not exceed 
a previously given vaiue. 


In the general case, to have the necessary accuracy of control of the quantity x(t), the holding of the fol- 
lowing condition should be provided, 


F (0), Rx > P, 


where P is a magnitude determined from the technological requirements on the accuracy of control, 


. The form of the function connecting the values of R,(0) and R,(T) in formule (9) depends both on the charac 
ter of the requirements on the qualitative indicators of the control process and on the method of creating the con- 
trolling device, 


Let us consider the choice of the form of the function in formula (9) for several definite problems. 


1, The case of the physical realization of the optimal transfer function with a leading link, The problem 
of providing a given quality of control may be posed as the problem of obtaining a definite (not less than some 
given) coefficient of mutual relationship (P,) between values of the output quantity separated by time Tr. In other 
words, the following conditions must hold: 
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R(t). 
> P, (0 <P; <A). (10) 


Indeed, holding of inequality (10) solves the problem posed, since the greater is the value of the coefficient 
of correlation, the less will be the error of the leading link and, consequently, the higher will be the accuracy of 
control of the system to be designed. 


2. Use is made of standard universal controlling devices. With optimal tuning of the selected regulator 
which is implied here), the mean-square error of the controlled quantity will, essentially, be determined by the 
average variation of the controlled quantity during a time equal to the lag time, since the regulator action will 
lag the changes in the controlled quantity in the object by time r. 


The mean-square variation of the controlled quantity during time r, L;, can be computed from the values 
of the correlation function of the given quantity, Indeed, by definition, 


Le VM +t) —2 (OP) = VM (t+ 2) — 22 +27) = (11) 
4- M - - 2M [x (t) x(t + 


where M is the mathematical expectation, 


By taking into account that 
(O) = M [2? M (t + (12) 
Rx = M [x (13) 


and by substituting in Equation (11) the values of the mathematical expectations given in expressions (12) and (13), 
we obtain, finally, 


L, = V (0) — (14) 


Thus, to have a bounded average variation of the controlled quantity time 1, the following condition must 
hold: 


1 
V2z[R, (0) > Ps. (15) 


In the problems considered, conditions (10) and (15), for guaranteeing a given quality, are concrete forms 
of the more general expression (9), Thus, if condition (9) is met, one is guaranteed that the operation of the con- 
trol system will provide a given average accuracy with optimal tuning of the chosen regulating device, Con- 
sequently, inequality (9) may serve as the condition for guaranteeing the quality of the control system to be de- 


signed, 


We now consider the design of control schemes for objects with lags if the given condition for guaranteed 
quality (and, all the more so, the existence condition) does not hold, In order to get the possibility of creating 
a control system of given quality, it is necessary either to decrease the lags in the object or to change the form 
of the correlation function of the quantity being considered, Decreasing the lag is, in practice, unfeasible; it is, 
however, possible to propose several methods toobtain the necessary changes in the form of the correlation func- 
tion, For the holding of the condition, guaranteeing quality [ef., its concrete forms in (10) and (15) ] for a given 
T, the rate of decay of the correlation function should be decreased, This requirement may be met both by sta- 
bilizing the highest-frequency disturbances whose amplitudes are significant in the given process, and by control- 
ling some intermediate quantity in the object which is stochastically related to the output quantity being consid- 
ered, which has a significantly smaller delay time for which, in its turn, the existence condition holds,* 


* In case one finds an intermediate quantity v(t) which is related, not stochastically, but functionally to the out- 
put quantity x(t),ie., v = f(x), then it can itself be taken as the output quantity and, consequently, the problem 
of controlling the quantity x(t) disappears in practice, 
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Control schemes for quantities for which the condition of 
rity guaranteed quality do not hold, i.e., for which it is impossible 
i) . to construct control systems of the structure shown in Fig. 1, 
ordinarily consist of at least two control systems: a basic system 
and a subsidiary one, In the overwhelming majority of the cases, 
such systems have the form shown in Fig. 2. In variants a and b, 
the subsidiary control system stabilizes a chosen disturbing fac- 
tor, m(t) or n(t), In variant c, the subsidiary system controls the 
’ intermediate quantity v(t). In variant a, the subsidiary control- 
Fig. 1. ling device acts on an individual controlling organ, while in 
variants b and c there is a two-stage system, which permits two 
quantities to act on one controlling organ, The finding of the disturbing factor or intermediate quantity which 
requires stabilization, as with the choosing of the corresponding scheme, is carried out on the basis of a concrete 
analysis of the operation of the object to be automated, 


W,ip) 


The design of the subsidiary control system is carried out with the purpose of lowering the frequency and 
speed of variation of the output parameter. This purpose also determines the method of estimating the quality 
of operation of the designed subsidiary system, namely, the degree of variation of the form of the correlation func- 
tion of the output quantity, Indeed, as a result of controlling a definite disturbing factor or intermediate quantity 
it is necessary that the rate of decay of the output quantity's correlation function should be decreased enough to 
render it possible to construct a control system of the desired quality for it. 


. b miti c mt} 
a(t) ft} 
v(t) 
2 2 2 2 
2 


Fig. 2. 1) Object; 2) controller. 


Consequently, the degree to which the condition for guaranteed quality, (9), of the output quantity holds, 
given the operation of the subsidiary control system, 


(0), Rz(t)} —P=N (16) 


is an estimate of the operation of the designed subsidiary control system, The estimate just cited allows one to 
determine the correctness of the choice of a concrete intermediate quantity or disturbing factor to be regulated, 
This estimate can also be used as a measure of the accuracy of operation of the subsidiary control system. 


If, for one object, several variations of subsidiary control system have been constructed, these systems dif- 
fering by the choice of controlled quantity or by the type of controlling device, then comparisons of their opera- 
tion and the choice of the best variation must also be made by means of the estimate cited. If the design pecull- 
arities of the constructed variants are not taken into account (and these may have an essential influence on the 
choice), then the best subsidiary control system should be that one for whose operation the quantity N in expres- 
sion (16) has the greatest value. 


If, with optimal operation of the chosen subsidiary control system, the quantity N in expression (16) is nega- 
tive, f.e., the condition of guaranteed quality of the output quantity still does not hold, then one must consider the 
question of constructing one, or several, subsidiary systems for controlling other disturbing factors or intermediate 
quantities, In this case, one must continue operating with sets of several subsidiary systems until, in estimate (16) 
the quantity N goes positive. 
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The method just described for designing control systems for objects with lags was employed for developing 
an automatic,control scheme for rotating cement-roasting kilns, Without going into all the questions connected 
with the automation of kilns (cf., for example [4]), we shall consider the design of a scheme for the basic out- 
put quantity— the qualitative indicator of kiln operation — the weight of a liter of clinker. The weight of a liter 
of clinker (WLC) is determined automatically asthe material is discharged from the rotating kiln. The basic factor 
affecting the quantity considered is the supply of fuel. Variation of the fuel supply leads to a-variation of the 
temperature field of the material along the kiln, which influences the course of the chemical reactions and, con- 
sequently, the quality of the final product, The transfer function of the object with respect to the quantity con- 
sidered is approximated by a single~capacitance link with a lag: 


Twi p+i 


Kins (17) 


where X,,, ~ is the value of the output controlled quantity in grams, and Xyp is the value of the controlling quan- 
tity in percent. 


For the concrete object considered, kw; = 9 grams, Tw c = 25 minutes and Twic ='10 minutes. (Here 
and in the sequel, all experimental data are given for a No, 2 rotating kiln of the Leningrad Cement Factory.) 


As the regulator of the quantity considered, we chose a universal intermittent controlling device. In ac- 
cordance with the technological requirements on clinker quality, the mean-square variation of the quantity (ex- 
cluding the regulator's dead zone) must lie within the limits + 20 grams, Thus, in accordance with (15), the con- 
dition of guaranteed quality for the control system to be created has the form: 


i -1 
> 
O— 


(18) 


Determination, from experimental recordings, of the correlation function* of the WLC (Cf,, the curve for 
Rwci(t) ou Fig. 3) showed that not only was the condition of guaranteed quality, (18), not met, but that condi- 
tion (7) also held, i.e,, the condition for existence was not met, (Indeed for Tw ;c = 25 minutes, we find from 
the curve of RwWLC;(t) that t) = 20 minutes.) Experiments reiterated the impossibility of controlling the output 
quantity, 


In accordance with the method we have presented, in order to design a system for WLC control, it is nec- 
essary to implement a subsidiary control system, Analysis of rotating kiln operation showed that the basic dis- 
turbing factors are the variations in the chemical and physical properties of the raw material charged into the 
kiln. Until the independent problem of automating the preparation of the raw material in preparatory crackles 
(slurry preparation) is solved, it is not possible to carry out rigid stabilization of the fundamental disturbing fac- 
tors, An experimental study of the various quantities characterizing the process in a rotating kiln showed that 
one of these quantities — the total temperature in the clinkering zone (the hottest zone in the kiln), determined 
by a radiation pyrometer sighted from the kiln head — can be an intermediate quantity, In fact, this temperature 
is stochastically related to the output quantity, the WLC (their coefficient of correlation is approximately 0,5). 
At the same time, the delay in variations of this temperature after variations in the fuel supply is practically non- 
existent, and the transfer function is well approximated by a single-capacitance link with a time constant on the 
order of three minutes, 


After construction of a control system for the temperature in the clinkering zone, which acted on the fuel 
supply mechanism (as a regulating device we used an isodrome universal VTI controller of type ER-III), the 


* The WLC correlation functions were determined from recordings of actual industrial aggregates, which did not 
permit the carrying out of the experiments in a pure form, All the variations of the WLC on the recording tape 
lie within comparatively narrow limits, This is occasioned by the operator's manual intervention in the process 
when the output's WLC goes beyond the limits of normal quality. In view of this, there were certain distortions 
in the form of the correlation function, and a noncorrespondence of ¥ Rwy with the mean-square error sus- 
tained by the WLC, which would have occurred in the absence of the manual intervention, 
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Fig. 3. Correlation functions of the kiln's output quantity, 
the WLC, for manual operation of the aggregate, Rw) c1(t), 
and for temperature regulation in the clinkering zone, 


RwLc2(t). 


correlation function of the output parameter, Rwi.C2(t), took the form shown in Fig. 3. In this case, the condition 
for existence, (8), held, establishing the possibility of constructing a stable control system (ty ~ 37 minutes), As 
shown by a substitution of the experimental data in formula (18), it is possible to construct a system with a given 
average accuracy, Indeed, RwLC2(0) = 227 gram? and RwLC2(25) = 86 gram’; hence, 


V 21 o(26)] = 17g. 


The complete control scheme for the WLC, as designed, corresponds to Fig. 2, c. 


Testing of the given control scheme on aggregates showed that it provided the given average control ac- 
curacy, Thus, during 48 hours of operation, about 75% of the time the controlled quantity was found within the 
limits given by the sum of the controller's dead zone and the average variation of the quantity during time r 
(17 grams), and only for about 25% of the time was the controlled quantity somewhat beyond these limits. 


SUMMARY 


1, In creating a control system for an object with lags acted upon by random disturbances, one should com- 
pare the decay time of the correlation function of the quantity to be controlled (ty) with the duration of the trans- 
port lag(r) in the object. The inequality t, > r is the existence condition, since its holding determines the pos- 
sibility of constructing a control system for the quantity under consideration, 


2. For the construction of a control system with a given average accuracy, it is necessary that the condi- 
tion of guaranteed quality, (9), hold. 


3, If the condition of guaranteed quality (and, even more so, the existence condition) does not hold, it is 
necessary to construct a more complex control system in which, together with the basic control system, there is 4 
subsidiary system controlling the highest- frequency disturbing factors with amplitudes which are significant for the 
given process, or intermediate quantities of the object, 


4, The degree to which the condition of guaranteed quality for the basic controlled quantity is met, given 
operation of the subsidiary system, 


N =F {Ry (0). Re 


is the measure estimating the accuracy of operation of the constructed subsidiary control system. 
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5. Using the method presented, we gave the construction of a control system for one of the output quan- 
tities of a rotating cement-roasting kiln, for which the existence condition was not met. The operating results 
of the control system constructed on industrial aggregates support the advantages of using the method presented 
here for designing control systems for objects with pure (transport) lags. 
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STEADY-STATE PROCESSES IN THE SIMPLEST DISCRETE 
EXTREMAL SYSTEM WITH RANDOM NOISE PRESENT 


A. A. Fel'dbaum 
(Moscow) 


Investigation of the processes in discrete extremal systems with random noise 
present , leads to the consideration of equivalent Markoff chains, Based on this con- 
sideration, we determine the steady-state error and the optimal size of the step (the 
quantum size) in the simplest system. Certain generalizations of the problem posed 
are considered, 


Posing of the Problem 


One of the most important classes of self-adjusting systems is that of automatic optimization systems, The 
task of these systems consists in searching automatically for the extrema of functions of, generally speaking, many 
variables, the set of whose values is limited to some definite domain. In the simplest case, we consider a search 
for an absolute extremal without taking limitations into account, Such systems are called extremal systems, The 
simplest extremal system finds the absolute extremal of a function Y of one variable X, with the condition that the 
corresponding value of X is not on the boundary of the range of variation of X. 


uliieetientientianstiontiontintibns 4 In automatic optimization systems, and in particular, in ex- 
| teh | tremal systems, random factors play a significant role. In spite of 
1 Zola - g oy this, there are still few works which take these factors into account, 
! , / ' In [1] there is considered the probability of a false step using the 
! ' correlation method of measurement. Work [2] gives a method for 
designing extremal systems, with an estimate of the tracking error 
win) C form in the absence of random noise, and the probability of a false step 
i nee with this noise ‘present. Those parameters are chosen for which the 


probability of a false step would be minimal for a given tracking 
error, In [3] there is considered the influence of random noise on 
Fig. 1 search time in the transient response for several types of simplest 
discrete extremal systems, 


The task of this work is the investigation of the steady-state processes in the simplest discrete extremal sys- 
tem, the block schematic of which is shown in Fig. 1. All quantities in this system are quantized in time, i.e., 
they are considered only at discrete moments of time t = Ton, where n = 0, 1,2,3,... By means of a linear trans- 
formation of the relationship between the input, X[n], and output, Y[n], of the object 0, we can go over to the 
relative quantities x{n] and y(n}, where the minimum of the curve y = f (x) will coincide with the origin of co- 
ordinates, ‘Je first consider the case when x(n] and y(n] are related by the function 


y(n] (1) 


This functional relationship is shown in Fig. 2, In addition, at the object's input there acts the input noise 
Ze{n] and, at its output, output noise z[n}, The system's task consists in obtaining the minimal value of y(n], whic 
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is attained with the condition that x[n] = 0, or u{n] = There- 
fore, the controlling stimulus u[{n] must “follow” the noise z,[n] and 
compensate it. The effect of the noise z9{n] generally consists of 

a comparatively slow translation of the curve y = f(u) along the axis 
of abscissas. We shall now assume thatzo[n] = 0. 


We shall consider as subject to measurement only those quan- 
tities lying outside of the dashed rectangle on Fig. 1, the block sche- 
matic of the object controlled as a whole, At the output of the ob- 
ject, there is added to the quantity y(n] the noise z[n}. Thus, the 
measured quantity is 


= y[n} + (2) 


For different values of n, the noise increments, z[n + 1] —z[n], will be considered as statistically independ- 
ent random variables with zero mean and given symmetric distribution laws, The quantity w[n] is applied to 
controlling device C, whose output, u{n], acts on the controlled object, 


Controlling device C acts in the following manner, Let there be some initial value, u = ug, with u = ska, 
and k = 0,1, 2, 3,..., and a= const, The system makes the initial trial step Au = +a, and the value of u becomes 
equal to us +a; with this, the quantity w = w, is measured. Then a second trial step, Au = —a, is made, and the 
value of u becomes uy — a; now the quantity w = ws is measured, The difference, Aw = wy — Wy, is used for the 
determination of the operational step in accordance with the formula 


Au [n] = — asign Aw [n]. (3) 


We shall agree that each measurement of the quantity w requires some sub-interval of time, but that the 


carrying out of a step requires no time, Then we shall understand by the timeT, the time required to make the 
two trial steps, This time may be taken equal to unity. 


By complicating the logic of the controlling portion, one can make repeated measurements of the quan- 
tities w, and we, i.e., m pairs of trial steps with subsequent averaging of the results of measuring Aw. This de- 
creases the probability of making a false step, i.e., a step in the wrong direction caused by the presence of the 
noise z{n}. This probability is the probability of the event that the sign of Aw[n] does not coincide with the 
sign of Ay(n}. With u # 0, the probability of a false step for the curve shown in Fig. 2 is everywhere constant. 
We denote it by the letter p, and the probability of a correct step by the letter q; then 


q=1-p, (4) 


For up = 0, the trial steps + a give identical values of y; consequently, Ay = 0 and the probability of a false 
step equals 0.5, 


In the steady state, the values of u(n], x(n], and y(n) vary by jumps, and the average value of y(n] is different 
from zero, i.e., errors appear. The task of this work is to determine the average value of the error for the quan- 
tity y[n} and the optimal value of the step a, for which this error is minimal, 


System Equations and Method of Investigation 


The process in the system being considered is a discrete Markoff process or, as it is sometimes called, a dis- 
crete Markoff chain, This is a process in a system with a finite or countable number of states and discrete time, 
t=0,1,2,..., n, where the probability, pj(n], of finding the system in state i at time n depends only on the pro~- 
babilities, pjlo-1), of the different states at the previous moment of time, n- 1, 


We denote by p,[n] the probability of the value u[n] = ia. If the set of possible values of i is infinite, then 
i=—o, ...,—-2,—1, 0,41, +2, ...., #0. A transition to the valueu[n] = ia, as is obvious from the description 
of system operation, is possible only from the values (i — 1)a or (i + 1)a, the probabilities of which, for time (n—1), 
equal p;_i{n—1] and 1). Therefore, 
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where Pj-;, { is the probability of a transition from (i—1)a toia and Pj 4;, { is the probability of a transition from 
(i + 1)a to ia, Since the movement to a minimum has transition probability equal to q, and the movement from 


the minimum has transition probability p, the probability of state pj{n] may be computed by formulas deriving 
from the control method adopted: 


p_s{n] = —1)q + 
= p.sin —1]q+piln— 4] p, 


pain] 
Po {nj Piln (6) 
= —1) > + paln—1l 9, 


Po(n) = pi (nm — 1) p +- — 1] q, 
Ps = pa {nm — 1] p+ 1) q, 


Thus, the process in the closed system being considered is a discrete Markoff chain, Since, at any moment 
of time, the system lies in one of its possible states, we must add the following condition to Equations (6): 


>) Pin) =1. (7) 
i=—oo 


If the initial probability distribution of pj[O] is given (as will be obvious in the sequel, this is not required 
for the solution of the problem posed above), then pj[n) may be found for any moment of time n, and we may 
compute the mathematical expectation of the quantity y: 


My ly] ++ + 2ap_y [n] + ap_, + 
+ apy + ap; [n] + 2ap, [n] + [n] +--- (8) 


The coefficient of pj[n] is called the weight of the corresponding state, For states with i # 0, the weight 
equals the arithmetic average of the values of y for two trials, i.e., |ija. For i = 0, the weight equals a, since 
both trials give the value y = a, 


If n + ow, we obtain the limiting value M,,[y], if it exists, in the steady state. It follows from formula (8) 


that 
24 = 3p_s [00] + [oo] + 1 [oo] + 
+ 1 py [oo] + 1p; [oo] + 2p, [oo] +- [oo] + 
where 
P, [9°] = lim p; [n). (10) 


It is possible to take 2n as the measure of the error in the steady state, Other criteria for the steady-state 
error are also possible (for example, the probability of the transition of y over some threshold €), However, 
the criterion given above has practical meaning, and, as will be shown below, makes it possible to get 
comparatively simple computational formulas, 


It may happen that the probability p(n), in the steady state, varies periodically, or quasi-periodically. 
Then, the limits in (9) and (10) do not exist. In such cases, one may use for the criterion the time average of 
the steady-state values of the mathematical expectation of y. The corresponding formula has the form: 
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‘be 
= lim >; Mall} (11) 
n=) 


If the limit p;[oo] exists, formula (11) converges to (9). 


The problem stated above is statistical by nature, However, after the way in which it was formulated, the 
problem's solution leads to the investigation of a finite, or infinite, set of linear. finite-difference equations, (6), 
Equations in finite differences are studied in the literature of Markoff processes [4-7], in works on the theory of 
finite-difference equations [8, 9] and also in works on the theory of sampled-data systems [10, 11]. For system 
investigation, it is useful to construct an equivalent schema of the discrete Markoff process (Cf., for example, [12]). 
We sketch the schema corresponding to Equations (6) (Fig. 3,a), The thick points on it denote the various states 
of the schema, corresponding to the various values of u = ia (i= —o, ..., ~2, —1, 0, +1, +2,..., +0). The 
transition paths from one state to another are shown by arrows with the given probabilities of transition. The 
weight of each state is shown above it within square brackets, 


For investigation of Equations (6), we may forget about the manner of production of the quantities pj[n] and 
simply consider them as the values of the signals at the i'th point of the scheme in Fig. 3, at the moment of 
time t=n. Then in Equations(6), the coefficients of pj[n —1] may be considered as the gains of the corresponding 
portions of the system. Since, in these equations, pj[n] depends on the values of pin ~1)(j = i-1, 1+ 1), each 
portion of the scheme must contain an element with a delay of one interval, As the result of such a replacement, 
we obtain the equivalent schema, shown in Fig. 3,b, The operator z in each link symbolizes a single-interval 
delay. It is easily seen that the process in the discrete schema of Fig. 3,b is characterized by Equations (6), Giv- 


ing the initial distribution of the pj[0] may be implemented by, for example, pulses applied at the points A,, these 
pulses satisfying condition (7). 


The method of Z-transformation is frequently used in the investigation of processes in discrete systems. 
The possibility of using this method was looked into for solving the problem posed here, However, very compli- 
cated expressions were thereby obtained. 


It is possible to obtain the limits sought by a comparatively simple route, without going to transformations, 
We consider any one of Equations (6), If the limit exists 


r, = lim p, [n] = p, [oo], (12) 
then by letting n go to infinity, we obtain, instead of the equation, for example for i = 2 


Pp, = p,_,(n—1) p +P, (13) 


the following equation: 


(14) 


As a result, we obtain, in the general case, an infinite system of equations in a infinite number of unknowns 
rj. With this, the expression Mg)[y], the sum of the rj multiplied by their weights, is meaningful if each of the rj 
tends to zero as the number of equations tends to infinity, 


If the limit in (12) does not exist, an equation of the type of (14) can be obtained by adding Equation (13) 
forn=1,2,... , N, dividing by N and letting N tend to infinity, In this case, ry and Mgo[y] are average values 
in the sense of definition (11). It should be mentioned that, in the case considexed, the limit in (12) exists for 
p> 


Steady-State Error and Optimal Step 


We consider initially a system with a finite number (2k + 1) of states (k > 2) and with equations 
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= p, (mn — 1) p+ p,_, (nm — 1] p. 


Fig. 3 


Here, the extreme values, u, and u.,, in case of a false trial outcome, remain invariable. The schema of 
this process is shown in Fig. 4. 


Fig, 4 
By adding the right and left members of the equations termwise, we find that 
im+k i=+k 
> Pp, (n) > P,{n — 4}. (16) 
i=—k 
Since 5) p,[0] = 1, it follows that 
p,(n] = 1. (17) 
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b gz q L 


(15) 
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(16) 


(17) 


If Equations (15) and (17) are considered together, then one of Equations (15) — for example, the last one — 
is redundant, 

By replacing the p;[n] by the limiting values, rj, we find the limiting system of equations, obtained in the 
same way as Equation (14) was derived from (13): 


Py = 


1 


To (18) 


4 
Te + rap, 


The last equation, for ry, 1s discarded here, Instead of (17) we find that 


i=+k 
1. (19) 


t=—k 
In Equations (18) and (19), the subscript i may be replaced by —1 without changing the system, Consequently, 


(20) 
We may therefore limit ourselves to the abridged system of equations: 
r, = 2r,9, 
i 
(21) 
To 
= "eI 
with the added condition that 
k 
i=] 
It is necessary to find the magnitude of 
) 
int thre = + Sir, (23) 
The relative error sought is 
M. 
co [yl (24) 


Px (p) 2k 


We may find the unknowns r; from Equations (21) and (22). Initially, we substitute the value of rg from (22) 
in the first equation of (21), We also eliminate rp from the first two equations of (21), Then the system of equa- 


tions takes the form: 
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Te ="sQ+rip, 
+ Tap, 


(25) 


From the second equation of (25) we find that 


(26) 


where B, = p/q. From the third equation we find that 


P 
2 
=r,—! =7,(4) = (27) 


where By = (p/q)*. Generally, if we have the equation 


(28) 
then 
i 


For v = 2, 3,...,j—1, let the equation 


be valid. Then 
Thus, 
r, =», (32) 


where X = p/q. By substituting expressions (32) in formula (23), and taking into account the first equation of (21), 
we find that . 


(a+ i?) r(a+ se), (33) 
j=1 j=1 


After substitution of (32), we find from the first equation of (25) that 


ry = (34) 


k 
a+ 
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(26) 


(27) 


(28) 


(29) 


(30) 


(31) 


(32) 


yn of (21), 


(33) 


(34) 


It follows from this that 


’ (35) 
where 


We now let k tend to infinity. In the limit, 


q+ 
(p) = (p) = 


(37) 
Since 
co 4 4 
j=1 j=1 
then 
g(i—az+i 
= (39) 
If we substitute the value of then, after some transformation, we arrive at the formula 
e(p)=14+5 (40) 
This is the basic design formula, The error sought is computed from the expression 
Moly] = (p) 4. (41) 


The magnitude of the probability, p, of a false step is itself a function of a, Since the magnitude of a trial 
step, i.e., the difference between the values of u + a and u— a, equals 2a, then 


—2a 


—co 


where F is a function of the probability distribution of the noise and W(x) is the density of this probability dis- 
tribution, For example, for a Gaussian distribution with dispersion 07, the probability of a false step is 


Then, 
oo [y] ag (p) (44) 
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Now we introduce a relative step, defining it as 


2a 
o = —. (45) 


With this 


M 
M,= co 


1 ®(— 9) 
= 9{1 + (46) 
The curve for My = f(#) can be constructed from this formula. Analysis of the formula shows that the rela- 
tive error, Mo, increases with increasing relative step 9. Consequently, in the given case, the error is the smaller, 


the smaller the step is (Cf,, below, as well), For 9 -+ 0, the minimal error, as may be computed by formula (46) 
(for this, @[—9] + 0.5), equals 


(Mo) 0125+ (47) 


If repeated measurements of the quantity Aw are made, then for m repeated measurements, the probability 
of a false step equals, as is well known, 


pi(m, 9] = (48) 


It is easily found from (46) and (48) that, for $— 0, the minimal error is 


49 
(Mo) in (49) 


’ Thus, in the given case, (a),»; = 0, and the error decreases monotonically with increasing number of re- 
peated trials, As should be expected, this repetition is desirable in the steady-state, But the optimal number of 
repetitions is infinitely large only in the case when the noise at the input, z9{n], equals zero, If there is a shift 
of the curve y =f (x) along the axis of abscissas then, with too large a number of repetitions, tracking will be 
carried out with a large error, Consequently, if we take into account that z9[n] # 0, the optimal number of re- 
petitions must be finite, 


We now consider the case when larger errors are more probable than smaller ones, and the probability den- 
sity, W(x), for the error has the form of the curve given in Fig. 5, a, The function F(x) is shown in Fig. 5,b. 


The probability of a false step is shown on Fig. 6 for several values of y, where 


| (50) 


2a P (+) 


If we sabstitute the linear relationship p = 0.5 — y_ », we may easily show that dg /dy > 0, i.e., the error 
increases as the step increases, This 0 and d*p /du* > 0, as occurs, for 
example, with a Gaussian distribution, 
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7 
lz on 
u=0 
and 
(51) 
The quantity A is shown on Fig, 5,a, Formula (41) for the error may be written here in the form 


45) 


(46) 


aller, 
46) 


(47) 


tlity 


(48) 


(50) 


(51) 


By formula (52), for the curves p(#) shown in Fig, 6, we con- 
Wis} structed the graphs € = &(#), (Fig. 7). It is clear from these graphs 
that, for small values of y 3, the size of the optimal step will be greater 
2 than zero, For example, for y, = 4/16, the magnitude of pope = 

= opr / 2A 0.85. As yg increases, the size of the optimal step de- 

creases, Further, for y, = 0.3, the errors for p = 0.7 and p = 0 are 
identical and minimal, and for y, > 0.3, the minimal error corres- 
ponds to the value p = 0, i.e., the optimal step equals zero, 


We now explain the physical meaning of the results obtained, 
An increase of the step, a, has a two-fold effect on the error. On the 
one hand, an increase of the step increases the hunting about the mi- 
nimum, thanks to which the error must increase, On the other hand, 
an increase in the step decreases the probability of a false step, which 
is capable of decreasing the mathematical expectation of the quan- 
tity, y, i.e., of decreasing the error, For large values of the step 4, 
Fig, 5 the second factor plays no role; it is therefore obvious that too largea 
step is harmful and must be decreased. However, as the step is de- 
creased, both factors begin to have effect, acting against one another, Only computation allows us to know which 
of them will predominate. In the case when the probability of large noise is significantly greater than the proba- 
bility of little noise (for example, for y = 1/16 on Fig. 6), too small a step “drown#*in the noise. Computations 
show that, in this case the optimal step must not be too small, 


a2 


a2 as 6 a7 typ 
Fig. 6 


If dW/dx < 0 for x > 0, i.e., a larger noise is less probable than a smaller noise then, as was found above, 
the tendency for the error to decrease as the step decreases predominates, Therefore, the optimal step may be as 
small as desired. 


The results just cited are valid only for zp{n} = 0. If the input noise differs from zero, then too small a step 
increases the error, For example, if ze{n] is:a monotonic function then, for too small a step, the magnitude of un} 
cannot follow the variations of ze[n] and the error increases, Therefore, even for small zo[n} = 0, the curves of 
Fig. 7 must be changed, and for small p the magnitude ofdg /dp must be negative, 


Certain Generalizations 


The problem considered above, and the methods of solving it, can be generalized in various directions: to 
large numbers of input variables, xi, (i= 1,2, ... , 0), to other forms of logic of the relationship y =f (x), to other 
forms of logic of the controlling device C, to other forms of noise and other probabilistic characteristics and method 
of application of the noise, Finally, additional conditions and limitations may ensue, for example, limitatious on 
the domain of variation of the variables x; to a certain surface in the n-dimensional space of these variables, 
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For example, let the output quantity, y, of object 0 depend on n input quantities, x;, in the following way: 


y= Dalal, (53) 
i=1 


where zoi[n] = 0, z[n ] # 0, and the trial steps, aj are carried out separately and in turn for each of the variables, 
after which the operational steps are carried out simultaneously. Then, the probability pot a false step, with 
respect to each of the variables xj, will differ from variable to variable: 


p = p® (2a,a,). (54) 


In this case, the steady-state error is defined by the expression 


n 
= >) (Pp), (55) 
i=1 


where the expression for ¢(p) is somewhat different from (40), due to the changes in the weights, since a trial is 
made during two cycles and then, during 2(n — 1) cycles, the quantity x; remains invariant, For example, for the 
state i= 0, the weight will equal 


1+4+4+2(n—1)0 


i 
n’ 
It 1s also possible to consider other types of the relationship y = f(x), for example 


y [n] = 2* (56) 


In this case, for the trial steps + along the axis of abscissas and on both sides of the values uj = ia, the quan- 


tity y changes by 
Ay = 4ia’, 


Therefore, the probabilities of false and correct steps will depend on i: 
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(56) 


> quan- 


—4 a? 
p= p = pi) =F = W (2) dz. (58) 


The corresponding schema of the Markoff process in which the transition probabilities, formerly constant, 
are now functions of the state's ordinal number, is shown in Fig, 8, for i= 0, p = 0.5, It should also be borne in 
mind that, in this case, the weight of the i'th state equals a%(i? + 1). 


Fig. 8 


Here, the method of investigation presented above is applicable, but the computations are significantly 
complicated. Let z»{n] be a regular function, varying linearly as a function of n, If, during one time interval, 
the translation of the minimum of the curve u = f(u) equals a/m, where m = 2, 3, 4, ... , then the analysis may 
be brought to the investigation of a system of finite-difference equations with constant coefficients. 


For simplicity, we consider the case of m = 2, where it is assumed that the function z{n] remains invariant 
during the trial steps (otherwise the computations of the state weights and the probabilities of false steps would be 
complicated), and then is changed by +a/ 2 before the following trial steps. In this case, it is convenient to begin 
with the different possible values of the quantity x: 


[n] = zo[n] + u[n}. (59) 


Since zo[n] = ak/ 2, where k is an integer, the possible values of x are also expressed by the formula ak/ 2, 
We shall number these values ...., ~2,— 3/2, —1, —1/2, 0, +1/2, +1, + 3/2, +2, ... as shown in Fig. 9. This 
figure shows the equivalent schema of the corresponding Markoff process, If the function z9{n] were equal to zero, 
then the magnitude of x(n + 1] would be different from x[n] either by +a or by —a, depending on the initial trial, 
However, due to the variation of zg[n + 1] by a/ 2 as compared with zo[n], the quantity x[n + 1] experiences an 
additional increase of a/2. Therefore, x{n + 1] differs from x[{n] either by 3a/2 or by —a/2, which is shown on 
the schema of Fig. 9. 


Here the probabilities py and qp for i = 41/2 differ from p and q, since 
Po= F [—a], (60) 


while p = F[-2a]. 


In the general case of an arbitrary regular function zo{n), the problem reduces to the investigation of a sys- 
tem of linear difference equations with variable coefficients. If zg{n] is a periodic function of n, with period T,, 
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then the system of equations obtained has periodic coefficients. The solution of such a system, generally speak- 
ing, can be found (Cf., for example, [9,11]. The steady-state solution is also a periodic function of n, In the 
given case, the weight y,(n] of state p(n] is a periodic function with the same period, T,, as the function zon]. 
The desired mean error may be computed by formula (11): 


N ,i==+co i=+oo T; 
= 2y = lim in} » a PP (61) 


Here, py[n] is the state probability in the steady-state, The right member of Equation (61) is written for the 
case when these probabilities are periodic functions with period T, (generally speaking, they may have periods 
which are multiples of T,). 


For definiteness, we consider the simple case when 


0 


[a] = T 
—a, n= Ti 


and the period equals T;, i.e., function (62) repeats periodically, Then, for the first half-period, the equivalent 
Markoff process schema has the same form as in Fig. 3 (Fig. 10,a). For the second half-period, the image of the 
process is shifted by one step to the right (Fig. 10,b), Omitting the superscript oo, we consider, for the steady- 
state process, any of Equations (6) or (15), for example, 


Pl) = P44, p+ —2). (63) 


We write this equation for n= 1, ... , T;/ 2, and then add the equations obtained. As the result, we find 
that 


Ts TM 4 
2 2 2 
> > (64) 
n=] nr n=0 
or | 
2 2 T 
>) Pill = Pigs + (0) — + 
n=1 n=] 
T 
> + (0) — . (65) 
We introduce the notation a 
T 
n=] 


Then, from expression (65) we find that 


Di=P Dig tPA (67) 


An analogous equation can be obtained for the second half-period: 


War te (68) 
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(64) 


(65) 


(66) 


(67) 


(68) 


where p* and q‘ are the corresponding probabilities for the second half-period (they are not equal to p and q for 
all states) and 


(69) 


Fig. 10 


Thus, for >}, and }), we obtain equations of the same type as Equation (18) for 4, but with added terms 
A;. Moreover, equations of the type of (19) are valid for 5}, and >); , specifically 


By knowing Aj we may find 5), and 5);, and from there we can find the error from expression (61), The 


values of Aj can be found from Equations (15) for the first and second half-periods. We write these equations in 
matrix form: 


pla] =Pp[n—1], p’ = P’ p’ [»—1)]. (71) 


Here P and P" are the transition probability matrices for the first and second half-periods respectively, and 
p(n] and p*[n] are vectors with coordinates pj[n) and p;‘[n] respectively. 


We find from Equations (71) that 
If we now add the periodicity conditions 
(13) 
we can find p[0] from Equations (72) and (73), by taking condition (17) into account, and then obtain p[T, /2] and 
4;. By using the matrix notation we can find solutions in general form: however, the concrete calculations are 


attended with horrendous computation, 
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ON APPROXIMATE SYNTHESIS OF OPTIMUM LINEAR 
DETECTING SYSTEMS 


Yu. P, Leonov 


(Moscow) 


The paper deals with the approximate solution of the integral equation for 
the weighting function of the optimum system. 


1, We consider the following task, We are given a signal y(t) = g(t) + n(t), where g(t) has the form: 
g(t) = Sache), (1) 


the a, are unknown but the #;(t) are given functions, and n(t) is a random function with known correlation func- 
tions Rp(ty, ts), 0 < ty, tg < @, where M[n(t)]= 0. It is necessary to determine the weight function [1], k(t,r), of 
the linear filter at whose output, with the signal y(t) applied to its input, will be obtained an estimate, g *(t), of 
the mathematical expectation* * of the random function y(t), satisfying the two conditions; 


a) M {g*(t)}=g(t), t>T; 
t 
b) Dy, = M {n, (t)}* = min, where n, (t) = \ k(t, t) n(t) de. 
t—T 
Thus, we are talking of obtaining the best, in the sense of minimum dispersion (condition b), unbiased. 
estimate (condition a) for the class of linear filters (in general case, with variable parameters), 


If we consider the class of input signals in (1), and remove, ina the condition that t> T then, by writing 
the signal at the system's output in the form * * * 


t 
(eae, 


we obtain, on the basis of a the following system of moment equations: 


* This work was presented at a seminar on the applications of probabilistic methods to automatic control, held 
at the IAT AN SSSR on February 17, 1958, 

** The problem is easily extended to the case when it is necessary to obtain an estimate, not of the function 
g(t) itself, but of an arbitrary linear operator on g(t) [1]. 

*** This formula is valid for the case when k(t, 7) satisfies the conditions of physical realizability k(t, r) = 0 
fort< r and, in addition, g(t) = 0 fort < 0. 


Mos - 
0 
| 
i 
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However, the class of signals in (1) is very limited, It can be exapended if we introduce systems with finite 
memories into the consideration. 


If one takes into account the conditions of physical realizability and the finite memory, then the weight 
function of such systems will satisfy the following conditions: 


k(t,t)=0 foro O<St<+00, tctrandt<t—T (4a) 


k(t,t)4-0 for O<t<+o0, (4b) 


where T is called the system's memory. 


On the basis of (2), the signal at the output of such a system is written as 


t 


0 
t 

gi(t)= \ k(t (edt, 6) 
‘—T 


It is clear from formula (6) that if one adds the requirement that t > T to the condition of unbiasedness, then 
the representation of the signal in the form of.(1) for a system with a finite memory is only required on the memory 
interval (t — T,t) and not for all t, as it is in the case of a system with an infinite memory, as in (2). 


Thus, in the case of a system with a finite memory, it is required to know only the local behavior of the 
function g(t) on the interval of the system's memory. The function g(t) need not be represented in the form of (1) 
for the entire interval of time considered, In this consists the possibility of broadening the class of input signals 
to include the case of a system with a finite memory, or the possibility of decreasing the a prioridata on the input 
signal g(t). Precisely for this reason was the requirement t > T introduced in condition a) and, for obtaining un- 
biased estimates, we consider systems with finite memories. 


For the case of a stationary function n(t), when R,(t;, tg) = Rn (ty — tg) = Ry (7), the problem is solved, and 
an expression for k(t, T) is found in closed form* [1, 2]. For the case of a non-stationary function n(t), a solu- 
tion in closed form is obtained [1] for the class of processes n(t) which satisfy a linear differential equation with 
variable coefficients, of the form 


= (t), (1) 
where 
Mm = bm (2), 
Lt = 0, + +++ ai(t), 
and the coefficients ag(t), ax(t), ... An(t), bolt), by(t), ... p(t), are given, 


* In this case, the optimal system for a sufficiently broad class of problems is described by a system of linear 
differential equations with constant coefficients, Consequently, the weight function, k(t, rT), depends on the one 
variable, (t- T). 
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Let €(t) be a random process with independent increments, It is obvious with this that € ‘(t) is stationary 
"white noise", i.e., the correlation function and spectral density of the process € ‘(t) are defined, respectively, 
by the formulas R¢ *(t) = c? 6 (t) and G(w”) = c* = const, where 5(t) is a delta function, The exact solution for k 
(t, T) is convenient for theoretical computations, and defines the limiting accuracy attainable in linear systems. 


The class of nonstationary random processes described by Equation (7) is seemingly sufficiently large to 
inciude a great number of the cases arising in applied problems, However, it is frequently the process's correla- 
tion function, R,(ty, tg), which is given, and it is then necessary for the finding of k(t, r) to determine the opera- 
tors My and L; which, in many cases, is not easily done, On the other hand, as is well known, the weight functions 
of optimal systems may contain delta-functions and their derivatives and are therefore impossible to use directly 


for the construction of optimal filters, In this case it is necessary to approximate the exact solution by asymptotic 
(or otherwise approximating) functions, 


The difficulties mentioned above may be by-passed if approximate methods are used to determine the 
optimal system's weight function, 


The suggested approximate solution is based on the methods of functional analysis, Application of these 
methods to the problem of determining the weight function of an optimal system is, so far as the author knows, 
made here for the first time. This application is related to the choice of the Hilbert space in which the solution 
exists, to the choice of the formula for succcessive approximations and with the proof of the convergence of the 
successive approximations in the metric of the selected space, 


The sole approximate method known to the author is in [2], where the solution is obtained in the form of a 
series which converges to the exact solution, The approximate solution can be obtained by limiting oneself to a 
finite number of the series’ terms, The method is quite convenient for an approximate determination of the weight 
functions of linear systems with variable parameters, However, the amount of computation connected with this 
approximate method of [2] is about double that required by the method presented in this paper, This amount of 
computation is determined, first of all, by the calculations connected with obtaining the coordinate functions of 
the canonical expansions and, secondly, by the calculations connected with obtaining the terms of the series for 


the weight function, Moreover, the method of [2] creates some additional difficulty in the case of systems with 
constant parameters, 


2, We now prove that the problems of extrapolation, interpolation, and filtering reduce to the determina- 
tion of an element which gives a minimum value to a certain quadratic functional, * 


If we write conditions a and b of section 1 in expanded form, then the functional 


where k(t, T), for any ¢ € 7’, is an element of Hilbert space Hy ** , must have a minimum with the additional 
conditions 


t 
t)de(t)dt (t) 2). (9) 
t—T 


The problem of the conditional minimum of the functional I,(k) leads to the problem of the unconditional 
minimum of the functional 


t t t 


t—T t—T k=1 t—T 


* A functional, I(k) is said to be quadratic if ( k + &z) is a quadratic trinomial in €, where k and z are arbitrary 
elements of the Hilbert space H, 


** A definition of the space Hy, will be given below, 


(5) 

(8) 

s, then 
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where the unknown Lagrange multipliers, y(t) are defined by the additional conditions (9), and t may be consid- 
ered as a parameter, 


Thus, the problem of an unbiased estimate of the mathematical expectation, (1), of the random function 
y(t) leads to the determination of the k(t, r) which makes for a minimum of the functional Kk), The best (in the 
sense of a minimum of the second-order moment of the random error) biased estimate also leads to a minimum of 
a quadratic functional, 


Let us, for example, consider the case when the coefficients a, in (1) are random variables and the mixed 
second-order moments, yj j = M{ aja;} (i,j =1,2,...,m), are known, The optimal system must carry out filter- 
ing, i.e., its output signal must be so close to the function g(t) that the quantity 


I,(k) = M{ k(t, g (uy) 


t—T 


has aminimal value, If we expand the expression for (k) and discard the terms which do not depend on k(t, r) 
we get 


t—T 
(12) 
k(t, *) du (x) ae. 
k, j=1 


In obtaining formula (12), it was assumed that the coefficients a, are independent of the random function 
n(t), With this, the correlation function Ry(t,, tz) has the form: 


Ry (ti, ta) = Pe(tr) (te) + Rn (t1, ta). (13) 
k, jm 


Thus, the problem of determining the weight function, k(t, r), for the optimal system in this case also leads 
to the search for a k(t, r) which gives a minimum value to the functional (k) which has the same form as in the 
case of the unbiased estimate [cf., (10) ]. 


It is not difficult to show that the problems of extrapolation (leading) and interpolation also lead to the de~- 
termination of a k(t, tT) which gives a minimum value to functionals analogous to (10) and (12). For this it is 
only necessary to use, as the “desired” signal, g(t + A), where A is the lead, In the problem of interpolation A< 0, 
while in the filtering problem A = 0, Consequently, the problems of leading (extrapolation), filtering, and inter- 
polation lead to the search for the minimum of some quadratic functionals, We therefore consider a method of suc- 
cessive approximations for finding the function k(t, r) which gives a minimum to some quadratic functional of 
general form, We rewrite (12) in the form 


I (k) = (Ak, k) — (14) 


Here, ¢ is an element of the functional Hilbert space H, The space H is obtained by defining, on the linear 
set M of functions (linear) [3], in a general way, a scalar product, (x, y), of functions x and y, and by then adding 
to the elements of the set M all limiting elements, x* , of sequences { x,} such that Lim m ll Xn —x*|| +0. With 


this, all terms of the sequences lie in the set M, but the limit elements may not lie ia - 


Thus, H is the set M plus its limit elements, and, consequently is a complete space, We note that, for the 
problems of leading and filtering, we may consider H to be an incomplete space, i.e., we may consider that EM, 
where M is a linear set of functions with a scalar product defined on M, In the sequel, it will suffice to take M 
as the set of continuous functions. 
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For the case of (10), the operator A obviously has the form: 


on t 
n the Ak = \ k(t, Ra(t, y)dt,, *—-T<t<t (15) 
1m of t—T 
xed and the function g(t, r) may be written in the following manner: 
lter- n 
P(t, t) = (t) de (t)- (16) 
(11) By virtue of the symmetry and positive definiteness of the correlation function R,(t;, ts), the operator A ts 
symmetric and positive, i.e., 
7), 
(Az, y) = (x, Ay) (17) 
and 
(Az, z) >0 (18) 
(12) 
forx,y 2, 
It is now necessary to construct a Hilbertspace H, which is known to contain an element k*, which gives 
a minimum to the functional in (14), For this we use the already considered linear set of continuous functions, 
tion M, and we define a scalar product and norm on it by the formulas 
(z, y)a = (Az, y), zy € M, (19) 
(13) lz, = (Az, 2). (20) 
It is easily verified that the scalar product and norm thus defined, by virtue of the symmetry and positiveness 
leads of A, satisfy the usual axioms for norms and scalar products. 
n the We now add to the set M all the limit elements of the sequences {xp} which converge to elements x* in 
the sense that: 
he de- lim 
0, 
inter- The set M, thus completed, forms the Hilbert space Ha. It is obvious that the space Ha contains elements 
d of suc- which lie neither in the linear manifold M nor in the space H, In particular, if M ts the set of continuous func- 
il of tions then H, obviously, is the space L, (i.e., the space of function vectors, the sums of the squares of the com- 
ponents of which converge); with this, Hy can contain elements lying neither in M nor in Lg, It is easily shown 
tht MCHC H,. 
14 
a8 Let the functional ( ¢, k) of (14) be bounded * for p€ A and &k€H,. We may then use the theorem of 
 Meeat Riesz, according to which there exists a unique element k” € H4,such that 
adding 
). With (p, ==(h, (21) 
On the basis of (14) and (21), we obtain 
for the 
at (k) == (k, k)a —2(k, (22) 
ke M 


* If the function (¢, k) is unbounded in Hg, the functional K(k) is obviously unbounded below, and the problem of 
its minimum is meaningless. Indeed, if (¢, k) is unbounded in Hg then, for any integer n, an element k, may be 
found such that || k,|| 4 =1 but | (¢, k)|> a. 
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It follows from this that the minimum is attained for 


k = k’, (24) 


= — (25) 


Thus, if the functional Kk) has a minimum, it is attained for k = k*, with k*@ H,. Certainly, it may hap- 
pen that &° € H, but this is not the general case. 


Consequently, there are three cases. In the most general case, the element minimizing Kk) lies in the space 
Ha; in the second particular case, A” € 17: finally, in the third case, the functional cannot in general have a mini- 
mum, if(¢, k) is unbounded in Hy. The most important is the general case, Therefore, further consideration is 
given to an approximate method for determining k°€ H». For an approximate determination of the element k * 
which minimizes the functional in (14), we use two direct methods for constructing minimizing sequences {kp}, 
for Kk), A sequence {k,} 1s called a minimizing sequence for a functional I(k) if 


Jim I (kp) = (26) 


We shall now use the two following sequences as minimizing ones: 


(Lk,,, ) 
hats = ka — n=0,1,..., (27) 
Lky = Alin — (28) 
and the sequence 
= In — aban, (29) 
where a is some positive number, 


The elements kp and x, in formulas (27) and (29) are chosen arbitrarily in H, In particular, for the case 
of the functional in (10), the quantities entering into formulas (27) and (29) have the following forms: 


Lk, (t, t) = \ ky (t, 1) t) —(t, *), (28a) 
t—T 


where the function g(t, r) is defined by formula (16), 


t 
(Lktny = Lky (ty 2) Litn(t, 2) de (30) 
t—T 
and 
t 
ALky (t, t) \ Lkn(t, Ra(t, dt. (31) 
t—T 
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(23) 


(24) 


(25) 
y hap- 


space 

tion is 

ntk* 


Kn}, 


(26) 


(27) 
(28) 


(29) 


(28a) 


(30) 


(31) 


In functional analysis, formula (27) is related to the method of steepest descent [4] and formula (29) is 
a somewhat altered formula of the gradient method, with the condition that the number a is choren in the re- 
quired manner. The problem now consists in showing that both sequences, (27) and (29), are minimizing, in the 


sense of (26), for the functional of (14). We prove this first for the sequenc in (27). By taking (14) and (27) into 
account, we easily obtain 


Lk,)* 


(kn) —1 (knits) = (32) 


It follows from (32) that the sequence { kn)} is monotonic and decreasing, But, on the other hand, in ac- 
cordance with (23) the functional Kk) is bounded below by the number|| k *||®, . Under these conditions, the se- 
quence { I(k,)} must converge to the unique limit 1,,;,(k*). Consequently, we have 


I (kn) (kn41)—>0 for oo. (33) 


It follows immediately from (23) and (33) that 
k, — k’, Ak, — | 0 (34) 


Thus, the convergence of the sequence { I(k,,)} necessarily implies the convergence of the elements kp 
in the metric of space Hy. 


It is essential to note that in the case when the limit element k*€ H, and k* €H, the sequence { k,} 
diverges in the metric of space H, In this case, k* is a generalized solution of the problem, At the same time 
we note that from (32), (33) and the boundedness of operator A, it follows that 


(Lk,,, Lk,)* 


(ALR —0 for (35) 

ALky) [a(Lkp, Lkn), (36) 
(Lk, Lk,)® 

(Lkn, Lkn) S Thy (37) 


Lin) —>0 for 00. 


Thus, the limit element, k*, is simultaneously the solution of the linear operator equation 


We now prove that the sequence in (29) is also minimizing for the functional of (14) with the proper choice 
of theparameter a, By using (14) and (29), we may easily obtain the following formula; * 


I (Ken) — (Kings) = Lien — fa. (40) 
By taking into account the fact that operator A is bounded and self-adjoint, we shall have 


— a? | >0, 
if we choose o < 2/|| . 


* In formula (40), the symbol || || denotes the norm in space H, and the symbol || || 4 denotes the norm in 
space Hy. The definitions of these norms were given previously, 


- 
| 
Whence, 
(38) 
ase 
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Under these conditions, the sequence { K(k,)} will also be monotonic decreasing and, by an argument anal- 
ogous to the foregoing, we arrive at the result that the sequence in (29) is also minimizing for the functional I(k), 
Thus, we can use formulas (27) and (29) for the approximate determination of k*(t, Tr). 


It is necessary to mention that formula (27) is everywhere more rapidly convergent than formula (29). How- 
ever, the computations connected with (27) are more involved than those needed for the use of formula (29). In 
a large number of cases of practical application, it turns out that the difference in speed of convergence ts small, 
and that it is therefore more efficacious to use formula (29), We also note that in practice {1(k,)} converges very 
rapidly, and not more than one or two approximations are necessary. 


For the determination of the dispersion, D,;, of the random function at the optimal system's output, one 
may use formula (23). By using the fact that Dp; = || k*|| 4 we obtain 


Da, = —T mak’): 


For the solution of Equation (39), the question of existence and uniqueness of the solution are essential, As 
follows from (15), in our problem (39) is a first-order integral equation, The actual necessary and sufficient con- 
dition for existence and uniqueness of the solution is contained implicitly in the possibility of using the Riesz 
theorem in writing Equation (21), specifically, in the complete (by definition) space Ha, the function (¢, k) must 
bounded for all k(: 14. If the system of characteristic elements (eigenvectors) of operator A lies completely in 
H,*, the condition for the boundedness of the functional (¢, k) in Hg can be given in a different form, On the 
basis of (21) and (39) we may write 


= (Ak, k) = = >} (42) 
j=1 j=1 


where bj = (k, kj) and a; =(¢, ky) with Aj and kj being, respectively, the eigenvalues and eigenvectors of the 
equation 


Ak = Ak (43) 
The equation 
(Ak, k) = 3} BA, (44) 
j=t 


is the condition that the eigenfunctions lie completely in space Ha, Consequently, on the basis of (42), if the 
system of eigenfunctions is completely in Ha, then the condition which is both necessary and sufficient for the 
existence and uniqueness of the solution comes to the convergence of the series 


(45) 
7 


This assertion is analogous to the requirement of the well-known theorem of Picard [5] for the existence and 
uniqueness of the solution of a first-order integral equation in the space H = Ly this requirement reduces to the 
convergence of the series 


(48) 


j 


> 


* For an incomplete system of characteristic elements of A, the conditions for the existence of a solution are to 
be found in [2], However, in this case, there may be an infinite set of solutions. 
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Since the positive A, tend to zero for j -» a, the requirement in (45) is less stringent than the one in (46), 
This should be expected, - 2 H CH ,. We note that the completeness of the system of characteristic func- 


tions (elements) in Hy is an essential condition, not only for the uniqueness, but also for the existence, of the so- 
lution, 


If the system {kj} is incomplete then, even for a continuous function ¢ M, a solution may not exist as, 
for example, in the case of the equation 


k(t) cos (t — t) dt = const. 


Therefore, the question of the existence of a solution must be investigated in each individual case, 


In conclusion, we shall consider the following simple example of the approximate determination of an 
optimal system's weight function, Let the input signal have the form: 


g(t) = ao + art, 
where a» and a; are unknown coefficients, 


The random function n(t) (noise) has the correlation function 
Ry (try te) = 


It is necessary to determine approximately the system's weight function, k* , which provides the best un- 
biased filtering. In this case, operator equation (39) has the form: 


k(t) — t)dt = Ao + Aut, 


owns 


where the numbers \» and A, are defined by the conditions of lack of bias. 
T T 
0 0 


The solution of this equation has the form [1]: 
k* (t) = Ag + + (t) + Did (tT), 


where 5 (t) and 5(t— T) are, respectively, delta~functions at the points t= 0 andt = T, aot Ag, Ay, ky, and Dy 
are cetain coefficients which it is possible to find. 


We now determine the approximate value, xp, for k*, For this, we use formula (29), setting x(t) = 1, We 
set a in formula (29) equal to one-half. It iseasily verified that, with the numerical values of the parameters, 6, 
y and T taken in the example, the conditiona < 2/|| A|| holds, 


On the basis of (29), a second approximation, x,(t), was found, For simplicity of the computations, the 
values of 8, y, and T in the formula for the correlation function were taken as 8 = 1, y = 20 and T=1, Then, 
the following expression was obtained for xq(t): 


q(t) = 4.029 — 6.088 t + 0.1322 — 0.0454 — 0.0125 + 0.0125 
0<t<i. 


Since both weight functions, the exact one k* and the approximate one x,(t), provide an unbiased estimate, 
in order tocompare them we computed the dispersion of the noise at the output of each of the two systems, It 
turned out that for the system with k*(t), 
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Dn = 0.316 
and for the system with the weight function x(t), 
Das = 0.375 


The example considered shows that the approximate solution found, x9(t), for the optimal system's weight 
function is close to the exact solution k* (t), The system with weight function x(t) is obviously simpler to im- 
plement than the one with k®* (t). 


SUMMARY 


The approximative method considered in this work is quite simple and, in practice, contains no essential 
limitations, The direct search for a function which minimizes a quadratic functional allows one to avoid the 
complex formal procedures of searching for an exact solution, 


In the general case, the problem's exact solution, k*, is a generalized* one and, consequently, this gener- 
alized solution is determined approximately, With this, any approximation of the solution lies in the function 
space H = L, (the exact solution lies in Ha) and is always simpler to realize. The closeness of the exact solution 
to the approximate one is defined in the metric of Hy. 


For obtaining an approximate solution, other approximating methods may also be used (Ritz, Galerkin, 
the method of least squares, etc,). The essential difference of our method from these is that the approximation 
does not take a form chosen a priori, but rather a form determined by the solution itself, This allows a more 
rapid convergence of the approximate solution to the exact one to be obtained, 


The algorithms defined by formula (27) and (29) may be used either for direct computation, or in electronic 
analog computers, 


The authors wishes to express his great appreciation to V, S, Pugachev for many helpful comments and for 
participating in discussions of the results of this work. 
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GENERALIZATION OF THE SHAPING FILTER METHOD TO 
INCLUDE NONSTATIONARY RANDOM PROCESSES 


A. M. Batkov 
(Moscow) 


The paper considers the basicproperties of the impulsive responses of linear 
systems with variable parameters of general form, and the correlation functions of 
the random processes at their outputs when acted upon by “white noise”. The prob- 
lem of determining the characteristics of the shaping filters for this class of non- 
stationary processes is solved, and a methodology is suggested for using them in 
simulation problems in analyzing the dynamic accuracy of systems, 


1, Introduction 


One of the basic tasks arising in the study of automatic control system operation with random stimuli is the 
problem of defining the system's dynamic accuracy, i.e., the characteristic error, which is a random function of 


time. 


Today, several methods of solving this problem are in wide use, They may be divided into analytical and 
experimental methods, 


Analytical methods consist of finding exact or approximate solutions of the differential equations describ- 
ing the processes in the system, and in determining analytic expressions for the dispersion, or the probability dis- 
tribution function, of the error. 


One such method is based on the canonical representation of the random function at the system's input, i.e., 
its representation in the form of asumofnoncorrelated elementary components [1], In particular, for the inves- 
tigation of the stationary random processes at the outputs of linear automatic control systems with constant pa- 
rameters, the methods of analyzing dynamic accuracy which are based on the relationship between the spectral 
densities of the input and output signals [2] are widely used, However, their use in practical problems is quite 
limited due to the lack of methods for solving a broad class of differential equations, and also to the complexity 
of determining an explicit form of the random function's canonical representation, 


Experimental methods are based on the statistical processing of the random processes at the output of the 
system being investigated. In the case when these processes are stationary, the problem of analysis may be solved 
by the use of correlators or spectral analyzers [2]. If the input processes are nonstationary , or if the system 
changes its parameters in time, it is necessary, in solving the problem to carry out an averaging over the ensemble 
of a number of realizations. The disadvantages of the experimental methods are the necessity for large quantities 
of data to obtain sufficiently accurate results, and also the complexity of the apparatus for processing these data. 


Recently, several other approaches to the solution of this problem have been noted for those cases when the 
automatic control system under investigation is linear (or a linearized model of a nonlinear system) [3]. The 
idea consists of this: that by starting from the characteristics of the stimuli and of the system being investigated, 
one obtains analytical relationships between the statistical characteristics of the input and output processes in a 
form which is convenient for solution on computers, by means of which the statistical characteristics of the output 
signal are also determined, 
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This method leads, from the practical point of view, to sufficiently simple operations in the case of linear 
systems described by linear differential equations with variable parameters, with stimuli in the form of "white 
noise”. If the input random signal is different from “white noise* then the problem of analysis based on this 
method reduces to the determination of a system by means of which a random process with the given correlation 
function is obtained when “white noise” acts on it, Such systems have received the name of shaping filters [3]. 


The problem of determining the shaping filter is only solved for the case when one is given the correlation 
function of a stationary random process with a rational-fractional spectral density, Below, we present a method 
for determining a shaping filter for a broad class of nonstationary random processes. This allows the develop- 
ment of a single methodology for analyzing the accuracy of linear systems with random stimuli by means of 
analog devices, 


3. Posing of the Problem 


We consider a linear dynamic system with variable parameters and an impulsive response , k (t, T), which 
is an unknown function of the two variables t and r. Let this system be described by a linear differential equation 
with variable coefficients, of the form 


Lip, p=. (a) 


We assume that the system is first found in its unstimulated state and that, at time t = ts (in particular, ty = 
= —oo) “white noise" was applied to its input, In Equation (1), x(t) and f(t) are, respectively, the output and input 
functions of the system, and L(p, t) and M(p, t) are unknown linear differential operators of the form 


Lip.t) =} a @) 
=0 


M (p.t)= bald) 
u==0 


where the coefficients aj(t) and by(t) have the necessary number of derivatives, and n and m are integers such 
thatn > m, 


By our assumption, the mathematical expectation of the f (t) equals zero, and its correlation function is 
5 -function® , i.e., 
M (t) = 0, 
Ry (t,t) = M [7 (t) (*)] = 8 —*). (4) 


The problem consists of thist that from the given analytical expression for the correlation function, R(t, T) 
of the nonstationary process x(t) one must determine the impulsive response, k(t, T), in the regiont > T, or 
the operators in (2) and (3). 


It is obvious that if the given correlation function of the process x(t) is Gaussian and is in fact obtained as 
a result of a linear operation on “white noise", then the random process at the output of the shaping filter will 
have the same probability distribution function as the given one. 


If R(t, 7) corresponds to a process which is obtained from a nonlinear transformation of “white noise” then 
the shaping filter will have an output process whose second moment equals, or is near in the case of approxima- 
tions to R(t, T), that of the given process, but with different higher moments. 


To solve the problem posed, we consider certain properties of the impulsive response of the system in (1) 
in the general case, when the operator M(p, t) is not constant, and also the properties of the correlation functions 


* Obviously, if the mathematical expectation of f(t) is different from zero and is known, then it may be compen- 
sated for at the input or output of the shaping filter, thanks to the linearity of the filters being considered, 
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of nonstationary processes at the outputs of linear systems of the form (1) with variable parameters and with 


aad "white noise” at their inputs, 
3, Properties of the Impulsive Responses of General Linear Systems with Variable 
Parameters 
| As is known [4], the impulsive response, G(t, rT) of a linear system of the form given in (1) with M(p, t) #1, 
my defined as the solution of the equation 
L(p. t)G(t,t) = 8(t (5) 
with zero initial conditions, has the following properties, 
1, It is a solution of the equation 
ich 
ation L(p,t)G(t, 7) =0 (6) 
with the conditions 
0, a=0,1,....2—2, 
(1) G (t, t) -| 
at* = a=n—i. (1) 
input 
2. If G*(t,r) is defined as the solution of the equation 
L’ (p, t) G* (t,t) = —8(t —*) (8) 
(2) 
with zero initial conditions® , i.e., 
L* (p, =0, (9) 
+) 
" then, on the basis of Lagrange's formula [4], 
G (t,t) = t). (11) 
(4) 
3. It follows from (11) that 
L* (p, t)G(x, t) = —8(t —2), (12) 


i.e, in the case that M(p, t) = 1, the impulsive response, as a function of the second argument, satisfied the ad- 
d as joint differential equation, 


‘i 4, It is clear from (6) and (12) that, fort > r, the function G(t, r) can be given in the form 

(1) where the ¢(t) and g*,(t) are fundamental systems of the homogeneous differential equations (6) and (9). 

— In the case of the general linear system with variable parameters, described by Equation (1) with the operator 

mpen 


* In formula (8) and in the sequel, the superscript star denotes the adjoint operator, 
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M(p, t) defined by expression (3), one may establish the following properties of the impulsive response, k(t, Tr), 
which is the solution of the differential equation 


L (p, t) k (t,t) = M (p,t)8(t—*) (14) 
with zero initial conditions, 
1, As is known [4], 


k (t,x) = M"(p,t)G(t,), (15) 
2. If we denote by k*(t, rT) the solution of the equation 
L* (p,t) (t,t) = — M*(p, t)8(¢ —*) (16) 
with zero initial conditions, it then follows from (15) that 
(t,t) = M (t, 2). (17) 
In addition, based on (11) we have that 
M (p, t) k(t, t) = M* (p, t) t). (18) 
In particular, if the operators M*(p, t) commute, then in analogy to (11), 
k (t,t) = k(x, t). (19) 


3. It can be shown [5] that the impulsive response, k(t,r), as a function of the second argument, T, satis- 
fies the differential equation 


Rp, x) k(t, t) = (p, 2) 8 (t —*), (20) 
where the operators R(p, t) and Q(p, t) have the form: 


RP.) = 


v=0 
and the coefficients tp (t) and q,,(t) are determined from the system of algebraic equations in the known coef- 


ficients of the operators in (2) and (3) (cf., [5]): 


— Ayers (2) (0) + (— 1)! gu (0) ae = 0 (23) 
(A=0,1,...,2 +m), 
al 
where C= 
4, The impulsive response corresponding to Equation (1) satisfies the homogeneous equation 


L(p, t)k (t,t) =0 (24) 
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(15) 


(16) 


(17) 


(18) 


(23) 


(24) 


with the conditions [5] 


(a =0,1,...,2—m— 2), 
| 3 (25) 
v=u 
k(t,*) 
r=n—m—} =t 


il 
where cj = 


5. Analogously, k(t, T), as a function of the second argument r is the solution of the homogeneous equa- 
tion 


R (p, 2) k(t, =0 


with initial conditions which are obtained from (25) and (26) by replacement, respectively, of coefficients a; 
by byt, and q,,. 


We note that the values of the derivatives,“ *‘":*)| can be determined also from the coefficients 


of the operators in (2) and (3) from the formula [5] 


(u>0) 
(a=n—m—i1, —m,...,n—4), 


where the derivative with respect to r of the impulsive response of the operator L(p, t) equal to G(t, Tr), for 
t= rT is determined from the relationship 


G (t, 7) i a'G (t, +) n—i—1 k+i—l) 


(i>0) 


(4 =0,1,...) 


6. It follows from (24) and (27) that the impulsive response has 2n-1 continuous derivatives* with respect 


to the variables t and r in the region where t # rT, It , together with n— m — 2 derivatives, is continuous for 
ter, 


The higher derivatives of k(t, r) have discontinuities of the first type on the diagonal t = r, which may be 
determined from formulas (25) and (26). 


Moreover, k(t, T) can always be presented in the form 


where y(t) and » ,(t) form fundamental systems of Equations (24) and (27). 


* It is assumed that the coefficients of the operators: in (2), (3), (21) and (22) can be differentiated n times. 
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(14) 
| | 
= 
(19) 4 
4 
(20) q 
q 
4 
i 
q 


4. Properties of the Correlation Functions of Nonstationary Processes at the Out- 
put of Linear Systems Acted on by "White Noise" 


Basing ourselves on the assumption of section 2 and the results of section 3, we now establish the properties 


of the correlation function, R(t, r), of the nonstationary random process at the output of a system described by 
Equation (1) with “white noise” at its input. 


Since “white noise” is applied to the system's input, 


k(t, \) k(t, A) dh, t<t, 


ree 


(29) 
R (t,t) = 
i, 
From which, by applying the operator L(p, t) to (29) and (30), we obtain* 
L(p,t) R(t,t) =0, t>+, (31) 
L(p,t) R(t,t) = M(p,t)k(z,t), (32) 
From the symmetry property of R(t, T) it also follows that 
Lip, t)R(t,t) =90, (33) 
L(p,t) R (t,t) = M (p,c)k(t,t), (34) 


From these relationships we can obtain the following properties of R(t, rT). 


1, Fort >r andt<1r, R(t, r) has, like the impulsive response, 2n-1 continuous derivatives with respect 
to t and to Tr. 


2, 2(n—-m-— 1) derivatives of R(t,r) are continuous for t = Tr. 


3, Starting with the (2n — 2m — 1) *th,the derivatives of R(t, r) with respect to both variables have dis- 
continuities for t = r, the magnitude of which can be determined, 


In particular (5), for k = 2n- 2m—1, 


k—n+m 
t=n—m—1 l=n—m—1 ar! at toot 


From Equations (321) and (35) we have 


R (t,t) = os (36) 


t<t, 
“1 


where ¢;(t) form a fundamental system of Equation (31), and the 7 ;(t) are particular solutions of Equation (32). 


* These equations, with the proper initial conditions added, can be used for the experimental determination of 
the correlation function of the class of processes being considered (Cf., [5]. 


erties 
d by 


(29) 


(30) 


(35) 


(36) 


on of 


5. If the operators L{p, t), M(p, t) and M*(p, t) commute, and k(t, r) is identically equal to zero for tb, 
T < t) then, on the basis of (31) and (32), it may be shown [5] that 


R (t,t) = M (p,t)M'(p,t)H (t,t), t,t >to, (37) 
where H(t, is the Green's function®, i.e., 


a) H(t, r) is continuous, along with its derivative up to the (2n — 2)'nd order in the entire t, r domain; 
b) the (2n — 1)'st derivative of H(t, r) has a jump att = r equal to (-12-/ aa T) 
c) fort # r, H(t, 7) satisfies the homogeneous self-adjoint differential equation 


L(p, t) L* (Pp, t) H (t, t) = 0. 


5. Determination of the Shaping Filter's Differential Equation 


We now consider the problem of determining the operators in (2) and (3) for the shaping filters being sought 
in the case when the correlation function is given, or is presented in the form of (36). 


Since, by assumption, the operators have the forms of (2) and (3), the problem reduces to the determination 
of the numbers n and m and the coefficients a;(t) (i = 0.1, ... , m) and by(t) (Hw = 0.1, ...., Mm). 


We first find the operator L(p, t). 


Based on (36), the correlation function R(t, 7) may, in the region t > 1, be presented in the form of a finite 
sum of products of the functions g(t) multiplied by the y,(r). Therefore, by putting R(t, r) in the form of (36) 
for t > rT, we can determine the fundamental system of homogeneous Equation (31), With this, the order n of 
operator L(p, t), is determined by the number of linearly independent functions of the variable t in this expansion, 


We know from the theory of differential equations [6] that the operator L(p, t) can be determined from a 
fundamental system of solutions by expanding the determinant: 


L(p,t)z(t) = (38) 


By cancelling out the common factor, we obtain the form of the operator L(p, t), and also the coefficients 
a, (t) (i = i. n). 


For the determination of operator M(p, t), we rewrite (32), on the basis of (15), in the form: 
L (p,t) (t,t) = M (p,t) 7). (39) 


As is known [4], the impulsiveresponseG(t, T) corresponding to operator L{p, t) can be determined from the 
functions g(t): 
Gilt) a(t)... Pn(t) 
i(t) 


(— 


CL) = 


where W (rT) is the Wronskian for Equation (6), 


* Sometimes H(t,1) is called a one-sided Green's function [4]. 
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|| 
(33) 
(34) | 
Pa-++Pn 
q 
(20). 


From Equation (39) we can find the result of the action of the product of the adjoint operators M(p, t) M* 
(p, t) on the known function G(r, t). From this we determine the product, and then decompose this into its ad- 
joint factors, thus obtaining M(p, t). 


This method turns out to be very useful if M(p, t) has a simple expression and, in particular, does not con- 
tain the differentiation operator, Otherwise, the problem consisting of the decomposition of a differential self- 
adjoint operator with variable coefficients into its adjoint (conjugate) parts, is difficult. 


Another way of determining M(p, t) consists of using expression (35), defining the jumps of the derivatives 
of R(t, r) fort=r. Indeed, by knowing n and a,(t) (i = 0.1, ... ,m), we can determine from this formula the or- 
der, m, of operator M(p, t), analogously to the way this was done in section 5. The coefficients b pit) (# = 0.1, 
-+» » M) can be successively found from relationship (35), rewritten in the form 


u—l i 

a” k(t, a, (t) [ 1 k(t, 
= Bute t,t) — X 

x at [ at" n—m-—i—1 | + (41) 
b,, 
+ G a! a, (t) n m n m n ) 


By determining, from (41), the relationships between the derivatives of the impulsive response with respect 
to tand r for t = r, and then expressing them in terms of the known coefficients aj(t) and the unknown by(t), 
we may, from formulas (26), (26a) and (26b), successively obtain m algebraic equations in the by(t)*. With 
this, each coefficient, b pl» is determined in terms of the coefficients for the higher derivatives, by + ;(t), 
by + x(t), ... » bm(t). In particular, it follows from these formulas that 


bm (t) = an (t) V(—1 Ron (t, - (42) 


Thus, the problem of determining the differential equation of the shaping filter from a given correlation 
function reduces to the operations of differentiation, computing determinants, and solving second-order algebraic 
equations, With this, on the basis of formulas (41), (26),(26a) and (26b), one can set up a table of the coefficients 
b,,(t) which permits one to quite rapidly determine these quantities in cases of higher-order equations. 


Example, Let the correlation function R(t, +) of the nonstationary process x(t) ( t > 0) be given in the form 


3t 
» 


3t 


R(t, t) = 


In this case, n = 1 and g(t) = 1/ *. Consequently, by expanding the determinant of(38) and cancelling out 
the common factor we obtain 


dz 
L(p, t)x(t)=t + 2z. 
* Formula (26) may be used for the determination of the coefficients in the right member of differential Equation 


(1) from the known impulsive response k(t, 7) in the problem considered in [7]. For this, (26) can be presented 
in the form 


ark(t =) n—r—m—1 
r=n—m o=0 


m 
> (— (t) (p = m, m—i,..., 1). 
ver 
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From formula (40) we have that 


Git, t>t. 
Hence relationship (39) takes the form: 
ot t 
M(p, t) th a, t<t 


from which it follows that 


M (p, t) M* (p, t)f (t) = 9f (#) 
and 


M (p, t) = 3. 


Consequently, the shaping filter sought is described by the differential equation 
dz 
t + 22 = 


and can be realized by an amplifier with gain of 3/ 2 and an integrating RC-filter whose resistance varies pro- 
portionally to t/ 2, 


6. Application of the Method of Shaping Filters to the Analysis of Linear Systems 


The method presented for determining the impulsive response* and differential equation of shaping filters 
allows one to reduce the problem of analyzing the dynamic accuracy of linear automatic control systems with a 
broad class of nonstationary stimuli to the task of investigating a variant of a system with “white noise” at its 
input**, This provides the possibility of developing a single methodology for analyzing linear dynamic systems 
with stationary and with nonstationary stimuli, based on the use of analog (simulating) devices, 


With this, in order to determine the dispersion of a random process at the output of a linear dynamic sys- 
tem described by a differential equation 


A(p,t)y(t) = B(p, t) x(t), 


where A(p, t) and B(p, t) are operators of the form given in (2) and (3) and x(t) is the input random function of 
the system with a known correlation function R(t,r), it is necessary to proceed as follows: 


1, From the given function R(t, rT), we determine on the basis of section 5 and the Appendix, the impulsive 
response or the differential Equation (1) of the shaping filter ***. 


2. Since the problem of determining the dispersion of the system's output signal, when the impulsive re- 
sponse is k(t, T) and “white noise” acts on the input, leads, when an analog device is used, to the computation 
of the integral 


t 


D(t) = (t, t)dr, (44) 


* Cf,, the Appendix. 

** The methodology of employing shaping filters in the stationary case is considered in detail in (3). 

*** The problem of determining the differential equation from the impulsive response of the shaping filter 
may be solved on the basis of the results of [7]. 
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it is necessary to determine the differential equation which satisfies the impulsive response k(t, Tr) as a function of 
of the second argument, rT, for given 1 [8]. The function k(t, r) is the impulsive response of the series-connected 
shaping filter and investigated system with impulsive responses, respectively, of kg (t, T) and k(t, T). 


As is well known, for the case of a series connection of a system, 


t 
k(t, = Akg (hy (45) 


Consequently, by determining from formula (23) the operators in (21) and (22) from the equation found for 
the shaping filter, and applying them to (45), we obtain, based on (20), 
R(p, t) k(t, t) = Q(p, t) ki (t, (46) 


with the condition that the impulsive response of the shaping filter, kg(t, 7), as a function of r, satisfies the dif- 
ferential equation 


R(p,2)k (t, =) = Q(p,2)8(t — +). 


Here, the impulsive response of the sytem being investigated, kj(t, r), enters as a function of the second 
argument, T, for a given t. Therefore, by applying (23) to Equation (43), we find the operator equation 


C (p, t) ky, (t,t) = D(p, t)8(t (47) 


which defines k;(t, 7) as a function of r. 


In Equations (47) and (46) we now introduce the new variable, r; = t— T(r, >0), for determining k(t, r) 
in a physically realizable domain [3]. 


If we now, in accordance with Equation (46), im- 


input at) a plement on an analog device the series connection of the 
Fa } ~ 2 7 circuitry representing Equation (47), having to do with 
wl the system being investigated, and the circuitry for the 


equation obtained from the equation which satisfies 

k¢(t, T) as a function of r, replacing r by t— Ty, where 

T, is the time"experienced" by the model, we obtain the 
impulsive response of the whole system, k(t, t— 7,). The 
transformation of the original system to the system for solu- 
tion on the analog device is shown in Figs. 1 and 2, 


Fig, 1. The initial system, 1) is the shaping 
filter, L(p, t)x(t) = M(p, t)f(t); 2) is the sys- 
tem being investigated, A(p, t)y(t) = B(p,t)x(t). 


a(tt-G) met Dit} 


Fig. 2. The analog (simulation) scheme, 1) is the initial 
conditions for r; = 0; 2) is the model of the system being 
investigated, C(p, r,)kj(t, t— T 4) = 0, 3) is the model of 
the shaping filter, k(p, y)K(t, t— Ty) = Q(p, 7); 
4) is a squaring device. 
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u- 


We note that, for the determination of k(t, rT), it is more advantageous to solve on an analog device the 
homogeneous differential equation of the type of (27), rather than Equation (47), L.e., 


C (Pp, ky (t ,t-— 7) = 0, (48) 
with the initial conditions for r,; = 0, defined by formulas (25, (26) or (26a). 


3. The final step in solving the problem consists in squaring k(t, t—7,) and computing the dispersion at 
a definite moment of time by formula (44): 


t—t, 
D(t) = \ (t,t — (49) 


By repeating the process of solution for different t, we find the curve of temporal variation of the disper- 
sion of the output random process. 


We note that in [8] there was presented another method of solving the analysis of linear systems, based on 
treating the correlation function of the input nonstationary random signal R¢(t, T) in the region t > r as the im- 
pulsive response of a linear system with variable parameters, 


However, the method presented in [8] turns out to be more complicated to realize than the method of solv- 
ing this problem which was suggested in this paper. Indeed, for simulating the system shown in Fig. 6 of [8], it 
is necessary to implement a scheme which transforms the impulsive response of the inverse system,* (t,t—T), 
as a function of r to the function of W*(t, rT), and then to apply this to the input of the direct system, Such a 
transformation may be realized by recording W*(t, t— T) and then tracing the reproduced recording in the op- 
posite direction, which is, as a rule rather complicated to do. 


With the use of shaping filters, as shown, it is only necessary to square the output of the inverse system and 
integrate the result, which may be implemented by means of standard blocks of analog devices, 


In conclusion, it is necessary to mention that the method suggested here for determining shaping filters also 
permits determining of the correlation function of the output process on an analog device, as was shown in [5]. 


APPENDIX 


Determination ef the Impulsive Response of a Shaping Filter 


We now consider the problem, posed in section 2, of determining the impulsive response, k(t, T), of shap~- 
ing filter from a given correlation function, R(t, T), We shall assume that R(t, T) has the properties considered 


in section 4, i.e., that it corresponds to a nonstationary random process obtained from “white noise” by subject- 
ing it to a linear transformation, 


Mathematically, this problem reduces to solving Equation (29) [or (30) ] with respect to the function k(t, Tr). 


The method of solution consists in reducing this equation to a Volterra integral equation of the second type, 
and in solving this by means of a resolving kernel* , 


To do this, we proceed as follows, 


We differentiate Equation (29) n — m times with respect to each of the variables t and r, where m and n 
are determined by the orders of the operators of the shaping filter sought, (2) and (3). 


On the basis of (25) we have, for ty = 0. 


dd + 


* The author turned his attention to this means of solving the problem after becoming acquainted with (9). 
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’ rent 


Or with the notation 
a" (t, +) 


(t,2) a" (4,2) 


y(t, t) = 


(51) 


(8, 2) 
0,A)= 5 
B (0, ») > (52) 


f 


we obtain the Volterra integral equation of the second type, with kernel 8(t, A), which is the unknown function 
(cf., 


t 
a(t. 2) = BU, A+B (53) 
0 
If we denote by a(t, u) the resolving kernel of 8(t, u), i.e., if we define a(t, u) by the equation*® [10] 


uw) = a(t, wu) + \B (A, u)a(t, A) dr, (54) 


eco * 


we may then show that a(t, u) is the solution of the Fredholm integral equation of the second type 
t 
a(t, u)=(t, u) — u)a(t, A)dA, t>u. (55) 
0 
In fact, since on the basis of (30), in the region of t>r n(t, T) is defined by the expression 


0 


then by dividing the interval of integration in (55) from 0. to u and from u to t, by substituting (56) and taking 
(54) into account, we arrive at an identity. 


After determining a(t, u) from (5), one can, using the properties of the resolving kernel, determine 8 (t, u) 
in the ordinary way. However, for solving the problem it is necessary to find the impulsive response k(t, rT). 


Therefore, returning to Equation (29) and differentiating itn— m times with respect to t, we obtain 
t 


with the notation 


* The lower limit here equals u, since 8( A, u) = 0 for A < u, which corresponds to the physical realizability of 
k(A, u). 
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a" —™ Rit, 


(t, tT) = . 58 


Amt 


Since the resolvent a(t, A) of the kernel A(t, A) is known, then we have, from (57), 


0 


Thus, the problem of determining the impulsive response of a shaping filter reduces to the solving of in- 
tegral equation (55) for a(t, u) and the computation of k( 1, t) from expression (59), 
It is obvious that, to solve (55) and compute (59), it is necessary to know the difference n— m and the mag- 


n—m—1 k x 
nitude of : ee er These quantities may be determined from the given correlation function 


R(t, T). 


Indeed, on the basis of section 4, the correlation function of a random process at the output of a system 
described by Equation (1) has 2n — 2m — 2 continuous derivatives with respect to t fort = r, Consequently, if the 
k'th derivative of R(t, r) with respect to t has a jump for t= Tr, then n—m =(k+ 1)/2. 


Further, the magnitude of the jump in the (2n — 2m — 1)'st derivative of R(t, 7) with respect tot, fort=T, 
equals, on the basis of (35), 


— 6-9) =] eo" ™—1 (60) 
Since [5] 

then 
£70, d) 
|, — om — 1(0,0)- (62) 


We can thus establish the following order for determining the impulsive response, k(t, T), of the shap- 
ing filter: 


1. By differentiating the given function R(t, rT) in regions t > r and t< rT with respect to t, and by deter~- 
mining the jumps of these derivatives for t = r, we find the difference, n— m = (k + 1)/ 2, where k is the order 
of the first derivative whose jump (saltus) for t = r is different from zero, 

2, By multiplying the magnitude of this saltus by (—1) "-™~* and extracting the square root § we deter- 


a 


. Then from formula (51), we determine n(t, rT) for t< rT and 


t=! 


termine the magnitude of 


t> 


* It is obvious that the indeterminacy at this point does not affect the final result, since both impulsive responses 
will correspond to one and the same function, R(t, rT). 


1061 


3. By solving either analytically or approximately the Fredholm equation of the second type in (55), we 
find the resolvent a(t, u). 


4, By computing y(t, T) from formula (58), we find the desired impulsive response from formula (59). 
Thus, based on an application of the method presented, the problem of determining k(t, T) from Equations 
(29) and (30) leads to a simpler problem, namely, that of solving a Fredholm integral equation of the second type, 
which may be done by one of the known methods [10]. 
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THE EQUATION AND CERTAIN PROPERTIES OF AN AUTOMATIC 
CONTROL SYSTEM'S ROOT LOCUS 


N. N. Mikhailov 
(Moscow) 


The paper considers the root locus equation for the general case of linear 
automatic control systems, as well as the root locus properties deriving from this 
equation, 


The root locus of a linear automatic control system is the set of trajectories of the roots of the system's 
characteristic equation when some parameter or other of the system varies, 


If we denote the system parameter to be varied by w then, in the general case, when the characteristic 
equation is solved for this parameter, it takes the form: 
P™ + P™ + ap’ +--+ +a,p+a, 


(1) 


where k, br-1, --- » Dos An-1, --- » &9 are Constant coefficients which do not depend on w, As wvaries, the roots 
of Equation (1) are translated in the complex plane of the roots (in p plane) along the arms of the root locus, If, 
with this, the parameter w is also presented in the form of a point on a complex plane (the w plane) then, since 

w is a real quantity, this point will be translated along the real axis of the wplane, In the general case, the para- 
meter w may vary within the limits of —o to +a, Consequently,the root locus may be defined as the mapping 
of the w plane's real axis into the p plane by means of function (1). 


The expression in the right member of (1), a rational-fractional function of p, we denote by F(p) for brevity: 


t be (p — 21) (p — 2a) ..-(P—Sm) 
q” + an—ip™ *+... + + a (pP—%)(P— 92) (2) 


Z(p) _ 
Here, Z;, Zz, ... , aNd qy, Gg, ... are, respectively, the zeroes and poles of F(p). 
Thus, 


w = F (p). (3) 


The root locus corresponding to the function F(p) with m zeroes and n poles, will be denoted by T(m, n). 
The parameter w, while varying, remains on the real axis of the wplane, This means that the function F(p) is 
real, i.e., the argument of F(p) is real, i.e., the argument of F(p) must equal to zero,sm, where s is an integer. 
Consequently, the equation of the root locus whixh corresponds to Equation (1) will, with varying w , have the form: 


arg [F (p)] = sx. (4) 
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Equation (4) is equivalent to the equation 


tg [arg F (p)] = 0. (5) 


By bearing in mind that the argument of a rational- fractional huitton, F(p) equals the sum of the angles, 


6zj, of the vectors to the point p from all the zeroes, zj, minus the sum of the angles, 6q,,, of the vectors from 
all the poles, q,,, of the function F(p) (Fig. 1): 


arg (p) = 2 arg(p— si) — (6) 
and by using the well-known trigonometric formula 
tg(e+y) = 


Yeh [1] transformed Equation (5) to the form 


m+n 
> >} (Il,tg = 0, (7) 
B=1 Y=1 


Me where 6, is the argument of a vector to the point p from 
a a zero, z; , or from a pole, qy. Arguments of the vectors 
Fig. 1 from the poles enter the formula with reversed signs, 


The application of this formula in obtaining the 
equation of a concrete system's root locus might be illustrated by the example, shown in Fig. 1, of a system whose 


function F(p) contains four poles (qy = 0, +jwy, Gg = jwy, and and one zero (z,). 
We write this function in the form 


pat p+ bo 
(Pp — 91) — 92) (P — 98) (P — 9) tap 


In this case, Equation (7) can be written as: 


(tg — tg — tg — tg 0,,—tg 0,4)—(tg8,, tg tg + te tg + 
+ tg 6, tg 0, tg + tg tgO,, + tg tg + tg 6, tg %,, — 
— tg tg tg 0,, — tg 0,, tg tg —tg tg —tg tg tg, )+ 
+ tg,8,, tg tg = 0. 
\ 
If we substitute in this the expression for the tangents of the angles, 


we then, after uncomplicated but tedious transformations, obtain the final equation of the root locus of the sys- 
tem considered 


w + 2 + ag) + + 2g) + + (aby — + 
+ + 2(—2by +5) + — 0. 
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In an analogous way, the equations of the root loci for other functions F(p) can be derived. The table pro- 
vides the final expression for various T(m, n) which were thus obtained. 


A comparison of these equations suggests that they may all be described by one general formula of the form: 


(2) Q (0)—Z (0) Q _ 
(a) Q*(e) Z() (a) ZY) (a) Q (a) ) 


(8) 


4. 272) Q* (0) __ 2° (0) Q" (9) (0) Q™ _ 
3! 2! il 4! 


Here, Z(o) and Q(a) are the polynomials in the numerator and denominator of the function F(p) in which, 
instead of the variable p, one substitutes 0, i.e., the abscissa of the variable p- 


On the basis of induction, we may assume that formula (8) does not apply only to the cases cited in the table, but, 
is a general equation of the root loci for all systems of the type considered. 


In succinct form, this formula is given as: 


¥ sin (—1)'C# ZY (a) (2) | =0. (9) 


If the function F(p) has no zeroes, i.e., F(p) = 1/ Q(p), then formula (9) is simplified, and takes the form: 


m+n 
ve wQ’(c) wQ"(c) wQY) (a) 
sin | Al —... =O, (10) 


G. A, Bendrikov [2] was-led to the same equation by an investigation of the root loci of systems without 
zeroes, This serves to strengthen the assertion made above that formula (9) is the general equation of root loci, 


Use of expression (9) instead of formula (7) significantly reduces the arduousness of the operations of set- 
ting up the locus equations. However, the fundamental value of formula (9) consists in the fact that it establishes 
a direct relationship between the current coordinates w and o of the root locus, and may therefore be used for 
investigating the general properties of such types of loci, Formula (7) does not give this capability. 


It follows from formula (8) that the locus equation is not changed if the polynomials in the numerator and 
denominator of F(p) are interchanged. This means that, for example, the root locus will be one and the same for 


the function F(p) = 1/Q(p) and for the function F(p) = Q(p). Therefore, the root locus equations in the table were 
only given for n> m, 


The character of formula (9) is such that the root loci will be the same for the function F(p) = ¢(p) and for 
the function F(p) = — ¢(p), i.e., the form of the root locus does not depend on the sign preceding the function F(p). 


The general Equation (9) decomposes into two equations: 


o=0, (11) 
sin (—1)'C, (0) = 0. (12) 
ver] i=0 


It follows from these expressions that the root locus always contains the real axis ( w = 0), on both sides of 
which are symmetrically distributed the locus arms defined by Equation (12). Some of the arms, to be specific 
| 2(n — m— 1)| of them go to infinity, Thus, the number of arms which tend to infinity is, including the real 
axis, 2(n ~ m). 
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| 


but, 


P). 


By using the ordinary methods of mathematical analysis, we can show that the locus arms, which tend to 
infinity have linear asymptotes which intersect in one common point, the center of asymptotes o¢... The center 
of asymptotes lies on theo axis, and its abscissa equals [3] 


= i= 


Thus, 0¢ a, would be the center of gravity of the poles and zeroes of F(p) if the poles were considered to 
have unit positive mass and the zeroes to have unit negative mass, 


The slope, ¥,, of the asymptote to the real axis equals 


fa = [i= 1, 2,...,(n—m— 4)]. (14) 


The asymptotes, together with the real axis, divide the angle 24 into 2(n — m) equal portions, 


By using general Equation (9), we can find the points of intersection of the root locus with the imaginary 
coordinate axis of the p plane, at which points the roots of the automatic control system move into the right half- 
plane, These points are found as the real roots of the equation 


[orsin2Z 53 1) = 0, (15) 


i=0 
obtained from (9) by setting o = 0 in it, 


We are ordinarily interested in the values of parameter w for which the roots of the characteristic equation 
lie on the imaginary axis. Computation of these values of w is carried out by formula (3) in which, instead of p, 
the real roots of Equation (15) are substituted, 


To find the points of intersection of the locus arms with the real coordinate axis, the branch points, it is 
necessary to set w = 0 in Equation (12). Then, the desired branch points are determined as the real roots of the 
equation thus obtained: 


Z' (a) Q(s) — Z(s) Q’ = 0. (16) 


It is obvious that the maximum number of branch points equals m+ n —1, 


It should be noted that a branch point is nothing other than a real zero of the function Fp), A complex 
zero of the function F'(p) is a branch point of the root locus only if the point whose image it is lies on the real 
axis of the wplane. If this condition does not hold, then the locus does not have complex branch points and, con- 
sequently, its arms intersect nowhere, in addition to not intersecting on the real axis of the P plane. 


The trajectories of the roots of Equation (3), corresponding to the changes in w, will be modified if some 
other parameter, for example, a zero or pole of F(p), changes its value, 


Figure 2 shows how the root locus of the equation 


p—z 
p(p -+ 4) + Gp + 13)’ 


is deformed if the zero, z, is translated along the real axis. As is obvious from this figure, the root locus may 
undergo a qualitative change as the parameter z passes through certain characteristic values of the given equation, 
Thus, if z = —4.61 (Fig. 2, c), two arms of the locus touch the real axis, thus forming a new double branch point. 
With a further shift of z to the right, this double branch point is bifurcated, thanks to which the number of arms of 
the locus is increased. ° 
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Such a characteristic point is also z = —4,.37, As is clear from Fig. 2,e, in this case two branch points 
merge into one double one, which disappears with a further shiftof z to the right, due to which the number of 
locus arms again changes, 


Fig. 2 


The point q, should also be considered a characteristic point. As is clear from Fig, 2,c, g and h, as the 
zero zZ passes through this point part of the locus arms disappear. The same may be also said, finally, of the point 
qi, with this sole difference, that as the zero z passes through it, new locus arms appear (Fig. 2, j). 


We call such characteristic values of the parameter being varied “critical points". It is obvious that a 
critical point is no other than a value of the parameter being varied for which a double branch point appears on 
the locus, 


We now find the conditions satisfied by critical points, Let the parameter to be varied be a zero, z;, of 
the function F(p), as in the case of Fig. 2. We factor out the term containing this zero from the polynomial Z(p), 


presenting the polynomial in the form Z(p) = (p — 2;)P(p). 
With this way of writing the polynomial Z(p), Equation (16) for the branch points takes the form: 


(o — 2) — — 1) P’ (2) + P (2) Q(2) = 0. a7) 


As z; varies, the real roots of Equation (17), the branch points of the original Equation (3), will be translated 
along the real axis. As tiie zero zj passes through a critical point, two real roots will merge into a double root, 
and will then separate into two complex conjugate roots, Thus, to find the critical points, it suffices to construct 
the trajectories of the roots of Equation (17) as a function of z;, to find the branch points of these trajectories, and 
to determine the values of z; corresponding to the branch points found, These values are also the critical points 
of the varied parameter z; of the original Equation (3). 


It is convenient to write Equation (17) in the form 
| P (a) Q’ (a) — P’ (2) — P (2) 8) 
P (8) Q’ (6) — P’ (2) Q (a) P 


The equation of the bfanch points of the root locus corresponding to Equation (18) is set up from formula (16). 
After providing such terms, the desired equation for the branch points takes the form: 


Q (0) {P (c) [P (a) (2) — (0) (2) — P* (2) + (2)}*) = 0. (19) 
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As follows from (19), the branch points are the real roots of the equations 


Q(c) =0 (20) 


and 


P (2) [P (2) (a) — (9) Q’ (2) — P* (2) Q + (2) [P’ = 0. (21) 


One easily convinces oneself that the left member of Equation (21) is the numerator of the function [Q(o)/ 
/P(o))". Consequently, the branch points of the root locus of Equation (18) are the real zeroes of the functions 
Q(a) and [Q(o)/P(o)]”. 


In order to find the critical points, it is necessary to substitute, into Equation (18), the values of the branch 
points that have been found. Substitution of the real roots, q,,, of Equation (20) gives 


(22) 


(ay) — P’ Q(9,) 


Since Q(qy ) = 0 [since q,, is a zero of polynomial Q(c)], Equa~ 
tion (22) reduces to the form 


Zi = (23) 


Thus, the critical points are all the real poles of the original 
function F(p). In addition to the real poles of F(p), the set of critical 
points contains the values of z; obtained from (18) when, instead of 
o , one substitutes the real zeroes of the function [Q(c) /P(c)]”. 


We now give an example of the use of the formulas derived for 
finding the critical points, whereby we use the equation 


P+ 4p + 5p 
Fig. 3 


in which the parameter to be varied is the zero, z. It is obvious that 
P(o) = 1 and Q(o) = —(0° + 407 + 50). In this case, Equation (18) is written in the form: 


0Q’(c)—Q(s) + 40? 
Ve) (4) 


and Equation (21) reduces to the form Q*(a) = 0 or 60 + 8 = 0, from whence o, = — 4/3, Substitution of this 
value of o, in (24) gives the first critical point of the varied parameter, zc, = — 64/9, The second critical point 
is the real zero, q; = 0, of polynomial Q(o) (the other two zeroes of this polynomial are complex), The root 
locus of the initial equation for z = z,, = — 64/9 is shown in Fig. 3, 


Ordinarily, in designing a system by the root locus 


dianeae " plane method, it is necessary to know the value of the parameter 
root w which corresponds to any point p; lying on the locus, 
ocus Computation of the magnitude of w sought is carried out 
= / by formula (3) with the substitution p = p,. 


The converse problem, i.e., citing thepoints on the 
locus which correspond to a given w, is more difficult since, 
for this, it is necessary to solve Equation (1). 


bj The giving of a real value of we is defined by the 
™ b point of intersection of the real axis of the w platie with a 
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line parallel to the imaginary axis of the w plane and separated from it by the distance wo (Fig. 4,b, line ab), In 
this case, the points of the root locus which correspond to given values of wy will be defined as the points of in- 
tersection of the locus with curves which are the images of the parallel line ab (points «', a", and 8 on Fig. 4,a), 
With a translation of the line ab parallel to the imaginary axis of the wplane, there is a corresponding translation 
of its image on the p plane which will intersect the locus in new points, Each position of the image curve on the 
Pp plane correspond to a definite value of wo. 


Certainly, constructing the images of line ab by means of the corresponding computations is a no less com- 
plicated task than the solution of Equation (1). However, if this computational task is mechanized, the method 
considered may be successfully [4] used for spotting values of w along the locus. 
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ON A TRACKING DRIVE WITH GYROSCOPES ON AN 
OSCILLATING PLATFORM 


S.G. Kolesnichenko 


(Moscow) 


A device (a direction finder) is investigated which is mounted on an oscillating 
platform with a power gyroscope drive, automatically tracking heavenly bodies and 
measuring the absolute angular velocity of its own line of sight, The false constant 
component of angular velocity is discovered. 


With the launching of the Sputniks there arose the problem of automatically tracking them and other heavenly 
bodies, Devices for this already exist. One might cite as an example the sun-following system (coelostat) set up 
in the Crimean Observatory, A description of a similar system is given in [1]. The existing coelostats are not 
capable of operating onoscillating foundations (platforms) such as, for example, ships or airplanes, This is almost 
always due to the poor quality of the electric power drive which turns the coelostat's mirror. 


With oscillating platforms, electric motors may not “have time" (due to lags) to change the relative posi- 
tion of the coelostat's mirror and, in this case, the mirror loses the heavenly body from its field of view. 


In the present paper we suggest the use, as the power drive for slewing the coelostat mirror of a two-axis 
gyrostabilizer of the classical type. The design and equations of motion of such gyrostabilizers are given frequently 
in the literature, for example, in [2, 3]. It is known that these power drives, with comparatively small dimen- 
sions, have high speeds of response in comparison with electric motors, This permits us to look forward to the 
possibility of using coelostats even on oscillating platforms, in particular, on board ships, 


1. Kinematic Relationships of the Coordinates 


We consider three orthogonal right-handed coordinate systems (Fig. 1). The first system {OX, OY, OZ} 
is attached to the platform on which the gyrostabilizer is mounted, We denote the vector of the platform's ab- 
solute angular velocity by 2 ={Qx, Ay, 2} where Qx, Qy, Mz are the vector's components along the X, Y, and 
Z axes respectively. 


The second and third coordinate systems are related respectively, to the exterior frame of the Cardan joint 
{ OXp, OYR,"OZR} and the gyro-nodes {OXg, OY¢, OZG}. 


We denote as follows the vectors of the absolute angular velocities of the frame and the gyro-nodes relative 
to their own coordinate systems 


= OR:} and @ = {Wz, Wy, 


If we introduce into consideration the relative dihedral angle A between the { OX, OY} and { OXp, OYR} 
planes, then the kinematic relationship of the first and second coordinate systems is determined by the relation- 
ships between the components of the absolute angular velocities: 
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= 2, + A, 
Ory = Q, sin A + Q, cos A, 
= 2, cos A — Q,sin A. 


ail 
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| 
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Fig. 1. Gyrostabilizer design. 


If we consider the relative dihedral angle B between the planes of the frame {OXp, OYp}and the gyro-nodes 
{ 5X, O¥G} , then the kinematic relationship of the second and third coordinate systems is determined by the 
relationships: 


= — Sin B + wr, cos B, 


Wy = Op, + B, (2) 
@, = Wp, cos B + sin B. 


The abosolute angular velocities of the gyroscope with respect to the axes of precession, 0,5; and O45, are 
denoted, respectively, by: 


= — Wy, B abs 


2. Equations of Motion of the Gyroscope 


For setting up the gyroscope's equations of motion relative to the precession axes, OS; and O,S,, it suffices 
to use Newton's second law and the well-known rule of Foucault which takes into account the action of the gyro- 
scopic moments. 


If we denote by p and o the so-called controlling moments of the external force, of which we shall have 
more to say later, we may present the equations sought in the following form: 


—idap, — m(a) — Ho, — Hoe = 0,” 
( 


o — —n(b) — Ho, — = 0, 


id 

‘gol a> 
= 
A wy 
aN = 
4G 

7 
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where the kinetic moments of the gyroscope are denoted by H, the moments of dry friction at the bearing of axes 


0,5, and OS, are denoted by ~m(&) and —n(8), and the moments of inertia relative to these axes are denoted 
by {and j. 


3. Equations of Motion of the Gyro-Nodes 


We first investigate the kinematics of the discharging device with respect to measuring axis OY; of angle 
B. For this we turn to Figs. 1, 2, and 3, 


Let the discharging drive for the measuring 
axis of B provide a torque p directed along the mount- 
ing axis (Fig, 1,8), coincide with axis OY, Then, for 
A #0 and B # 0, this torque is transmitted, not only 
to gyro-node axis OY,, but also to axis OX,, which 
is an impediment to gyrostabilizer operation, To de- 
termine the components of p along axes OYg and 
OX, . we consider the spherical triangle A(II, 2, 3) 
on Fig. 2, Two of its sides are known: 3, 2 = 9/2, 
2, IL = A, and the angle between them is < (3, 2, II)= 
=a —B, 


We now determine the spherical angle < (2, 3, 
Il) and the side 3, Il, which we denote, respectively, 
x and x, 


By using the well-known rules for solving spheri- 
cal triangles, we find that; 


Fig. 2. Diagram of the angles. 
tg = sin Btg A, cosx = —cos Bsin A. (5) 


We can assume that all the constructions on Fig. 2 are made on a sphere of unit radius, 


In this case, we consider the external force applied to the yoke (7 on Fig. 1) via the mounting (8 on Fig. 1), 
lying in the plane tangent to the sphere at its point 3 (Figs. 2 and 3), At this point, the mounting 8 acts on the 
so-called yoke 7 (the cylindrical protuberance of the gyro-nodes in Figs, 1 and 3) with a force equal top sin x. 
We denote this force by the vector 3,4 (Figs. 2 and 3), Its projections on the mutually perpendicular gyro-node 
planes {OXg, OZ} and {OYg, OZG}, the vectors 3,5 and 3,6 respectively, define the torques acting relative 
to the measuring axes OY, and OX¢ (Fig. 2). 


In fact, it is easily seen that in the plane of triangle A(3, 4, 5), lying in the tangent plane, angle < (4, 3, 5) = 
=x and angle< (4, 5, 3) = 1/2. Consequently, leg 3,5 = 3,4 cos < (4,3,5) = usinxcos x. Since vector 3,5 lies 
in gyro-node plane {OX¢, OZg}, and 08 = 1, the discharging torque applied to the gyro-node and acting rela- 
tive to axis OYg may be found from the formula 


ty = x cos x. (6) 


In an analogous manner we find the component px of torque # applied to the gyto-node (via the yoke) and 
acting relative to axis OXG 


= psinesin y. (7) 


We denote the moments of inertia with respect to axes OX, O¥¢, OZ¢ by, respectively, Jy, Jy and Jy. 
The frame's moment of inertia with respect to measuring axis OX is denoted by Jp,. 


We now write the expression for the sum of the forces acting on axis OXG, these being the force of inertia, 
the gyroscopic force and the external force applied to the gyro-nodes via the mounting and yoke (Fig. 3) 


My = px + Haigh, — J xg. (8) 


? 
'G 
| 
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The similar total moment acting with respect to axis OZ¢ is expressed by the formula 
M, = Hag — J, — (9) 


In expression (9), the first and third terms in the right member are the gyroscopic moments acting on the 
gyro-nodes, 


If we denote by ~ M(A) the dry frictional torque applied to the gyro-nodes at the bearing of the measuring 
axis for angle A, then the equation of motion of the frame with respect to this same axis takes the form 


h+¢cosB—J,. wp, —M (A) + M,cos B+ M,sin B =0. (10) 


where o is the reaction of the torque transducer and ) is the discharging torque imposed by the drive for the 
measuring axis of angle A, 


The equation of motion of the gyro-node with respect to the measuring axis of angle B is found in the ana- 
logous manner 


— Jy — N (B) + HB abs= 0, (11) 
where N(b) is the dry frictional moment at the bearings of the measuring axis for angle B. 


4. The Controlling Device 


For imparting to the gyro-nodes the necessary orientation in space, or a forced rotational movement (fol- 
lowing), a controlling device (for example, a thermal direction finder [1]) is used. It operates as follows: 
Signals, proportional to some quantity D, are applied to the 
controlling device which is mounted on the gyro~-nodes, whose 
moments (p and o) act relative to the gyroscope’s precession 
axes (0,5, and O,5,), and are applied to the gyroscope casing 
(with reactions on the gyro-nodes). This quantity is a measure 
of the deviation of the gyro-node axis OZg, or of the optical axis 
of the ooelostat's mirror from a given direction, OZ (for ex- 
ample, from the line-of-sight axis to a heavenly body). 


For the sequel, it is convenient to introduce the fol- 
lowing notation (Fig. 4), 


Fig. 3, Kinematics of torque transmission. We shall call the half-line OZ, in the direction of 

which the gyro-node axis OZ, tracks, the sighting ray; we 
shall call the angle {OZ7, OZ} the defect, D > 0, and the dihedral angle {OXgZg, OZgZ7} the phase, ¢, 
of the sighting tay. 


In our notation, the dependence of the controlling torques, p and o, on the coordinates { D, ¢} of the sight- 
ing ray may be given in the following approximate form [1}: 


To -- p = HtDsin 9, To + =: HED cose, (12) 
where T is the time constant and € is a quantity called the quality of the controlling device and equal to the 
angular velocity of the gyro-nodes for D = 1, 


In the general case, the quality ¢ depends on the magnitude of the defect, D, Below we consider the case 
of this dependency when 


Bes & = const > 9 for D<D, (13) 
0 for D>D, 
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where D > 0 is the so-called visual angle of the controlling device, 


G The gyro-node axis OZg cannot follow the sighting ray OZ; when 
fad latter is outside of the visual angle of the controlling device, i.e., when 
D =D, 
y | Let w and & be the absolute values of the angular. velocity and ac- 
0 | celeration of the sighting ray relative to immovable space, If we ignore 
\ ” friction, gyroscope inertia (i, j) and gyro-node inertia, t.e., if we consider 
yO NI G an “ideal” gyroscope, then the necessary condition from tracking, as derived 
%G N from (12) takes the form 


Fig. 4. Direction finder co- 


To +o (14) 


This last inequality defines the class of sighting line motions which do not take it beyond the visual angle 
of the controlling device under the conditions of an ideal gyroscope. 


5. Range of Regulation Quality 


In the (— n) quality plane we now determined the region of asymptotic stability (after A, M, Lyapunov) 
under the most elementary assumptions, We first assume that angles A and B are quite small in absolute value, 
that the sighting ray OZ, is immovable, i.e., the angle of defect D = — A, and the phase g 0, In addition, we 
consider for clarity, ideal controlling and correcting devices, limiting ourselves to the simplest representations 
of them: 


p = HEA, = He, (15) 


where — and n are the constant quality coefficients of the controlling and correcting devices, 
Then, for an immovable platform, the contracted Equations (4) of the gyroscope motion take the form 


9 —ia —-ma— HA=0, (16) 


and for the external frame of the gyro-node Cardan joint, the form 
—MA4- Ha =0, (17) 


where J is the moment of inertia of the frame and gyro-nodes with respect to axis OX, and M is the dry frictional 
torque on the frame axis, 


For the system of (15)-(17), the characteristic equation is 


agp* + + + ap + a, = 0, (18) 
a, = Ji, a;=Jm Mi, a, = H* + Mm, 


a, = H%, a, = H%, n=E+ 


where 


All roots of this equation lie in the left half of the complex plane only if the Hurwitz determinant 
A i ( — 


This last inequality defines the region of stability, S, in the &, parameter plane, Its boundaries are de- 
fined by the €, n axes and the hyperbola 
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+ (1 + 2a)éq + an? — be — bn = 0, 
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For 0< € < b/a, this hyperbola may be approximated by the straight line segment 
(20) 


since = Bie Mm, which occurs in practice. 


With the tuning corresponding to the points of the segment in (20) there appear oscillations with the iden- 

tical angular frequency 
o= ——. (21) 


Indeed, at these points 


a Hy ‘m M 


6. Forced Gyro-node Oscillations 


We now consider the phenomenon of forced gyro-node oscillations which appear when the platform oscil- 
lates with respect to the horizontal axes, 


We assume beforehand that angles A and B are small in absolute magnitude and that the controlling device 
is switched off, i.e., p= o= 0. 


For sufficiently-small oscillation frequencies we may ignore the inertial torques, since Equations (4), (10) 
and (11) take the form: 


m (a) + Hw, + Hox = 0, (A) + Mz =0, (22) 
+ Hoy + Hof = 0, py —N(B) + HBaps=0, 
where My = Habs, Hy = My = B. 
If we limit ourselves tothe first powers of small terms, we may write Equations (22) in the following form: 
m (a) + Hw, =0, \X—M(A)+ Haads = 0, 


n (6) + Ho, = 0, (B) + abs= 0. (28) 
If the ‘platform oscillations are described by the law 
= {Q,, 0, 0} (24) 
it follows from relationships (1) that 
A= —Q,. (25) 


Equations (23) can be given in the form: 
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Thus, Equations (23) decompose into two independent systems, related to different, but similar, channels, 
It therefore suffices to investigate one channel for condition (24), 


From conditions (25), (26) and (3), and bearing in mind that n(B) « A, we find that 
M (Q,) +2 
From whence, taking the second of conditions (15) into account* , we obtain 
° M (Q,) 
a+ a= 
This equation, together with the equation 


m (a) 
completely describes the gyro-motion, 


It is clear from these last equations that the amplitude of the forced gyro-node oscillations depends only 
on the dry friction in the gyroscope's precession axes, and not on the friction in the exterior frame axes. 


From the functional relationship of the angular velocity and the friction w, = —m(d)/H, we find the ap- 
proximate value of the amplitude of the frame's forced oscillation 


where mg, is the maximum frictional torque on the gyroscope's precession axes. 


7. Gyro-node Discharge for Platform Oscillations with Respect to the Vertical Axis 


We now consider the gyro-node discharge for platform oscillations with respect to vertical axis OZ¢(Q = 
= {0,0, 2,}). 


The gyro-node angular velocity can be found from the formula 
w, = Q, cos B,, 


where Bp is the initial angle of deviation of the gyro-node. 
Angular velocities A and Q, have the same algebraic sign, while Mg(A) has the opposite sign. 
The gyroscope's angular velocity with respect to the precession axes equals 


@,, cosB>Ofor 2,>0, 
—, cosB for Q,<0, 


where ws, = Mg,/ H and Mg, is the maximum frictional torque on the frame axis, 


* For simplicity, the lag is considered to be zero, 


m (a) 
(26) 
M (A)—> N 
) 
) 
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The algebraic sign of the steady-state gyroscope rotation angle with respect to the precession axes, a = 
=- “tat cos is opposite to the sign of Mz. From whence the gyroscopic torque acting on the gyroscope's 
precession axes, M = Hwzsina, is not greater than zero, 


If the smallness of a is taken into account, and we substitute the values of a and wz, we obtain the expres- 
sion for the angular discharge velocity of the external frame 


(Q 
Asbs = — 008" B, <0. (27) 


If Q22max = 2% f Zp, where z» equals the amplitude and f the frequency of platform oscillation, then the average 
value of irame discharge velocity is found from the formula 


4w 
Aabs av = — —" cos* B, (28) 


SUMMARY 


1, With oscillations of the platform with respect to the horizontal axes, the amplitude of the forced gyro- 
node oscillations depends only on thedry friction in the gyroscope's precession axes, and not on the friction in the 
measuring axes (the Cardan joint axes), 


2. If the platform oscillations with respect to the vertical axis have amplitude z». and frequency f , and the 
slope of the gyro-nodes with respect to the external frame is the angle B = By > 0 ( A = 0 for t = 0), there is a gyro- 
node discharge by axis A with an average absolute angular velocity given by (28). 


The author expresses his deep gratitude to A, M, Letov and B, V. Shirokorad for their scientific direction 
under which the present work was carried out. 
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NOISE STABILITY OF CONTINUOUS PROGRAM- CONTROLLED 
SYSTEMS WITH MAGNETIC RECORDING 


V.N. Shadrin 


(Moscow) 


The potential noise stability is considered for multi-channel program-control- 
led systems with magnetic recording when phase modulation is used, The problem 
is solved for systems with weak fluctuating noise, resulting from disturbances in the 
magnetic recording channels, The expression for noise intensity at an ideal receiver's 
output is found for the cases of spatial (for multi-track recording), frequency and tem- 
poral separation of the channels. A comparison is made of the noise stabilities of systems with 
separated channels for identical effective voltage and unit width of recording carrier, 


1. Introduction 


A magnetic recording channel plus a program-controlled reproducing device may be considered as some sys- 
tem for transmitting information, In this case (Fig. 1), the transmitter is the recording block RB, the communi- 
cation channel is the recording carrier RC, and the receiver is the block for reproduction, BR, The noise source 
NS may include a source of fluctuating noise in the recording carrier itself and a source of external noise applied 
to receiver BR. 


The principal cause of the recording carrier's inherent noise is its magnetic inhomogeneity, A nonuniform 
distribution of ferromagnetic material leads to the appearance of random fluxes on the carrier's surface, In pas- 
sing through the core of the reproducing head, the random fluxes induce in its winding an emf known as the in- 
herent noise of the recording carrier, 


In addition to the fluctuations induced by the inhomogeneity of the carrier's magnetic layer, noise is ob- 
served in magnetic recording from mechanical vibrations of the carrier as it moves, from nonconstancy of the 
contact between the carrier and the head, from a delayed field in the reproduction head, from random oscillation 
in the carrier's speed of motion, and others [3-5]. 


The fluctuating noises listed are not related to the signals recorded on the carrier and, by virtue of this, 
possess the property of additivity. Under this same heading come external noises of the fluctuating type. In the 
sequel, we shall only consider noise of this type. 


In the analysis of multi-channel magnetic recording noise stability ,it is assumed that the transient noise 
between channels (cross-talk) is small compared to the disturbances mentioned above, and will not be taken into 
account in the computations. 


When magnetic recording is used in continuous program-controlled systems, one has to do both with trans-~ 
mission of continuous functions of time (direct modulation systems) and with transmission of the instantaneous 
values of functions (pulse systems of modulation). 


The theory of potential noise stability of V. A. Kotel'nikov [1] allows one to determine the very best pos~ 
sible noise stability with any method for transmitting information, 


In the case of direct modulation systems, i.e., those in which the oscillation (information) to be transmit- 
ted, F(t), enters as some function directly into the signal's expression, the noise stability, in accordance with [1] 
can be determined from the expression 


* 
c= (1) 


where o is the specific value of noise at one cycle, Af(t) is the signal to be transmitted and F(t) is the transmit- 
ted oscillation (information ), 


In formula (1) 


where T is the time interval in which the signals considered lie. 
A continuous program-controlled system with phase modulation 
(PM) may be subsumed under the heading of a direct system, since the 
RB ea BR oscillation to the transmitted (the phase, which is a function of the 
: , translation of the piece being processed) enters directly into the ex- 
pression for the signal, 


In the present work multi-channel magnetic recording is anal- 
yzed for identical effective voltages, taken off from unit widths of 
recording carriers, In magnetic recordings, the dynamic range is a no 

Fig. 1 less important factor. However, in the final analysis, for program- 
controlled systems, all relationships reduce to the power, an indicator, 
which may in some sense, be provided by the effective value of voltage taken off the recording carrier. 


The phase deviation, Ag, which enters into the expression for the signal to be transmitted, depends in con- 
tinuous program-controlled systems, both on the receiving transducer (in the given case, on the rotatable trans- 
former) and on the degree of connection of the transducer with the translations of the piece being processed, There- 
fore, the question of the optimal value of the phase deviation, Ag, for a given system pass band, Aw, since it in- 
volves other characteristics of the program-controlled system, will not be analyzed in the present work. 


NS 


All the methods of separating the channels are compared for the identical signal phase deviation in each 
channel, According to [1], the effective voltages are calculated for the cases F(t) = 0 and F(t) = cosQt by the 
formulas 


—T/2 
U2= = + A*(0, t) dt, 
wth (3) 


—Tp 


2. Spatial Separation of the Channels 


With spatial separation of the channels, each signal is carried on a separate track of the magnetic carrier, 
i.e,, multi-track recording is used. With this, the signals in each channel are written and reproduced by indivi- 
dual magnetic heads, 


The frequency band taken up by the individual channels in spatial separation equals 


Aw = 2A9Q, (4) 
where Ag is the phase deviation for PM andQ is the maximum modulating frequency. 
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With spatial separation, the magnetic carrier is used inefficiently due to the necessity of having inter-track 
gaps equal, for narrow tracks, to the width of each track, The signal voltage in each channel will equal 


U 
(5) 


where Us is the voltage taken off the entire width of the carrier and N is the number of the channels (magnetic 
tracks), 


If (5) is taken into account, one can present the signal in each channel at the reproducing device's out- 
put in the form 


Ap (t) = c08 [wet + (6) 


By using (2) and (6) we find that 


dA, 


If we compute, by formulas (3), the effective voltage , Ue, taken off unit width of the carrier, we obtain 


(8) 


On the basis of formulas (1), (7) and (8) we now find the expression for the noise intensity at the output of 
the BR: 


 2No’ 
AeUe (9) 


With multi-channel magnetic recording, the fact should be taken into account thatthe carrier's noise volt- 
age is inversely proportional to the square root of the width of a magnetic track [5]. 

If o is the specific value of the noise at one cycle with unit recording carrier width, then the specific noise 
voltage o° in one track will equal 


= 10 
Yan 
By taking (9) and (10) into account we find that 
> 
(11) 


Formula (11) shows that, with spatial separation of the channels and phase modulation, the noise stability 
of magnetic recording deteriorates as YN with increasing number of channels. The noise stability is improved with 
increasing phase deviation, A¢, and increasing effective voltage taken off unit carrier width, Ue. 


3. Frequency Separation of the Channels 


With frequency separation of the channels, each signal is carried in a definite region of the entire pass 
band of the system. With this, the signals are written and reproduced by common magnetic heads. 


For identical channels we have (with PM) 


Aw = 2AgQ (12) 
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U =". (13) 
In a system with secondary modulation (PM-AM), the frequency band equals 


Aw = 2(w, + ApQ). (14) 


In a system with PM-AM SSB (single side band) 


Aw = + AgQ. (15) 


For the analysis of noise stability in these cases, one can use known methods [2]which give, for a system 
with PM frequency separation 


N 
Spm == (16) 
for a system with PM-AM frequency separation 
V3No 
1 
PM AM ~ 
and, for a system with PM-AM SSB frequency separation 
No 


‘As is obvious from the formulas given, the noise stability is identical for frequency separation with PM and 
with PM-AM SSB. 


For given A ¢ and Ug, the noise stability of a PM-AM system using two side bands is worse, by a factor of 
¥3, thanin the cases of PM and PM-AM SSB, 


4. Temporal Separation 


In this case, the signals in each channel follow each other successivley in time, The total time for signal 
transmission is 


(19) 


where h is the signal granulation coefficient, a is the coefficient of time usage per channel (in practice, a= 
= 0.5) and Ty is the length of the individual signal transmissions. 


With temporal separation, the frequency band used by the system is 
Aw = 2 +=. (20) 


By using the methods given in [2], we get 


OT (21) 
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», Comparison of the Spatial, Frequency, and Temporal Separations of the Channels 


In considering the expressions obtained for noise stability with various methods of channel separation, we 
can draw the conclusion that, for a given phase deviation, A ¢, the best noise stability must be reached by the 


spatial and temporal means of separating the channels, The worst noise stab ility is provided by PM-AM frequency 
separation, 


As an example, let us consider the noise stability of a multi-channel magnetic recording system with the 
following data: the number of channels, N = 4; the system's pass band is w = 2m - 10° radians/ second, the car- 
rier frequency is = 2m *500 radians/ second; the maximum modulating frequency is Q = 2m +25 radians/ second, 
the effective voltage taken from unit carrier widths is U, = 4 millivolts; the noise voltage is 0.04 millivolts in 


the band w = 2m - 10‘ radians/ second;and the phase deviation is Ag = 1 radian, On the basis of the formulas 
provided, we find that 


Fig. 2 The example shows that, for N = 4 and for identical phase 


deviation, Ag, the spatial and temporal methods of channel sep- 
aration. have the best noise stability, and PM-AM frequency modulation has the worst, 


However, practice shows that it is impossible to write on more than two channels on the magnetic carrier 
(for h = 4, the signal frequency is 500 cycles and the speed of carrier motion is 190.5 millimeters/ second), This 
occurs because of the essential interaction between the channels due to the finite width of the magnetic impress 
on the carrier, 


With frequency separation, due to the small probability of all signals coinciding in amplitude, it becomes 
possible, with the same number of channels, to improve noise stability by a factor of 1.2 to 1.5 by increasing 
the signal voltage in each channel, 


The use of PM-AM SSB and PM frequency separation is not recommended due to the large phase distortions 
for PM-AM SSB and the practical inconvenience of PM, 


Thanks to what has been presented here, we may recommend the practical usage of two systems; spatial se- 
paration, and frequency separation using PM-AM with the two side bands, 


Figure 2 shows the dependence of noise stability for spatial and temporal channel separation on the number 
of channels, The noise stability of frequency separation was computed with account taken of the fact that it can 
be improved in practice by a factor of 1.5 by increasing the signal voltage in each channel. 


SUMMARY 


Comparison of spatial and PM-AM frequency separation of the channels allows one to conclude that, for 
small numbers of channels (N 23), the noise stability of these systems may, for the same phase deviation, be 
considered equivalent. For a larger number of channels (N > 3), spatial separation of the channels may provide 
better noise stability. 


Frequency separation, wither with PM or PM-AM SSB, and temporal separation of the channels are not 
recommended for practical use, although they give theoretically better noise stability than PM-AM, 
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AUTOMATIC INSPECTION OF THE POSITION OF THE 


BOUNDARY SEPARATING TWO MEDIA BY MEANS OF THE 
PHASE METHOD USING ULTRA-SHORT WAVES 


V. B. Brodskti 


(Moscow) 


A method for the automatic inspection of liquid level, based on the reflection 
of radio waves, is extended to the case of inspection of the position of the boundary 
separating two media, It is shown that the method considered allows this inspection 
to be implemented without essential measurement error when the dielectric constant 
of one of the media changes within wide limits, 


In works [1, 2] I considered an ultra-short -wave phase method for inspecting liquid level, This method 
may be easily extended to the case of inspection of the position of the boundary separating two immiscible media, 
with the condition that the upper medium has small losses and an invariant dielectric constant, and the lower is 
a good conductor, Frequently, however, the dielectric constant of the upper medium varies within broad limits 
during the technological process, its exact value remaining unknown, Examples of this might be media which 
are mixtures of two liquids with low dielectric losses and with variable complements of their components. In 
such cases, direct application of the phase method presented in work [1] leads to large measurement errors, 


t, 


Fig. 1 


Therefore, for use of the phase method to inspect the position of the separation 
boundary between the transducer and the coupling line there is introduced an intermedi- 
ate transformer, matching the coupling line with the transducer for some average value 
of the dielectric constant of medium I (Fig. 1), The rest of the measuring scheme's 
block diagram corresponds to the scheme described in work [1]. We shall assume in the 
sequel that one medium displaces the other, so that the transducer is always filled. The 
quantities which relate to the upper medium with the low losses, and to the lower me- 
dium, the good conductor, will be denoted by the subscripts, respectively, II and I. 
With a deviation of the dielectric constant of medium II from its average value, there 
is a change in the phase velocity in medium II and, thanks to this, a displacement of 
the minimum in the measuring line, However, with this, the intermediate transformer 
becomes reflecting which, in its turn, leads to a displacement of the minimum in the 
measuring line, These displacements are of different sign, Under definite conditions, 
the displacement of the minimum in the measuring line due to variations in the dielec- 
tric constant of medium II can be brought to a sufficiently small value, and then the 
instrument's indicator will be insensitive to variations in the dielectric constant, We 
now turn to the analysis of the intermediate transformer's operation, 


We introduce the following notation: Zo, is the characteristic inpedance of the coupling cable, Zo; is the 
characteristic impedence of the empty tranducer,€ is the dielectric constant of medium II, 8, is the image phase 
constant of the coupling cable, 6, is the image phase constant of the transducer when filled with medium 1, k, 
¢; and ¢, are parameters of the intermediate transformer, 1, is the distance from the left terminals of the 
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intermediate transformer to the minimum current in the coupling cable, J, is the distance from the right terminals 
of the intermediate transformer to the position of the short circuit, and) is the generator wave length. 


We consider a quarter-wave transformer, which is a symmetric T-network (Fig. 2). The parameters of 
such an intermediate transformer have the following values: 


Za Z 


where Zp is the parallel impedance of the four-pole. For the series impedance, Zs, of the four-pole, there is the 
relationsip 


(2) 


The quantities B, 1 , and 6, /, are related by the equation of the S-curve [1]. A Taylor series expansion, 
up to third-order terms, of 6,2, as a function of 8,7, in the region of values of 691, close to Nw (where N is an 
integer) leads to the expression 


Bilp = Mn — kAV— 4 — (33 


where 


By setting N = 0 and taking (1) into account, we get 


i 
Bil) = Mr — 7, Pals — (i — (a0). (4) 


Thus, the linear portion of 6,2, does not, in general, depend on the dielectric constant of medium II, in 
order that the dielectric constant of medium II have a small effect on the measurements, it is necessary that the 
non-linear terms be small in comparison with the linear ones, To estimate the effect of the dielectric constant 
on the measurements of the level, we introduce the average value of the dielectric constant of medium Il: 


i+e 
(5) 


where € ma, is the maximum vaiue of the dielectric constant of 
medium Il, With this, we shall assume that the dielectric constant 
of medium II can vary from unity to €;ax during the measuring 
process, Further, we will sc choose Z, for the quarter-wave trans- 
former, that when the transducer is filled with the medium with 
dielectric constant equal to €9, the transformer will be matching. 


ve Then, for the definition of Z, we get the expression 
Zea Z 
|Z, P= (6) 


The value of the non-linear coefficient n(Cf., [1]) for the limiting value of € will also equal the maximum 
relative error of level measurement due to the variation of the dielectric constant, 


For n, we obtain from (4) that 
-_ 2 
k (1 — (= VeL), (7) 
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where L is the range of translation of the boundary separating the media, By using Expression (1) for k, on the 
basis of (7), substituting the maximum value of dielectric constant and using (6)* , we obtain 


h = (8) 


In comparing this expression with the analogous expression which can be obtained from formula (6) of work 
[1] for k » 1, we find that the wave length must be ¥2(€ max ‘) times larger for the same value of n than the 
wave length in the case of level measurement, We now compare the accuracy of measuring the position of the 
boundary separating the two media with the accuracy of level measurement by the phase method, The measure- 
ment error with the phase method is proportional to the ratio of the phase wave length to the range of translation 
of the standing wave's minimum, With this, it is indifferent whether this quantity is determined in the measur- 
ing or the coupling line. We denote this quantity by € u, and find it for the coupling line, The range of trans- 
lation of the standing wave's minimum in the coupling line is Vé, times larger than the translation of the bound- 
ary separating the media, Thus, by using expression (8), we obtain 


®max—1 4 
For the phase method of level measurement we have the expression [1] 


Va’ (10) 


where n denotes the greatest relative deviation from linearity. Thus, the error of the phase method of measur- 
ing the position of the boundary separating two media, with the condition that the relative error due to variation 
of the dielectric constant in the process of measurement not exceed 7, is increased not more than two times as 
compared to the error of the phase method of measuring levels. Analogous results may be obtained if, as medium 
I, instead of a conducting liquid one uses a liquid whose dielectric constant is large in comparison with the di- 
electric constant of medium II, With this, a matching absorber must be placed at the end of the transducer, 


CONCLUSIONS 


The phase method of automatic inspection of the position of the boundary between two media, one of 
which has small losses and the other of which is a conductor, may be implemented also when the dielectric con- 
stant of the medium with small losses varies within wide limits during the course of the measuring process, The 
error, due to the variations of the dielectric constant of the medium with small losses, may be made insignificant 
by means of connecting, between the coupling line and the transducer, a quarter-wave transformer which matches 
the coupling line to the transducer for the average value of the dielectric constant of the medium with small 
losses. 
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* If € ,,ax is 80 close to unity that) < 2L, expression (8) is useless, and it is necessary to choose L from consider- 
ation of single-valuedness (L < A/ 2) and to use expression (4) for estimating the error, 
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ASYNCHRONOUS MOTOR SPEED REGULATION BY A 
VALVE-CONTACT CASCADE CIRCUIT* 


I. B. Semenov and V, Z. Yarina 


(Moscow) 


This paper considers a valve-contact cascade circuit for controlling the speed 
of an asychronous motor with a phase rotor, The basic relationships are provided, and 
a comparison with a valve cascade circuit of ion-controlled valves is given. It is 
shown that, by its principle of operation, a semi-conducting diode in conjunction with 
a series-connected synchronously-operating contact is a controlled valve. 


The development of power germanium and silicon diodes with small losses, and also the creation of power- 
ful contact devices operating with high accuracies, permitted the design of power controlled valves. Such valves 
consist of series-connected semi-conducting diodes and synchronously-operating contacts, The latter play the 
role of the control grids [1], Such controlled valves are very profitably distinguished from the controlled ion valves 
(mercury-vapor tubes, thyratrons, etc.) used today, particularly for operation at low voltages, and may be used 
for all the standard circuits for rectifying and inverting current. Thus, if ion valves have internal voltage drops 
of 15 to 25 volts, the voltage drop in a germanium diode is about 0.3 to 0.4 volts and from one to 1.5 volts for 
a silicon diode, The basic mode of operation of such contact devices is the switching mode, 


Transfer of the contacts occurs during the non-conducting portion of diode operation, which makes the con- 
tact system very reliable, 


1, Contact Operation with Diodes in Series with It. 


We shall demonstrate the principles of operation of the controlled valve by the example of a single-phase 
two-half-period rectification circuit (Fig. 1). In such a circuit, control of the rectified voltage may be imple- 
mented by displacing, in time, the moment of contact closure, Such a displacement may be obtained by a rota~ 
tion of the synchronous motor's stator or by varying the excitation current if use is made of a special synchronous 
motor with longitudinal-transverse excitation, i.e., by mechanical or by electrical methods, It should however, 
be mentioned that, in contradistinction to circuits with controlled ion valves, where the commutation angle y 
and the overlap angle ) are equal, in the circuit mentioned above, as a rule, A > y [2]. 


In the-circuit to be considered, regulation of the rectified voltage may, in theory, be implemented either 
by variation of the contact time, t- or with tc = const. 


The contact deyice is more easily built with tc = const, In this case, A and, consequently, also t,, must 
be so chosen that the conditions for sparkless commutation always hold, i.e., 


(1) 


* This work is part of an investigation into electric machine control, being carried out under the direction of 
M. V. Meerov, 


jes 


ilves 


1) 


where Y max is the maximum commutation angle possible in the circuit and a is the contact's switching angle, 


The magnitude of y 4x can be determined if the parameters of the concrete circuit and the maximum 
allowable current, I, max, are known. If condition (1) holds, the contacts will transfer during the non-conduc- 
ting portion of diode operation, {.e,, for values of current which are practically zero, 


Uy Ug 
pe? 


Fig. 1, Block schematic and current and voltage curves, D, and Dy are 
semi-conducting diodes, A and a are contacts, SM is the synchronous motor 
and Ch is a smoothing choke. 


For the circuit of Fig, 1, the commutation angle y is defined by the expression 


V21,2, 
E 


cos [a + 7] = cosa — ' 


where I, is the load current, Em is the emf in the transformer's half-winding and x-y is the reflected inductance 
of the tFansformer's secondary half-winding, 


The contact time tc which is thus chosen must be increased 


~ cycte when the “play” of the contacts is taken into account, The dis- 
: advantage of regulation with tc = const is a certain amount of 
limitation on the range of regulation, 


A controlled valve, consisting of a series connection of 
synchronously-operating contacts and a diode, may be used in 
circuits for the rectification, inversion and,regulation of ac volt- 
ages. It may also be used for regulating the speed of electric 
machines, namely: a) for speed control of dc machines (the 
controlled valve operates in the rectification mode); b) for speed 
control of short-circuited asynchronous motors (the controlled 
valve operates in the mode wherein it controls the magnitude 
of an ac voltage); c) for speed control of asynchronous motors 
with phase rotors by a valve-contact cascade circuit for speeds 
less than the synchronous one (the controlled valve operates in 
the inverting mode), 


Fig. 2. Block schematic, AM is the 
asynchronous motor, Ch is the smooth- 
ing choke, Tr is the transformer and A, 
a, B, b, C and c are contact devices, 
2. Valve-Contact Cascades 


In this paper we give the results of an investigation into 
speed control of asynchronous motors with phase rotors by a valve-contact cascade circuit for speeds lower than 
the synchronous one (in the motor mode of operation), 


Figure 2 gives the schematic of the valve-contact cascade, In contradistinction to the circuit suggested 
and developed by E. L, Ettinger and F, I, Butaev [3], the ionic controlled valve is here replaced by a union of 
contacts and semi-conductors in the inverter portion and uncontrolled semi-conductors in the rotor circuit. 
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For the circuit of Fig, 2, built with ion valves, the general expression for the rotational speed of the motor 
cascade may be written in the form [3] 


3 
2,34 U, cos + 4AU + 2p 2, + Ry+ Ry) 


(2) 
— 


nm =m} 1— 


where n, is the rotor speed of the controlled motor, n, is the synchronous speed, Uy is the phase voltage of the 
transformer's secondary winding, U, is the rotors phase voltage for ny = 0, 8 is the angle of advance of the con- 
tacts, AU is the voltage drop across the valve, ly is the constant current, x7 is inductive impedence of the trans- 
former reflected into the secondary wind ing, x), is the motor's inductive impedance reflected into the rotor, Rr 

is the transformers active impedance reflected into the secondary winding, Rjy is the motor's active impedance 
reflected into the rotor, For the free-running state, with I, = 0, we obtain from (2) the formula 


U. 
=m (1— 0088 — (3) 


If we now go over to the circuit with semi-conductors and contacts, the quantity AU, due to its smallness, 
is omitted from expressions (2) and (3), This somewhat increases the rigidity of the cascade's mechanical charac- 
teristics, 


We now determine what conditions must be satisfied by the ratio of the rotor and transformer inverter 
voltages. 


The free running speed of the valve-contact cascade is 


Nr =n, cos). 


For 8 = 90° (here and in the sequel we refer to electrical degrees), n, = my, i,e,, the upper limit of regula- 
tions does not depend on For = 0, ny = UT /U;). 


In the case when Uy < U,, the lower limit of speed is raised by the amount n(Uy /,). In the case when 
UT > U,, the full range of speed regulation is attained but, with this, it is necessary for starting the motor to 
give a certain lead angle, B start, defined by the relationship 


U, 
1 = cos os cos start= 


If Uy = U,, the complete range of regulation is implemented, and most efficient use is made of the trans- 
former, since B start = 0. However, as will be shown below, motor starting for 8 = 0 is impossible under the com- 
mutation conditions, 


In the actual system ,commutation will be reliable, starting with some 


minimum value of lead angle, 8 Where 
Uy 
cos Brin = (4) 
1 | wt 
ae deme We note that, for inverter systems with controlled ion valves, the mag- 
4! | nitude of 6 4, Which provides reliable commutation is about 30° [4]. 
The power of the transformer is determined by the sliding energy trans- 
iad sia ' mitted to the network, i.e., by the depth of regulation, For the case of regula- 


tion from zero, it approximates the power of the motor being regulated. In 
Fig. 3 the case when the rotor voltage is approximately equal to the network voltage, 
the inverter transformer can be omitted, 
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‘Or In the valve-contact cascade circuit under consideration, reliability of operation is determined by the re- 
liability of the contact apparatus. Therefore, the question of implementing reliable sparkless commutation is the ) 
decisive question as far as circuit reliability is concerned, 


2) In [5] there was formulated the general condition for sparkless commutation,and ways were presented for 
causing this condition to be met for circuits with commutating chokes, This condition may be found from a con- 
sideration of Fig, 3, and is written in the form 


T+, (5) 

ns- 

. where x is the play of the contacts, 

” The duration, v , of the choke's closing action is determined by 
its magnetic polarity reversal time [2]. The duration of the closing 
action of a semi-conductor is determined by the duration of reversed 

(3) voltage action on the semi-conductor, Therefore, when a semi-con- 
ductor is connected in the contact circuit instead of a choke, condition 
(5) remains in force, but the quantity v is now defined by the duration 

wo of reversed voltage action, 


In the valve-contact cascade being considered, the contact ap- 
paratus operates on an inverter circuit with a controlled grid. We now 
make concrete the condition for sparkless commutation for such an in- 


; : verter with the condition for controlling angle 6. 
' 1 The duration of reversed voltage action on the semi-conductor 
(Fig. 4) is expressed by the relationship 
' 
= 
(6) 
la- where 
4max™ 
m—Anax By substituting (6) in (5), we obtain the condition for reliable 
commutation for the given inverter 
Fig. 4 


(7) 


The magnitude of x is determined by the design features of the contact apparatus, and is a constant when 


wih its velocity is constant. In modern contact apparatus x = 0.2 to 0.4° and so this quantity may be neglected in the 
sequel, 
In the valve-contact cascade circuit, the maximum value of angle 8 = #/2, and therefore, according to 
— the condition for commutation, the upper limit 6,,,, = ™/2 <1, is always satisfied. The lower limit is bounded 
by the overlap magnitude, But it follows from (7) that, in all cases, the lowest value of overlap must not exceed 
the maximum possible commutation angle, Y max, 1.€., for the uncontrolled overlap : 
(4) 
max max= 4 min (8) 
tog We now determine y pay im the valve-contact cascade circuit. The magnitude of the commutation angle 
y can be determined from the relationship [3] 
ans~- 
gula- [ 2) 
= B — arc cos| cos : 9 
tage, 
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As is easily shown, the magnitude of y increases with increasing Ip and decreases with increasing 6, But, 
as is clear from Fig. 5, as the angle 8 decreases below some value 8, the magnitude of the current I, also de- 
creases, It can be shown that the maximum value of the commutation angle, y max, Corresponds to the intersec- 
tion by the controlled motor's external characteristic of the axis of abscissas at the point Ip max, 1.¢., is dete:- 
mined by the values of ly max and 8x, specifically 


ymax = x — are cos| cos Br + 
T 


a From Equation (2), for np = 0 we get | 
=» i U 32M 
4 cos = 7 — + 2R; + (11) 
~ 
a Fan ae af As was stated above, the most advantageous relationship 
Pde as bon! ig ;™! between rotor and transformer voltage is U, = Uy. If we assume 
—s / 4 that the transformer and motor parameters, as well as the maxi- 
ee x. ge mum possible current in the circuit, Ig max, at known, we then 
pinay ~ - AN get the expression for determining 8,: 


ly 


Igma 


Fig. 5. By > By > 8, >84 >8min- cos Be = 1 — £4 2R;+2Ry)- (12) 


Fig. 6. Contact apparatus, 200 volts, 200 amperes, dc, 
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. By substituting 8, and Ig max in Equation (10), we get the 


} value of y may and, consequently, in view of (8), Amin and 
fA, 8B min. With this there is lost part of the lower range of speed 


regulation, defined by the magnitude of (Fig. 5). 


If the conditions of use require that regulation be introduced 
} tev ] h \ from ny = 0 for Ig = O, it is then necessary to change the ratio of 
h transformer and rotor voltages so that the equality U;/ Uy = 
= cos B min is obtained, 
Definitive corrections of the magnitudes of ymax and 
Xmin May be found by the method of successive approximations. 
“rev 3. Experimental Investigation 
my The experimental setup included the following elements, 
4 1, An asynchronous type MKMAR-14 motor of 5.3 kilo- 
watt power; stator voltage of 220/ 380 volts; stator current of 
» 22.7/13.1 amperes; rotary speed of 920 rpm; rotor voltage and 


current, respectively, 70 volts and 49 amperes; efficiency = 0,79. 
2. Type VG-50 germanium diodes, 


Fig. 7. a) Current in the contacts, 3. A 200 ampere, 200 volt contact apparatus, prepared at 
voltage on the contacts and reverse the IAT AN SSSR (Fig. 6); x = 5°. 

voltage on the diodes fe = 40°, 

Ic = 30 4, Atransformer formed by a circuit of three single- 
same but with = 90°, = 12°, phase transformers of 5 kilovoltamperes, for 380/ 127 volts. 

= 60 amperes; c) current and voltage The magnitude of the overlap A was determined by the 
curves at the transformer output for method presented above. For the parameters of the experimen- 

8 = 90° and ly = 60 amperes. tal circuit, x = 0.0442 ohms, Ry = 0.0582 ohms, x), = 0.24 ohms, 


Ry = 0.0841 ohms, the nominal constant current is 1,,, = 60am- 
peres, ly max = 3l,n = 180 amperes. From (12), taking into account the A/A connection of the transformer, we 
cos 6, = 0,143 and that By = 81°, From (10) we obtain y pax = 2°. With contact play being figured in, Amin = 
= Ymax + 2X = 12° and 84, = 12°. 


acteristics 


ae 


0 as {0 5 


Fig. 8. External characteristics of the valve-contact cascade, 


For obtaining the full range of regulation, it is necessary to choose the voltage relationship ¥3U,/Uy = 
= cosB j= 0.978. In our case, ¥3U,/ Uy = 0.955 and Bstart = 17°. Therefore, By = 62.5° and ymax = 7 
i.e., A = 12°, 


With such a setup, overlap of the contact apparatus occurred without sparking over the entire range of 
regulation. As increased, the lower portion of the regulation range was lost, 


Figure 7a, b show the oscillograms of current in the contacts, voltage on the diodes and voltage on the 
contacts during the commutation process and Fig, 7c shows the oscillograms of voltage and current yielded by 
the network, Figure 8 gives the external characteristics of the valve-contact cascade. 
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SUMMARY 


1, A semi-conducting diode in series with synchronously-operating contacts is, by its principle of operation, 
a controlled valve which possesses small losses, reliable operation and small dimensions, Such a valve may be 
used in rectifier and inverter circuits connected in low-frequency schemes, It may also be used for controlling 
de and ac electrical machines, 


2. A valve-contact cascade, as compared with a valve cascade, has more rigid external characteristics, 
better energy indicators and permits the use of series asynchronous motors (in the case of motors of small rotor 
voltages). 


3. Reliable commutation is achieved by the proper adjustment of the contact's overlap angle which, in 
operation, remains constant, The magnitude of this angle is computed from the circuit parameters. 
LITERATURE CITED 
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AN APPROXIMATE METHOD FOR SOLVING EQUATIONS WHOSE 
ORDERS EXCEED THE CAPABILITIES OF ANALOG COMPUTERS 


G. M. Ostrovskii 
(Moscow) 


Investigators working with analog computers frequently find themselves in the position where the order of 
a system of differential equations exceeds the maximum value admitted by the actual computer, 


This note presents a method for the approximate solution of such problems using the technique of succes- 
sive approximations, 


Let the differential equations of order m, supposed to be investigated on the analog computer, have the 
form: 


dz, 


At the same time, let the computer be capable of solving differential equations up to order m,< m, We 
divide system (1) into two parts, such that the order of the first part will be n, and the order of the second part 
will be m—n, We assume that n < m, and m—n< mj, i.e., that both parts of the system of equations can be 
simulated on the given computer, Although it is indifferent into which groups system (1) is decomposed, it is, 
however, from the point of view of decreasing the difficulties, advantageous to divide the system into those 
groups which are related to each other by a minimum number of variables. We shall assume , for simplicity, that 
one variable from the first group, x, enters into the second group, and that one variable from the second group, 
Xn + 4, enters into the first group. We note that the method is also applicable to the general case of variable 
distribution, 


The system of equations now takes the form: 


= bn (2 


m 
“a = Im . +09 
We shall assume that the system is stable. 


Solution of systems (2a), (2b) will be carried out in the following manner, We set up the first group of 
equations, (2a), on the computer, In this group, variable xq 4 , from the group in (2b) enters, If some function 
of time xn + 1 = X°n + 1(t), is given for it, then there will be n equations in n unknowns in the first group. By 
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solving this group of equations, we find the functions x;r, xs), occ 9 xnCXt). The superscipt (1) indicates 
that this is the first approximation, If the exact form of the function x, , ; were known, then the solution found 
would be an exact solution of systems (2a), (2b), But since the function xp 44 was chosen arbitrarily, the solution 
found is not a solution of the general systems (2a), (2b). 


After performing the eretions on the analog computer, we set up the second group of equations where, as 
X,(t), we use the function x; Xe), which was found in solving the first group. We thus obtain a system of n—m 
equations in n~ m unknowns, By solving these equations, we find the functions x, (1) a(t), Xn 
After this, we again ue the first group of equations, (2a), on the computer where, as the function xp , ;, we 
use the function xn 1(t) which was found in solving the second group of equations, (2b). 


Solution of System (2a) gives a second approximation: x;(*)(t), x,t), . . x(t). The function 
x,(")(t) just found is used for solving the second group ofequations, etc. Thus, the k’th approximation will be 
related to the (k-1)’st approximation by the following equations: 


see 


dz. 


i.e,, after the (k — 1)"st approximation, the variable x, , en) is retained and used as the given function of time 
in carrying out the k'th approximation. 


The convergence of the given process may be rigorously proved in a number of cases (for example, for sys- 
tems of linear equations). 
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TELEMETERING SYSTEMS DEVELOPED AT THE ELECTRICAL 
ENGINEERING RESEARCH INSTITUTE OF THE MINISTRY 
OF HEAVY INDUSTRY OF THE RUMANIAN PEOPLE'S REPUBLIC 


R, Stere, M, Duma and M, Syrbu 
(Bucharest) 


Created in Bucharest in 1950, in accordance with the plan for electrification of the country, the Electrical 
Engineering Research Institute of the Ministry of Heavy Industry of the Rumanian People's Republic carries out 
investigations in the domains of industrial electronics, means of automation and industrial inspection; work is 
also carried on in the domains of heavy currents and electro-technical materials, in particular, semi-conductors, 


The Electrical Engineering Research Institute (ETNII) carried out,and carries out, a number of investigations 
in remote control, We mention several of them (in chronological order). 


A tele-control device for street lighting (A. Kretsesku, 1954). 


A tele-voltmeter and tele-ammeter for close operation in a rectification system (A. Kretsesku, M, Kalushitsa, 
M, Duma, 1956), 


A tele-signaller (warning device) for electric transmission line icing (D, Mekhedints’, 1956). 


A pulse-frequency telemetry system for voltage, current, active and reactive power (M, Duma, M, Kalush- 
itsa, V. Roisman, M, Syrbu, 1956-1957), 


A controlled-current static device useful for summation (M. Duma, V, Roisman, 1957). 


A device for tele-dynamometry of petroleum equipment at great depths (A, Andrei, K, Nikulesku, 1957- 
1958). 


An adding device for pulse-frequency systems (M, Duma, M, Syrbu, V. Roisman, 1957-1958), 


Work began in 1958 on the development of remote control devices built of static elements (M, Syrbu, M, 
Duma, M. Pauker) and, in 1959, work began on the development of devices for the telemetry of non-electric 
quantities, 


Below are considered in more detail certain questions relating to the development of frequency ~pulse sys~ 
tems, adding devices and a device fortele-dynamometry, 


As a result of the analysis of the most widely used modern telemetry systems for remote operation, the pulse- 
frequency system was adopted, as the system giving great reliability, simplicity and economy of apparatus, and pos- 
sessing a number of known advantages, 
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Fig. 1 


In the transmitting devices, pulse generation is implemen- 
ted by a high-frequency generator (100 kilocycles, the mutual in- 
| duction between the loops of which is varied by being passed 
l through the notched disks of a transducer with an electric meter 
as the basic portion of its design. Periodic variation of the mutual 
rl induction leads to the periodic excitation and quenching of high- 
tc —, frequency oscillations; after rectification, pulses are obtained. The 
electronic receiving device accepts the pulses applied to the com- 
munication channel and transforms their frequency to a constant 
current which is applied to an indicating instrument. Separate trans- 
ducer types were developed for the telemetry of each electrical 
quantity (voltage, current, active and reactive power; the latter for 
— — 3 the cases of unitary and binary direction of overcurrent power), 


The transmitting and receiving devices do not depend on the quan- 
tities, subject to the telemetering; the indicating instruments dif- 
Fig, 2 fer only in the scale graduations. 


So far as we know, the literature on telemetry does not give 
sufficient answers to the essential question relating to the design of transducers based on meters, To increase the 
accuracy of meter-transducers, it is necessary that the rotational speed of their disks, and consequently, the pulse 
frequency, be as low as possible, while to guarantee accuracy of reception and rapidity of telemetering it is 
necessary that the pulse frequency be as high as possible. As is well known, this contradiction is usually solved 
technically by using photo-electric devices, permitting a higher number of pulses per second (above 20) to be 
obtained for a low nominal angular disk speed, The solution just described has many disadvantages: the neces- 
sity of creating special meter designs, the presence of an optical system, low reliability, etc. 


The other contactless transducers so far created,which are based on the modulation of mutual induction,are, 
as we know them, instruments with low pulse frequencies, 


In the apparatus developed at ETNII, there is a transducer which provides operational oscillations for the 
high-frequency generator by pulses of increased frequency, reaching 24 cycles, for low nominal disk rotation 
speeds, The frequency of 24 cycles was chosen on the basis of the operational conditions of the communication 
channel; frequencies of 30 to 40 cycles can be attained. The increase of the frequency was obtained by using 
a special, very simple, system of “magnetic optics". The principle of operation of such a system is that a mag- 
netic flux is compressed in the apertures of three baffles placed one over the other; the meter disk is the middle 
one of these baffles, and the moving one. Thus, the apertures are successfully made narrow and close to one 
another; with this, the number of apertures is increased. 


Meters are not ordinarily considered high-accuracy instruments, but investigation showed that, if they are 
used as transducers in telemetry, they are found under very favorable conditions, Thanks to the removal of the 
counting mechanism from the meter, and also to the introduction of an initial velocity, the frictional errors are 
decreased. Moreover, in using a meter as a power transducer, the error resulting from the measurement of power 
lower than the nominal, must be related to the nominal value of power, as is conventionally done with any in- 
dicating instrument, The solid line on Fig. 1 represents the error curve,related to the measured value of power, 
while the dashed line relates it to the nominal value; the axis of abscissas gives power values related to the nom- 
inal value, 


The curves (dashed), giving the error may, in a number of cases, be taken into account by the scale gradu- 
ation, The telemetry error will be then determined only by the instability (which is insignificant) of these curves, 
The complete role of the error curves is only retained in summation. In telemetry by call, when it is necessary 
to provide interchangeability of transducers, certain errors will appear due to the non-identity of the error curves 
for the various transducers, 


Thus, the basic error of a power transducer, based on a class 2.5 meter, is of the order of 0.5% for continuous 
telemetry and is lower than 1% for telemetry by call. The additional errors of such a transducer, when related 
to the nominal value of power, are also partially decreased in comparison with the meter errors. This pertains, 
for example, to the additional errors for cos g < 1. 


For telemetry of voltage and current, the dependence of the pulse 


R frequency on the value of the quantity being telemetered was taken to 
be a square law, since this allows the required accuracy to be obtained 
together with maximum simplicityof transducer construction. 
” Although it is ordinarily considered necessary to widen the upper 
The portion of the scale only in the case of a tele-voltmeter, with a linear 
= scale being recommended for a tele-ammeter, in this case a square scale 
was used. This was done because the tele-ammeters were set up for in- 
ene- spection of the maximum permissible currents in the line transformers; 
in practice, when there are tele -wattmeters and tele-voltmeters, the 
fee lower portion of the tele-ammeter scale is not used, 
Under these conditions, a convenient solution is to build transducers 
2 Fig. 2 with induction screening systems, The transducers are built by the pro- 


per modification ci single-phase standard meters; the measuring wind- 
ing current or voltage is mounted on the core of the meter's voltage wind- 
ing. By using the system, shown on Fig, 2, of short-circuited turns, a travelling magnetic field is obtained, this 


“ts field consisting of four magnetic fluxes, ,, &,, 3, and ®,, shifted with respect to each other, as is obvious from 

Ise the vector diagram given in Fig. 3. The core of the meter's current winding is replaced by a three-leg core with 
the same system of short-circuited turns. The voltage for the initial disk rotational velocity is supplied to the 

i winding of this core. In order to eliminate interaction of the two fields, the voltage for the initial velocity has 
a frequency of 100 cycles, which is obtained by the use of a magnetic doubler of the supply block's frequency. 

$- The accuracy of the transducers just described is high compared to that of the single-phase meters, This is de- 
termined by the conditions noted above for the case of power transducers, In addition, certainly, there are no ad- 
ditional errors for cos g < 1 for tele-voltmeters and tele-ammeters. 

are, 


re 

wer Fig. 4 

i- 

ty The circuit developed for the transmitting device is shown in Fig. 4, The tuned circuits A and B are posi- 

om tioned in the transducer on the two sides of the notched disk. Unsoldering is provided for at the output of the 
high-frequency transformer, to allow for post-rectification pulse transmission at the different levels necessary 

adu- under the conduction or frequency conditions of the communication channel, 

htex The design is very reliable; self-excitation and quenching of high-frequency oscillations occur for large 

wil variations in source voltage or circuit parameters, 

rves 


Tuning of the oscillatory loops to one and the same frequency (great accuracy in this is not necessary) pro- 
vides interchangeability between transducers and transmitting devices, so that telemetry by call can be imple- 


— mented by only one transmitting device by switching two high-frequency transducer outputs, 
ns, An investigation of the transmitting devices for pulse-frequency devices led to the conclusion that it was 


advisable to contruct static pulse-frequency apparatus, This is a very promising direction, particularly for the 
telemetering of non-electric quantities, and also for those electrical quantities which are easily transformed to 
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constant currents, but are subjected with difficulty to transformation to rotary motions (frequency, phase , total 
power). For telemetering of current, voltage and, particularly, power, it seems to us that instruments using rotary 
motion, as for example, the type described above, still remain simpler, more reliable, cheaper, more accurate 
and faster than static ones, 


The receiving device is a pulse frequency meter, based on the prevalent principle of measuring the average 
current of a charge -discharge capacitor, switched with the frequency to be measured. 


Fig. 5 


The circuit shown in Fig, 5 contains certain stages which are well-known, and are used in other similar 
circuits, We note certain details of the circuit, 


The use of a flip-flop state as an amplifier-limiter excludes the possibility of errors appearing for weak 
signals down to reception cut-off. Final pulse shaping is implemented in the output stage by limiting grid cur- 
rent, on the one hand, and by cutting off the tubes on the other. The lower pulse threshold,when type 6P6S or 
6P3S tetrodes are used in the output stage, does not exceed 10 volts,(even when the screen voltage is lowered to 
70 volts), is quite stable under oscillations in heating, and differs little from one realization of the same tube type 
to another, It is also necessary to mention that any change of this voltage must be referred, not to the value of 
the voltage being considered, but to the voltage drop across the charging condenser, which is of the order of 200 
volts. 


In the receiving device, a convenient method of verifying the fundamental calibrating quantity — the volt- 
age drop across the charging condenser — is provided. By pressing a telephone key in the receiving circuit, the 
pulses from the communication channel are disconnected, but the output of the indicating device is switched, 
via a shunt and an additional impedance, to the measurement of the voltage drop in charging. By means of a 
potentiometer, one can adjust the calibrating divisions on the instrument scale, with the purpose of providing 
interchangeability of the receiving devices at a definite pulse frequency and a definite portion of the scale, By 
means of another potentiometer, one can compensate for the differences, arising in instrument manufacture, 
between the actual and the calculated values of shielding condenser capacity and, by testing during the process 
of using the instrument, for the variation of the shielding condenser capacity due to aging. 


Let us now turn to the question of smoothing the output current of the receiving device. As is well known, 
this smoothing is limited by the required speed of telemetry and also by the accuracy and linearity of the scale, 
Meanwhile, smoothing is determined for the very lowest pulse frequency, As the pulse frequency increases, and 
consequently, the output current, the impedance of the core drops sharply and the time constant is decreased, 
With this, the smoothing conditions are met at the beginning of the scale and the response speed conditions are 
met for the remaining portions of the scale. 


The accuracy of the entire telemetered path (transducer, transmitting device, receiving device, indicating 
instrument) influenced the level of accuracy of well-telemetered systems. The fundamental errors do not exceed 
1.5% for insignificant addition errors, 


All the electronic stages are stable, the tubes were underloaded and slightly underheated, which provided 
a long service life. A supply block was developed from which all the necessary voltages for from one to four 
transducers and transmitting devices, or from one to four receiving devices, can be supplied. 
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The apparatus is mounted, together with the corresponding supply block, on a common stand, All electronic 
blocks are connected by means of 14-contact cannon plugs. With calibration, there is provided interchangeabil- 
ity for rapid replacement of defective blocks by spares. 


More than two years of testing verified the design characteristics of the apparatus developed, which has its 
first mass-production run this year, 


For the development of summing devices of a given accuracy, the aim was set of implementing a system 
with the least number of different instruments or, at least, with the least number of different theoretical solutions, 
to be used for inputs to the system's blocks and devices. In addition, it was decided to avoid, if possible, the use 
of specially built mechanisms, and mechanisms requiring complicated adjustments and regulation, 


In the consideration of the possible variations of power summation, such solutions were discarded as adding 
the secondary currents in current measuring transformer, summing the mmf's in adding transformers, or adding 
torques in special multi-disk meters with large numbers of measuring systems, The special features of electronic 
pulse-frequency receiving devices with rectifying bridges with grounded supply diagonal points complicate sum- 
mation on impedances, In addition, the solutions listed are partial ones, allowing summation either at the trans- 
mitting or the receiving points, with no possibility of relaying in the latter case, 


It was agreed to use summation of constant currents in devices which did not require galvanic connections 
between the circuits; such devices could only be magnetic amplifiers. 


To use analogous summing devices at the transmitting and receiving points, we developed a dc power trans- 
former, by means of which we simultaneously solved the problem of power telemetry of conductor communication 
channels, 


The adder was developed in two variants, The first variant was designed for transmitting constant currents 
for local reporting or for telemetry over conductor communication channels, and the second variant was designed 
for transmission by pulse-frequency systems, 


A transformer of power to dc was constructed on the basis of a standard three-phase induction meter with 
three systems from which the counting mechanism was removed. The active moment engendered by two of the 
induction systems was equalized by a contra-acting moment from the third system. By means of a two-tube 
circuit, the contra-acting moment is automatically regulated, providing an equilibrium state. Simultaneously, 
a dc current appears at the output which is strictly proportional to the contra-acting moment, 


The summer for transmission, along conductor communication channels, is a universal static balance tele- 
metering device (with controlled curgent), The communication channel is part of a circuit containing a total 
negative feedback path which shunts the entire scheme, 


The summer for pulse-frequency systems provides accurate transformation of dc sums into disk rotations 
of transducers constructed on the basis of single-phase meters; pulse generation is implemented by the same 
methods used in other pulse-frequency instruments, In making up the scheme, stages were used which are em- 
ployed in transformers of power to dc current and in static telemetering devices with controlled current. Thus, 
there are summing magnetic amplifiers in the scheme null-organs and electronic amplifier stages from control- 
led current devices and the electronic output stage from transformers of power to dc current, By means of the 
output stage, control of the transducer's induction system is obtained and, simultaneously, a dc current is ob- 
tained for the negative feedback path. 


All the summing devices have the same physical construction as in the pulse-frequency telemetering sys- 
tem, and are supplied from the same supply block, 


The device for tele-dynamometry of petroleum equipment at great depths provides the registering and 
transmission, by call, over a distance of diagrams which give the functional relationship between the stresses 
in the burnished piston and the translation of the piston, 


The diagrams are transmitted by wire over distances up to 3 kilometers to central dispatcher points, The 
device can service a large number of stands on a call basis. 


For the registering of the diagrams, measurements of two parameters, stress and translation, are simultane- 
ously transmitted over a distance, Both quantities are transformed into electrical signals by the corresponding 
transformers. 
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A specially developed stress transformer of the magneto-elastic type has high sensitivity, and contains no 
moving parts, which is particularly important in conditions of work where there is a corrosive medium. The trans- 
former core is constructed from sheets of material with sharply defined magneto-elastic properties. The trans- 
former is one of the arms of a bridge; another, neighboring, arm consists of an identical transformer on which 
the stresses do not act. This arm is used for compensation of phenomena which arise due to temperature vari- 
ations, The remaining two arms are active impedances, The bridge is supplied with voltage at a frequency of 
5000 cycles from a generator at the dispatcher point. The supply voltage is so regulated that the variations of 
the transformer's permeance occur on the linear portion of the curve expressing the permeance as a function of 
current, The bridge's unbalanced voltage, proportional to the stress in the piston, is amplified, rectified and trans- 
mitted in the form of a de voltage along the communication channel, At the dispatcher point it is again ampli- 
fied and applied to the vertical plate of a cathode-ray tube with a long afterglow. 


The position transformer is implemented in the form of a potentiometer which is supplied from a stabilized 
50-cycle voltage source, The voltage tapped off between one end of the potentiometer and the slider, which 
moves with the piston's controlling device, is proportional to the translation of the piston. The voltage is ap- 
plied to the communication channel, rectified in a voltage rectifier-doubler, and applied to the horizontal plate 
of the same cathode-ray tube. 


Under the action of the two voltages, which reproduce the stress and the translation, the cathode-ray tube 
presents the work-cycle diagram of the pumping equipment. 


The capability is provided, at the dispatcher point, of so regulating the amplification, by applying calib- 
rating voltages at the input, that an accurate correspondence may be established between the deviations of the 
electron beam on the screen and the measured stresses. 


Currents of three different types were used: for supplying the bridge, for the stress signal and for the trans- 
lation signal, This allowed interaction (cross-talk) in the communication channel to be avoided. 


A block schematic of the device is shown in Fig. 6. 


Fig. 6. 1) Transformer; 2) bridge, 3) amplifier-recti- 

fier; 4) dc amplifier; 5) the translation transformer 

6) rectifier. 


There is also registered, at the dispatcher point, the variations in current of the electric motor of any 
pumping device, the choice of which is made by call, For this, tele-ammeters of a rectified system are used, 
the constant currents of which are applied to the cathode-ray tube's vertical plate; as in the registering of the 
stress diagrams, the signals from the translation transducer are applied to the horizontal plate. 


The device for tele-dynamometry and electro-diagramming was subjected to testing for years, with satis- 
factory results. 
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A CONTROLLED -CURRENT STATIC TELEMETERING DEVICE 


M. Duma and V. Roisman 


(Bucharest) 


The theory of operation is given, and the circuit is described, for a controlled- 
current static telemetering device, developed at the Electrical Engineering Research { 
Institute of the Ministry of Heavy Industry of the Rumanian People's Republic, The 
operation of the scheme is based on the properties of an amplifier with a total feed- 
back path. An analysis of the action of the scheme, and its individual parts, is given 
and the results of testing a specimen static device are provided. 


For the telemetering of a number of electrical and, particularly, non-electrical quantities, the output power 
of the primary transformer is insufficient to allow _ its direct use as the transmitting device in a conductor tele- 
metry channel. For outputting into the communication channel, widespread usage is made of intermediate ampli- 
fiers (balance, compensation systems, or controlled-current systems), ordinarily constructed in accordance with 
the principle of automatic equalization of moments in duplex electromeasuring systems [1]. 


For satisfactory transformation accuracy, such instruments have a number of disadvantages, related to the 
presence of moving mechanical parts; one may also cite the difficulties attendant upon the necessity for summa- 
tion, etc, 


In the Electrical Engineering Research Institute of the Ministry of Heavy Industry of the Rumanian People's 
Republic there has been developed a controlled-current static device based on the properties of amplifiers with 
total negative feedback; The theoretical difference between the existing balance instruments and the static ones 
consists of the following. In balance instruments, within the limits of system sensitivity, the input and output 
quantities always equal each other (certainly, in reduced magnitude, in the steady state, and ignoring the twisting 
torques of the current-carrying springs), In static amplifiers with negative feedback, there will always exist a . 
definite difference between these two quantities which, in its turn, is proportional to the input quantity and tends 
to zero only for infinite gain. 


In practice there is no necessity to have infinite gain, The gain in an open-loop circuit is defined only by 
the necessary stability of transformation under the possible disturbing stimuli, and the statism thus obtained, is 
taken into account by graduation, In the following exposition, statism will be ignored, 


The comparison element (the null-organ)and the amplifier stage with feedback can be of different types. 
A rational choice is the conjunction of electronic and magnetic amplifiers to implement the tracking (following 
function of the null-organ [2]. ‘ 


A de from the primary transformer is applied to the magnetic amplifier's input, The magnetic amplifier's 
output voltage is amplified and rectified by an electron circuit whose dc output is both the current sent along the 
communication channel and the feedback current passing through the corresponding winding of the magnetic 
amplifier, Subtraction of two currents can also be carried out on a resistor, with only the difference current being 
applied to the magnetic amplifier's control winding [3} in the sequel, this case will not be considered. 


Thus, the circuit will always maintain such a value of output current which would correspond to the input, 
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Controlled-current static instruments possess the following advantages: 
a) In them there are no moving parts subject to wear; misalignment, and giving rise to frictional errors, 
b) The instruments are constructed of standard, mass-produced elements. 


c) Use of the instruments is quite simple; it is appropriate to note that it is not necessary to provide any 
regulating organ for use, 


d) A high dielectric stability between the input and output circuits can easily be implemented, this being 
necessary, for example, in the telemetry of dc voltages from tractive stands. 


e) For telemetry by call, the necessary scale transformations can be easily carried out by switching the 
number of turns in the magnetic amplifier's control winding. 


f) Summation is easily carried out by supplying the proper number of control windings for the magnetic 
amplifier. With this, the summing circuits are insulated from each other, there is no difficulty in obtaining an 
accurately weighed sum and the number of summands can be very large; the accuracy is not decreased, if one 
of the summands equals zero, as occurs when many other summation methods are used. 


The electrical scheme of the device developed is shown in Fig. 1. 


The magnetic amplifier is executed in the form of a bridge-type amplifier with toroidal Permalloy cores. 
During the transient response, the difference between the control and feedback ampere-turns may change sign. 
In order to avoid self-excitation of the circuit due to transition of negative feedback to positive, the phase of 
the magnetic amplifier's ac output voltage is changed when the sign of the residual control ampere-turns changes; 
the electronic amplifier reacts only to the correct phase. The magnetic amplifier is provided with a winding for 
the initial reading current, which is used to check the state of the device when no control signal is present, 


The four-stage electronic block is constructed of two 6N9S and 6N7S tubes, The first two stages amplify 
the magnetic amplifier's ac output voltage. The third stage is simultaneously an amplifier and a phase-sensitive 
rectifier, which control the operation of the output stage by varying the magnitude of the negative bias. The 
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be 
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ac anode voltage of the phase-sensitive stage is in phase with the magnetic amplifier's ac voltage. We note that 
the rectified voltage of the third stage, which is negative with respect to ground, eliminates the necessity of a 
separate source of biasing for cutting off the output stage. 


The de anode current from the output stage is applied to the communication channel and passes through 
the feedback winding and through other instruments for local recording. The communication channel is grounded, 
but (and this applies to all existing controlled-current instruments) it is necessary to provide protection against 
the increased voltage when an arm breaks. 


It should be mentioned that errors can appear if, in the magnetic amplifier's control winding, there are ac 
components of the current with frequencies which are even harmonics of the magnetic amplifier's supply voltage. 
The source of such ac components can be either the primary transformer or the magnetic amplifier itself, In the 
first case, the error depends on the amplitude and phase of the ac component; itis possible to neglect these for 
pulsations of less than 2%. As for the second case, it is necessary to provide an output impedance of the primary 
transformer's ac current which is higher than some critical value, or to graduate the device from a source which 
has an internal impedance of the same order as the primary transformer's impedance, 


The engineering data on the instrument are as follows: the output current is 5 milliamperes; the fundamen- 
tal error is the error of the indicating instrument; deviations of the characteristics from linearity do not exceed 
0.5%, and comprise the basis error in summation at the receiving point; with an unstabilized supply, a variation 
of + 15% in the supply voltage engenders an additional error which is less than the error due to variations in the 
initial indication current; as the impedance of the communication channel varies from zero up to 5 kilohms, 
the output current changes by less than 0.5%; the time to establish the output current does not exceed 2 seconds. 


‘ 


Fig. 2 
An overall view of the instrument is shown in Fig. 2. 
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LETTER TO THE EDITOR 


In [1] there was given the equation for the mechanical characteristic of an asynchronous machine with a 
cylindrical and massive rotor in the form: 


2+1 030_ 
+ (1) 
1,030 


At the same time, there was given in [2] (without a derivation) the equation, similar to (1), for the mecha- 
nical characteristic of an asynchronous clutch with massive armatures: 


24+ V2 (2) 


The fact that the numerical coefficients, 1.030 and ¥2, in Equations (1) and (2) are not identical requires 
further explanation, which is given below. 


Equation (1) was derived with the assumptions that the active impedance of the machine's stator was neglig- 
ibly small, i.e., ry = 0, and that the ratio of the reactive and active components of the equivalent, fully reduced 
massive rotor impedance was constant for all values of slipping, equalling x‘,/r', = 0.60. 


It can be shown that the equation of the mechanical characteristic of an asynchronous machine with a mas- 
sive cylindrical rotor, derived by the method of [1] under more general conditions, namely for r, # 0, and 
x',/r'y = a, has a following form: 


2(1+- 9) 


a+b 


where» 


Here, the plus sign refers to the motor, and the minus sign to the generator mode of machine operation, 
It is not difficult to obtain Equations (1) and (2) as special cases of expression (3). 


It was already mentioned that Equation (1) is obtained from the condition that r, = 0; i.e., for b = 0 and 
x",/r', = a = 0.60. 


Formula (2) is easily obtained from expression (3) with the conditions that r, = 0 and a= 1, 


It is well known that the assumption that a = lissues from the simplified theory of massive rotors in which 
account is not taken of the dependence of the rotor material's magnetic permeability on the magnetic field 
strength, which varies in the direction of permeation of the electromagnetic field, or of the hysteresis losses in the 
rotor’s steel. 
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At the same time, the validity of the assumption that a = 0.60 
M was proven theoretically and experimentally by L. R, Neiman [3] 
for all ferromagnetic materials used in practice, when the follow- 
ing conditions are met: 


\\ 


a) The ferromagnetic body (in the given cases, the clutch 

2L— “y armature or the asynchronous machine's rotor) must be massive, 

Zz i.e,, its cross-sectional diameter, and the least linear dimension 

of this section, must significantly exceed the depth of permeation 
4 of the electromagnetic waves, 


1, 


b) The case of a “strong” magnetic field must obtain, where - 
in the magnetic field strength at the surface of the massive ferro- 
Lis a magnetic body exceeds that value of field-strength for which the 


magnetic permeability is a maximum, 


It should be mentioned that the value of the coefficient a = 
= 0.60 is used for analogous computations by other authors in cer- 
tain works devoted to electromagnetic slipping clutches with mas-~- 
sive armatures, for example, in the work of V. S, Sharov [4]. The experimental data given in this work show that 
the magnitude of the coefficient a depends on the slipping, particularly for small slipping, but practical calcula~- 
tions show that this dependency can be ignored and that, for all values of slipping, one can take a = 0,60, 


0 id 


The following circumstances has a more essential influence on the mechanical characteristics defined by 
Equation (3). 


The massive armature (rotor) parameters which enter into the formulas for s,, [5] depend on the magnitude 
of the equivalent reduced current I', of the massive armature (rotor). The current I's, as is well known, is a 
function of the slipping s. It follows from this that to each value of s there will correspond a definite value of 8,,, 


1.€., Sor = f (8). 


It was shown in [1] that, by means of the parabolic approximation suggested by L, R, Neiman [3] for the 
fundamental magnetization curve (for strong fields) for materials from which the massive armatures (rotors) were 
prepared, it is possible totake into account, in a comparatively simple manner, the relationship scr = f(s) and, 
by using Equation (3), to obtain the universal relative mechanical characteristic M/ M,, = £ (8/8,,), which is ap- 
plicable to all asynchronous clutches with massive armatures and also for all asynchronous massive rotors (with 
= 0). 


On the figure there is shown the mechanical characteristic, 3, of an asynchronous machine, computed from 
formula (3) for a = 0.60 and with account being taken of the function s,, = f(s). The characteristic is computed 
for the machine used as an example in [6] (the critical slipping of this machine is found in the breaking mode 
of operation, i.e.,s_ > 1.0). For purposes of comparison, the figure also shows the mechanical characteristic, 
2, of this machine, calculated only on the basis of formula (2), Comparison of the computed characteristics 
with experimental results, shown by curve 1, (from the data in [6]) indicates that the use of formula (3), with 
account taken of the function scr = f(s), gives the more accurate results, 


We note in conclusion, that formulas (1)-(3) are particular cases of the generalized equations of the relative 
mechanical. characterisitc of an asynchronous machine [5}: 


2(1+ 9) 


where m_ is the exponent of the slipping in the formulas for the dependence of the parameters of the ae 
circuit of the replacement circuit for the asynchronous machine on the slipping, for example, r'g = r'g9/ s™ (for 


asynchronous machines with constant parameters, m = 1.0; for asynchronous machines with massive rotors, m = 
= 0.5; for hysteresis machines, m = 0), 


(4) 


V. M, Kutsevalov 
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CONFERENCE ON INVARIANCE THEORY AND ITS APPLICATIONS 
TO AUTOMATIC DEVICES 


From October 16 to 20, in Kiev, there was held a conference on invariance theory and its applications to 
automatic devices, The conference was organized by the Kiev City Seminar on Automatic Control Theory and 
by the Electrical Engineering Institute of the AN USSR at the initiative of the Office of the Technical Science 
Department of the AN USSR, 


Taking part in the work of the conference were more than 200 delegates from 38 research organizations, 
educationalinstitutions, and industrial enterprises in Moscow, Kiev, L'vov and other cities. At the conference, 
there were heard, and discussed, 25 papers given by representatives of 12 organizations, 


In opening the conference, A. Yu, Ishlinskii, director of the Kiev City Seminar and Academician of the 
AN USSR, stated that the idea of convening a conference,which would be devoted to a consideration of the theory 
and practice of automatic devices using partial, or complete, compensation of external disturbances, arose in 
the spring of 1958 at an enlarged meeting of the seminar in connection with an interesting paper of Academician 
V. S. Kulebakin on the K(D)-transformation proposed by him, and its applications to the analysis and synthesis 
of linear systems, 


The questions dealt with at the conference have a long-established history, and one can adduce a number 
of examples of automatic devices, starting with windmills and ending with modern synchronous drives, where 
the idea of compensating an external load is used to one degree or another, 


Certain theoretical developments in this direction were not always irreproachable, and were occasionally 
subject to bitter criticism, From this derived the attempts of the late G, V. Shchipanov to formulate the condi- 
tions for complete compensation of external disturbances in automatic control systems which act on the deviations 
of the controlled quantities. With this, even the criticisms themselves turned out to be insufficiently guarded, as- 
serting, for example, that in no actual physical system was complete compensation of external disturbances pos- 
sible, 


As always,happens in such cases, concrete examples turned out to be more useful than general considerations, 
particularly examples taken from practice of actual operating devices. 


One should distinguish at least two forms of implementation of the independence of the variations of one of 
the parameters characterizing the system from the law of variation of one of the external loads acting on it. 


To the first form belongs the decomposition of the set of differential equations ,describing the system's be- 
havior into two or several sets: the parameters of interest fall into one of these sets and the corresponding loads 
into another, With this, the processes which correspond toeach subset of equations occur, generally speaking, in 
the same elements of the system being considered; the conditions for decomposition are, essentially, the so-called 
conditions of invariance, 


The other form of implementation of the independence of a system parameter from the load consists in 
using additional information as to the variations of the load itself and, consequently, providing for two indepen- 
dent channels of system stimulation, 


In certain cases, complete compensation can ensue, if there is information to the effect that the load is the 
sum of harmonic or exponential functions with given parameters, 


In conclusion, A. Yu. Ishlinskii mentioned that the problems of automation have tremendous significance 
for our country. This was asserted by the resolutions of the twentieth session of the KPSS, The questions of 
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automation occupied a significant portion of the paper of N, S, Khrushchev at the special twenty-first session of 
the KPSS, It is therefore particularly important that all branches of the theory and practice of automatic devices 
be developed to the fullest degree and that further growth of domestic technology be made possible. 


Up till now, the theory of invariance and external disturbance compensation has received inadequate at- 
tention, In addition, the practical value of these questions in matters of creating perfected automatic devices 
is quite high, Therefore, the tasks of this conference are to supply the results of the many works of Soviet scien- 
tists in this direction, and to discuss the further paths of scientific and practical investigation. 


The first paper at the conference, with the subject, “Invariance as a fundamental principle of high-quality 
automatic regulation and control,” was given by Academician V. S, Kulebakin, The basic requirement im- 
posed on automatic control and regulation systems is that the deviation between the given and the true values of 
the quantity to be controlled, when acted upon by external disturbances, should be as small as possible, 


The regulation process is the better, the more fully compensated in the system are the effects of the ex- 
ternal disturbances, The physical essence of the control problem is the implemenation, to some degree or other, 
of compensation of the effect of disturbances on a selected coordinate and the attainment of a definite inde- 
pendence, or invariance, of it from the disturbance. 


The following forms of invariance are possible. 

a) Absolute invariance, 

b) Invariance to within <«, 

c) Selective or chosen invariance, 

d) Invariance tw within a free component. 

v) Invariance up to wie effec: of a derivative of the disturbing force. 


V.S, Kulebakin considered in detail the various enumerated forins of invariance and cited their domains 
of application, V,S. Kulebakin illustrated his paper by description of realized systems. 


In the paper, “The invariance problem for control systems based on deviations", A. I. Kukhtenko considered 
in detail the essence of the idea of compensation and the criticisms, in 1940-1941, of the opinions of G. V. 
Shchipanoy. The speaker went on to demonstrate the possibility of creating control systems based on deviations 
which satisfy the conditions of invariance to within € , having corroborated his conclusions by electro-simulation, 


B, N, Petrov, corresponding member of the AN SSSR, appeared before the conference with a paper on the 
subject, “On the criteria of applicability of the conditions for invariance in the synthesis of dynamic systems". 


In noting the positive value of invariance theory in the development of high-quality automatic control 
systems, the speaker. asserted that the principle of invariance can be applied and physically implemented if cer- 
tain definite conditions are met. 


B, N, Petrov formulated the following criterion: “The conditions of invariance may only be applied when 
the scheme contains at least two channels for transmitting stimuli between the point of application of the dis- 
turbing force and the point at which the quantity, in which we are interested, is measured". 


A number of examples were considered by way of illustrations of the criterion presented, 


The conference was opened for discussion of the papers of V. S, Kulebakin, A, I, Kukhtenko and B, N, Pet- 
rov, and in it participated, in addition to the speakers themselves, A, Yu. Ishlinskii, A. G, Ivakhnenko, E, P. 
Popov, G. S, Pospelov, O. M. Kryzhanovskii, N. M. Chumakov, V. N, Yavorskii, S, V. Krivonosov and others. 


In his concluding words of this dicussion, V. S. Kulebakin noted that none of the participants questioned 
the principle of invariance as such. This is an actual, physically realizable principle, In this connection, it is 
necessary to expand the theoretical and experimental work in this direction, Use of this principle does not mean 
that other principles will be discarded, On the contrary, it will be necessary to make use of them when advantage- 


The second day of the conference saw papers given by A. Yu. Ishlinskii, G, M. Ulanov, V. A. Bodner, G, S, 
Pospelov, A. G. Ivakhnenko and V, I, Nechiporenko, 
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In the paper, “Total compensation of disturbances arising in gyroscopic maneuvering systems,” A. Yu. 
Ishlinskii noted that when a variable disturbance acts on some mechanical or electrical system, the variations 
in the basic parameter of this latter depend, as a rule, not only on the current value of the disturbing quantity, 
but also on previous values of it, i.e., on the “disturbance's prehistory”. 


If the system's behavior is described by linear differential equations then, by artificial means, the influence 
of the “prehistory” can be eliminated, For this purpose, as is well known, one applies to the system's input, not 
only the basic disturbances, but also linear combinations of derivatives of the disturbances with respect to time, 
However, in the general case, these derivatives cannot, in theory, be reproduced accurately, and therefore such 
combinations can be implemented, as a rule, only approximately. 


In mechanical systems in which the disturbances are the inertial forces of system translational movements 
on the Earth's surface, complete compensation of the effect of these disturbances can be attained, Examples 
of such systems are the spatial gyroscopic compass of Gekkeler, Shuler-Bulgakov's physical pendulum, Baikov's 
integrator and certain other systems of the inertial navigation device type. 


In the paper of G.M, Ulanov on the subject, “Invariance to within € in combined non-linear automatic 
control systems with disturbances of bounded modulus,” there were formulated the conditions for small perturba~ 
bility of non-linear systems, In addition, the speaker formulated a criterion for dynamic accuracy of non-linear 
systems, and presented the results of a comparison of control systems based on deviations and combined systems, 


V. A. Bodner, in the paper "On the application of the principle of compensation to the construction of 
automatically stabilized distribution systems”, showed that disturbance of unidirectionality in systems carrying 
out the function of energy transformation is frequently the cause of disruption of stability and creation of auto- 
oscillations, For a definite class of systems, such a mode of operation is inadmissible and, therefore, must be 
eliminated by means of the techniques of automation. 


The dynamic properties of systems designed for energy transformation can be estimated on the basis of 
the properties of the quadripoles equivalent to them. The condition of system stability leads to the imposition 
of certain requirements on the coefficients of the quadripoles. These requirements may be met by connecting 
special compensating systems, constructed of passive elements, to the basic system, Their purpose reduces to 
the neutralization of the action of the internal positive feedback,which leads to instability. 


The speaker showed that the principle of compensation in complex distribution systems gives a greater 
stabilizing effect than the usual principle of stabilization based on deviation, 


G. S. Pospelov, in the paper"Autonomous and coupled regulation of an aircraft's parameters of motion,” 
presented the invariance conditions for the slipping angle of an aircraft maneuvering by course, An analysis of 
the longitudinal motion was used to explain the condition for invariance or autonomy of the control loops for the 
altitude and velocity of flight. 


In A. G, Ivakhnenko's paper, "Combined control as the general case of controlling state and quantity? 
there was given a general approach to combined systems, The speaker proposed renouncing the terms conven- 
tioanlly used today, “controlling quantity" and “basic disturbance,” and returning to the terms of classical auto- 
matic control theory. With the example of a concrete electrical system, A, G.Ivakhnenko showed the superiority 
of a combined control system as compared to other systems, In combined systems, tuning is rendered easier and 
its limits are extended and, in addition, the capabilities of using the conditions of invariance for improving sys- 
tem dynamics are significantly widened, In them, the coefficient of statism may assume a negative value with- 
out loss of static stability, which is used for improving the parallel operation of the object of regulation, In con- 


' clusion, the speaker dealt in detail with one of the varieties of invariance, suggested and named by him as a 


forcing device, which is used in program~-controlled systems, 


In the paper, “Invariance to within € in linear two-loop combined automatic control systems," V.1, Nechi- 
porenko showed that the idea of invariance to withine can be extended to the case of multiloop systems, in par~ 
ticular, to two-loop systems, 


In the paper were considered the conditions for absolute invariance, the approximate conditions for invari- 
ance as applied to slowly varying processes and invariance to within € on finite intervals of time, 


Conditions were given for broadening the domain of realization of the idea of invariance in combined two- 
loop automatic control systems, 
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Taking part in the discussion of the papers were N, I, Simonov, V. A, Besekerskii, P. I. Kuznetsov, L. V. 
Tsukernik, I, M, Krikheli, G, F. Zaitsev, V. G. Vasil’ev and others. The participants made a number of critical 
remarks anent the paper. 


Academician V. S. Kulebakin, who took part in the discussion, made a number of concrete remarks about 
the papers and gave an analysis of the foregoing discussion, 


In summing up the results of the two-day's work, conference chairman A, Yu. Ishlinskii stated that the time 
had come to present the results achieved in invariance theory, not only in the form of scientific communications 
of some degree of complexity, but also in a scientific-popular form, in order to make invariance theory reach 
a wide circle of scientific workers and engineers, A. Yu, Ishlinskii underlined the importance of obtaining infor- 
mation on disturbing factors. 


In continuing its work, the conference heard the paper of V. A. Besekerskii and S, M, Fedorov, "The design 
of combination control servo systems by the logarithmic frequency characteristic method", One of the methods 
of designing servo systems by using logarithmic frequency characteristics was presented in this paper, According 
to the author’s methodology, one must first determine the equivalent transfer function of the open-loop system 
with account taken of the regulation by the disturbing stimulus, and then construct the logarithmic frequency 
characteristic, Then, bearing in mind that the servo system's accuracy is affected by the form of the logarithmic 
frequency characteristic in the low frequency range, one proposes a compensating signal which takes only this 
range into account, 


For the standard forms of open loop transfer functions in the low-frequency domain, atable is provided giv- 
ing the transitions to the transfer functions with account being taken of the compensating signals. 


The paper of P, I, Dekhtyarenko, "On the application of forcing devices to servo systems with constant- 
speed executive elements,” the application to servo system was explained for forcing devices constructed with 
the use of detachable connections and permitting the approximate holding of the invariance conditions in systems 
with constant-speed executivemechanisms given step-wise excitations, The system has an independent biasing 
generator which is switched insimultaneously with the application of the disturbance, and provides a subsidiary 
signal, As a result, the executive mechanism is switched out with a lead such that the system proceeds to the 
required position in a free-coasting mode without repeated switching-in. 


In the paper, "Combined regulation of synchronous machine excitation on the basis of compounding and 
voltage correction," L, V. Tsukernik spoke of the development of systems realizing the principle of invariance 
with an accuracy up to the transient component of the quantity being controlled. The developed controllers have 
been put into industrial use. The speaker noted that the choice of the controller's block schematic and the intro- 
duction of additional compounding and correcting connections, significantly affect the operation and physical 
characteristics of the regulator, 


In G. K, Nechaev's paper, "A temperature regulator which reacts to variations in the disturbing stimulus," 
there was presented a method for designing regulators which use the condition of invariance for three working 
modes of operation, A regulator of such a type can be applied for special thermostats. 


V. N. Yavorskii gave a paper on the subject, "A new approach to the use of the invariance principle for 
designing compensated power servo drives," in which the author emphasized the role of differentiating reproducing 
and forcing devices in raising the quality of servo drive operation, and in which the author showed how the degree 
of compensation of a servo drive is determined. 


InV. M. Kuntsevich’s paper, "On the use of regulation by disturbance in extremal control systems,” it was 
shown that if, in extremal systems, one introduces additional connections for the disturbance and its derivatives, 
this significantly improves the control process, As an example, a system was adduced in which the absence of a 
compounding connection for the disturbance made it unworkable, 


In the paper of V, I, Kostyuk, “Correcting servo systems," it was shown that, if the circuit connections and 
the compounding connection coefficients are optimally chosen, the compounding being by basic disturbance, the 
corrector power will be minimal and the signal at the corrector input will vary in sign, This can be used for in- 
creasing the accuracy of combined systems by varying the corrector gain as a function of the input signal param- 
eters, 


In the paper, “An inductive-capacitive stabilizingdevice as an invariant system,” A. N. Milyakh and B, E. 
Kubyshin stated that on the basis of the principle of invariance, they had carried out a choice of parameters for 
a three-phase induction-capacitance transformer with antisymmetric reciprocity which possessed the properties 
of a gyrator, This device so transforms energy that the current at the output is proportional to the voltage at the 
input, and, conversely, the current at the input is proportional to the secondary voltage. The device developed 
is an invariant system in which significant disturbing forces, arising as the result of load impedance variations, 
only slightly affect the magnitude of the secondary current. 


The conference then considered questions in the theory of automatic control of coupled systems and the 
conditions of selective (autonomous) invariance as applied to such systems, 


In the paper of V, T. Morozovskii, "On the theory of mono-typic coupled automatic control systems with 
symmetric cross-coupling,” a method was considered for analyzing mono-typic multi-channel systems with sym- 
metric synchronized neutralizing coupling. It was shown that the neutralized system motion is invariant with res- 
pect to the signals passing through the synchronizing coupling. On the basis of the theory developed, an inves- 
tigation was carried out of the parallel operation of aircraft synchronous generators with automatic speed control, 
for frequency correction and for couplings equalizing the active load. 


The paper of L, V, Tsukernik was devoted to the investigation of stability of a symmetrically-coupled 
combined control system, Stability was considered in the "small", In a coupled dynamic system it is possible 
to isolate several mono-typic groups of complex objects, The motion in each group is described by equations 
in matrix form which may be replaced by operator equations in the generalized coordinate of the group and the 
generalized coordinates of all the remaining groups. Such a replacement allows Lyapunov's method to be used 
for estimating stability, In case of intra-group symmetry only two forms of system motion can be considered; 
group and intra-group, the stability conditions of which are invariant with respect to each other, The conclusions 
obtained were verified with the example of an energy system with several multi-aggregate stands, 


In the paper of Yu. G. Kornilov, "Autonomous regulation for parallel aggregate operation,” there was con- 
sidered I, N. Voznesenskii’s autonomy problem as applied to objects with one and with several loads under con- 
ditions of joint operation, The speaker showed that implementation of autonomous control for steam turbines 
with tapped-off steam allow one to attain, not only mutual independence of action of the speed regulator and the 
pressure regulator in the tap-off chambers, but also an independent distribution of the loads between the parallel- 
operating aggregates, which is important in practical usage. 


Other results regarding the effectiveness of autonomous control were obtained for steam boilers — objects with 
one load in the form of steam discharge, Autonomous systems had no explicit superiority over non-autonomous 
systems either in static or in dynamic modes of parallel boiler operation, By means of a concrete example, it 
was shown that, in this case, an autonomous system can give a qualitatively poorer transient response. 


A, I. Sud-Zlochevskii gave a paper on the subject, “on the application of the principle of invariance to 
systems with several controlled parameters”, In the speaker's opinion, systems with coupled control of several 
parameters, from the point of view of the disturbances acting on them, may be divided into two groups. In the 
first group belong those systems in which each disturbance has a direct effect only on *its” controlled parameter. In the 
second group are those systems in which each disturbance has a direct effect on all the controlled parameters, 


In systems of the first group it is impossible in principle to obtain complete invariance, To obtain invari-~ 
ance it is necessary to introduce couplings for the disturbances or for their derivatives, In systems of the second 
group, there is the possibility, in theory, of obtaining invariance of n-1 controlled parameters without introducing 
coupling by disturbance. However, a coupling providing invariance of n-1 controlled parameters is impossible 
to realize physically. 


In the paper of P. I, Chinaev, "On the question of synthesizing linear systems of coupled control,” three 
principles were considered which, in the speaker's opinion, can be considered fundamental to the entire synthesis. 
These are the principles of realizability, invariance and optimality. 


The author considers it reasonable to use the methods of matrix calculus for the synthesis of systems with 
many variables, since this permits the solution to be given compactly. The condition for realizability is decom- 
posed into two: mathematical and physical. By using equivalent initial conditions, establishing the connections 
between the system parameters and the initial values of velocity, acceleration, etc,, and by imposing limitations 
on speed and acceleration at the initial moment, one can obtain additional conditions for the synthesis, 
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The conditions for invariance (complete or selective) also impose limitations on the coefficients of the mat- 
rix equations, Finally, the optimality conditions are obtained as limitations imposed by the Euler equations on the 
integral relationships (functionals) which depend on the system parameters, The joint holding of these three con- 
ditions thus provides a basis for synthesizing systems with many parameters. 


The paper of N, G, Dyadyunov, "On certain questions of automatic flight control by means of discrete sig- 
nals,” was devoted to the choice of optimal speed, channel capacity and noise stability of discrete systems for con- 
trolling aircraft flight, The best results are given by combining a discrete system, operating with terrestrial in- 
formation, and an onboard extrapolating device. A self-correcting system is obtained. 


O. M. Kryzhanovskii gave a paper on “The quasi-invariance of the transient response in non-linear automa- 
tic control systems for mine hoists". 


The speaker introduced the new concept of quasi-invariant automatic control systems, assuming that this 
concept corresponds to the physical meaning of the technological realizability of a system. The speaker gave 
examples of systems which are almost invariant to individual disturbances, By a mathematical model, the various 
kinds of quasi-invariant transient responses in existing non-linear mine-lift devices were shown. 


In addition to the papers listed above, communications were made by: V. G, Vasil’ev, "On the invariance 
problem for linear reproducing systems," L. V. Karnyushin, "On the application of the principle of invariance in 
multi-motion electric drive systems” and G, F, Zeitsev, "On a broad-band differentiator”. 


After the papers and communications, there was a general discussion, participated in by Academician V. S. 
Kulebakin, V. A, Besekerskii, P. I. Kyznetsov, G. F. Zaitsev, V. P. Sigorskii, S. I. Totel"baum, A. I. Sud-Zlochev- 
skii and others. 


In the concluding portion, the participants discussed, and unanimously adopted, resolutions. 


N, A. Kachanova and P., I. Chinaey 


RESOLUTIONS OF THE CONFERENCE ON INVARIANCE THEORY 
AND ITS APPLICATION TO AUTOMATIC DEVICES 


(Kiev, October 16-20, 1958). 


Convened by the Engineering Department of the AN USSR, by the Kiey City Seminar on Automatic Control 
Theory and by the Electrical Engineering Institute of the AN USSR, the conference, as a result of its work, went 
on record as follows. 


1. Today, in the design and creation of automatic systems, both in the SSSR and beyond its borders an ever- 
greater use is being made of the condition of external disturbance compensation, or various forms of the principle 
of invariance, 


2. The condition of compensating external disturbances was first formulated in the SSSR by Prof. G. V. 
Shchipanov. This idea was the subject of mathematical investigations by Academician N. I. Luzin; by him, in 
conjunction with P, I, Kuznetsov, were formulated the criteria of absolute invariance and invariance to within €. 
Academician V. S, Kulebakin, and then A. G. Ivakhnenko, B. N, Petrov, A, Yu. Ishlinskii and others, proved the 
realizability of the principles of invariance in actual electrical and mechanical systems and automatic devices, 
in particular, in combined control systems, bridge circuits, gyroscopic systems, etc, 


3, It was by Soviet scientists that the theory of invariance was developed and automatic control systems 
were implemented on this basis. 


4, Papers published here and abroad, and also the papers and speeches of this conference's participants, at- 
test that the conditions of compensation and invariance have been further developed and sufficiently based on the 
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mathematical and physical sides; also developed has been a broad class of methods and technological means 
for compensating the effect of disturbing stimuli and for realizing the principle of invariance. 


5, It was established, based on the paper given in this conference, that the compensation of disturbances 
and the use of the principles of invariance allow the creation of automatic systems and various devices which 
are more highly perfected, from the point of view of quality of the control and regulation processes, stability, 
simplicity of implementation and reliability of operation, This is supported, not only by the theoretica! and 
experimental investigations, but also by the engineering developments of different automatic devices and measur- 
ing instruments, in wide use for more than 10 years in maritime matters, aviation, ordnance energy plants and 
other domains of technology, In particular, many boilers, turbines, generators and other apparatus of our heavy 


energy stations are regulated by combined methods; by disturbance or with the application of the principle 
of invariance, 


The conference considers it necessary to: 


1, Acknowledge the efficiency of and effectiveness of the methods of compensating disturbances and the 
principle of invariance for the development of, and application to, automatic systems. 


2. To recommend a broader development of the theoretical and, particularly, the experimental investiga- 
tions in the theory of automatic and other dynamic systems, acting on the deviations of the controlled quantity 
and on the disturbances (combined systems) and satisfying the conditions of complete or partial compensation. 


3. To acknowledge the desirability of presenting the theory of compensating the invariance, and also com- 
bined methods of control, in textbooks,handbooks and scientific-popular works on automation. 


4. To consider it essential to publish the works (papers and discussions) of this conference and, in connec- 
tion with this, to request the Presidium of the Academy of Science of the USSR to publish these works in the 
shortest possible time. 


5. To acknowledge the advantages of presenting papers on the theoretical and practical developments in 
automatic devices operating with additional stimuli from disturbances and on the basis of the principle of invari- 
ance at the forthcoming International Congress on Automatic Control to be held in 1960, 


6. To resolve, in the next few years, the following scientific -engineering problems: 


a) The development of invariance theory as applied to control and regulation with additional stimuli 
from external and forcing disturbances for linear systems with constant and with variable coefficients and lag- 
ging arguments, for non-linear systems, for systems with many controlled quantities, for self-adjusting and ex- 
tremal systems, and for systems with computing devices. 


b) The development of methods for analysis, synthesis and design of the systems cited above. 


c) The search for new, better methods of measuring disturbing stimuli and the development of the corres- 
ponding instruments, measuring disturbances and loads of various characters, 


d) The measurement of disturbances and loads in various productive objects and the compiling of statis- 
tical data on them. 


e) the sharpening of the terminology in the domain of the various methods and principles of regulation 
and control. 


1. In view of the fact that the decisions of the Presidium of the Academy of Science of the SSSR on April 1, 
1941, on the conclusions of the special commission on the work of Prof. G. V. Shchipanov delayed the develop- 
ment of the theory and practice of constructing automatic systems using compensation of disturbing stimuli and 
the various principles of invariance, and does not permit the proper progress today in this important branch of 
automation, we ask the Presidium of the Academy of Science of the SSSR to review its decisions regarding the 
erroneous conclusions of the afore-mentioned commission, 


8. The resolutions of this conference should be published in the periodic press (in academic and other 
scientific-engineering journals), 
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